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HOW NOT TO FIND AN OIL FIELD* 


JOHN SLOATt 


ABSTRACT 


This paper discusses causes and cases of seismograph misinterpretation, with ex- 
amples of errors resulting from certain spreads, from weathering computations, and 
from incorrect velocity assumptions. The question is raised as to whether subterranean 
weathering corrections, echo reflections, transverse waves or regional velocity gradients 
have caused serious mis-mapping of structure. The evaluation of results from dip shoot- 
ing is reviewed, with cases showing how oil fields have been missed. 

The psychological factor in picking records is stressed and the use of marked records 
in correlation is questioned. The need for a better criterion of correlation is expressed, 
and it is suggested that a more complete study of results might result in fewer dry holes. 


We geophysicists like to speak with pride of the correct correlations 
we have made and the oil fields we have found. We gloat with glee 
when the drill reveals a fault we have predicted, but are prone to 
forget the cases where the drill disproves a fault, or where the drill 
says up where we said down. We forget the case where we had a flat 
dip or two, failed to check them, and Company X found an oil field 
(this, of course, was the fault of the geologists—they didn’t want to 
spend the money to shoot continuous spreads), or the case where we 
showed a little closure from incomplete data or study of the records, 
and the drill found the “high” barren (this, too, was the geologists’ 
fault—they insisted on using that erratic dip to make closure). Then, 
too, how about the case where we ran a line between those two wells and 
put ina fault to tie, and it later turned out that therewas no fault but 
that our velocity chart was off a thousand feet at one of the wells? 

You may say that it is human nature to forgive and forget our 


* Read by title at the fourteenth annual meeting, Dallas, Texas, March 1944. 
t Union Oil Company of California, Bakersfield, California. 
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mistakes. Some may contend that otherwise the conscientious inter- 
preter will develop.a defeatist attitude and be afraid to put down a 
doubtful correlation which may prove to be correct. To some extent 
this may be true, but probably very few geophysicists have a defeatist 
complex. Most of us are conceited enough to think we know just a little 
more than the next fellow about mapping structure. It might do us all 
a lot of good, then, if we would look into our mistakes (whether we 
admit them or not) and search for the answers to our errors. In so 
doing, perhaps next time we can at least avoid the same pitfall, and 
rely on a fresh alibi to explain any new misinterpretations. 

The following list of mistakes which many of us have undoubtedly 
made with the seismograph is not complete, but may encourage the 
revelation of a few others and help us all to see why they were made. 
Miscomputations and other run-of-the-mill errors are not considered, 
as it is assumed the computers were awake that day and that the 
party chief was sober. The mistakes considered here are principally 
the ones which sneak up on us unawares and result either in our show- 
ing the wrong relief or in overlooking a structure. Perhaps you prefer 
to call these injudicious estimates, rather than mistakes. 


I. MISTAKES RESULTING FROM TYPE OF SPREAD USED OR HOLE 
SPACING 


You may remember that in the early days of reflection dip shooting 
many companies used single-end or end-on spreads, with holes 2000 
or 3000 feet apart. Each hole was shot in opposite directions along the 
line, and the dips on the opposite sides of the hole were found by sub- 
tracting the weathering correction and “normal’’ move-out from the 
record move-out, with the residual then being converted to a tangent. . 

When these resulting dips were plotted, the cross-section often 
appeared to have a syncline at every shot point as though the under- 
ground had been given a permanent wave. The 1935 cross-section 
shown in Fig. 1 is a fair example of this sort of thing. We soon learned 
that this effect was caused by errors in the assumed velocity used in 
constructing the normal move-out chart. It then became common 
practice to average either the move-outs or the tangents from op- 
posite sides of the same hole so that the cross-section showed the same 
dip for the two sides, thus (theoretically) averaging out the error in 
normal move-out. This, too, was dangerous because in some cases 
there was actually an anticline or syncline present which was erased. 
You may say that this is ancient history—that no one shoots like 
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that today. Nevertheless, many of us have gradually returned to a 
modification of this scheme in present continuous shooting. We place 
our holes 1000 to 2000 feet apart and shoot single-end or end-on pro- 
files continuously down the line, giving so-called continuous subsurface 
coverage. The profiles are computed in pairs with the move-out or 


Sea Level 


Fic. 1. Velocity error in “normal move-out” chart caused these single-end spreads to 
show a syncline at every hole. 


tangent being averaged to obtain the dip between holes. This works 
fine until the going gets tough and one side has reflections which do not 
appear on the reverse profile from the next hole. When this happens 
we are back where we were with the early work—the accuracy of such 
dips (assuming reliable picks) will depend on the accuracy of our 
velocity assumptions for the area. This is an important point, particu- 
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434 JOHN SLOAT 
larly in areas of poor results where reliance is placed primarily on dip 
projection in making the maps. It is important that we realize what 
type of spread was used before we use these unmatched dips in showing 
a nice closure, nose or fault. 

One way to avoid this trouble is to shoot split-spreads across the 
hole, or offset from the hole. But this will require twice as many holes 
if the same length spread is used; if the spread is doubled, this may 
require too much cable or result in too long a spread to get consistent 
picks. 

I recall one case where we were using continuous end-on spreads 
with some of the holes 2640 feet apart. The line computed to be gen- 
erally southwest dip, with dips between one pair of stations appearing 
rather flat. This was noticed but was interpreted to be what it ap- 
peared—merely flat. Unfortunately, Company B shot this line using 
holes goo feet apart—and not even continuous subsurface—and found 
a nice little anticline between our two flat stations. A small oil field 
resulted. Since our dips and depth points were obtained over 2640 foot 
intervals, we could not know that an intermediate point was higher, 
for intermediate traces were not computed because of absence of 
weathering control. Our alibi was that there was a water well between 
the 2640 foot stations that we didn’t want to buy. Company B bought 
a water well but got an oil well in return. 


II. WEATHERING ERRORS 


Here is a problem: Assume a line of holes shot in an area where the 
underground beds are known to be perfectly flat and the surface is 
likewise flat at an elevation of 60 feet above sea level, as sketched in 
Fig. 2. Let us say that the depth of weathering in hole #3 computes to 
be 60 feet deep and holes #2 and #4 compute 30 feet deep from first- 
break refraction intercept plots. Assume the velocity of the low ve- 
locity zone at’ 2000 feet per second and that of the unweathered at 
6000 feet per second, as shown from the refraction plots and up-hole 
information. A correlating reflection on hole #3 shows a record time 
of 3.000 seconds to a shot depth of 60 feet, with the same reflection 
coming in on holes #2 and #4 at a record time of 2.995 from 30 foot 
holes. 

Now compute the depth of the bed in each hole by referring the re- 
flection time to the depth of weathering in the usual manner. The time 
below datum in hole #3 will be 2.970 seconds and holes #2 and #4 will 
be 2.980. Using a run-of-the-mill velocity chart, we get a subsea 
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datum in hole #3 of 14,039—0= 14,039, while holes #2 and #4 show 
the bed to be 14,109 — 30 = 14,079 below sea level. 

Thus we show a nice 40 foot high on hole #3, yet the bed is actually 
flat! This is a familiar problem to many geophysicists, yet it always 
seems impossible at first sight because it seems this difference should 


ACTUAL REFLECTING BED IS FLAT 


WHAT WRONG WITH THIS PICTURE? ANS. RECOMPUTE 2&4 TO 60° on 200 rappers 


Fic. 2. What is wrong with this picture? Answer: Recom- 
pute 2 & 4 to 60 feet and see what happens. 


be taken care of by subtracting the different datum levels from the 
depths. The answer, of course, rests in the fact that we are, in effect, 
applying the difference in datum reflection times on a time-depth 
curve, not at 6000 feet per second where the actual difference occurred, 
but at 14,000 feet per second (in this case) which is the interval ve- 
locity of the curve between the times of 2.970 and 2.980. It turns out 
that if the ratio of the unweathered velocity to the instantaneous 


999 

_-#-- 

| 46000 

1 | 


436 JOHN SLOAT 


velocity of the time-depth curve at the observed reflection time is 3, 
then the difference in subsurface datums will always be equal to but 
in the opposite direction to the difference in computation datum levels. 
In other words, if you compute the same data to two different datum 
levels, you raise or lower the depth of the beds a corresponding amount 
in the opposite direction. 

This problem is extremely important, because it may ee 
gradually over a long distance rather than abruptly; in the former case 
it may result in a miss-tie to wells, and in the latter, it may result in 
showing a local high which is not present. For this reason, some com- 
panies now compute all results to a datum plane, either horizontal or 
inclined. The problem eventually becomes tied in with the question 
of whether the velocity surfaces are flat or inclined, but it is safe to say 
that in a generally flat area serious error can result from variable com- 
putational level. 

One thing that brings this matter to mind is a job some years ago 
where we were able to map, on times, a shallow reflection with con- 
siderable success. The surface was generally flat, and the bed of low 
relief, but the depth of weathering was very erratic. The area had 
been core-holed extensively to the bed being mapped, and just for fun 
we tried to see if we could duplicate the core-hole picture on 10-foot 
contour intervals. The geologists, being of a suspicious nature, refused 
to let us see the core-hole map until after the seismograph map was 
made. The comparison was not favorable, but when the reflection 
times were recomputed to a constant datum level, the two maps were 
virtually identical. Naturally, the geologists claimed this was fudging, 
but it did show the source of our error. 


III. VELOCITY ERRORS 


The error just described brings to mind a boner pulled in New 
Mexico in the early days of reflection shooting. We were mapping a 
bed several thousand feet deep on correlation. Over a certain area the 
reflection showed short times, resulting in a closed high on the con- 
tour map. A well was drilled, finding nothing but a high head of salt 
water in what appeared to be a nice low. It turned out that a thick 
salt lens was present where our presumed high was mapped, and the 
increased velocity through the salt had made the reflection times to 
the deeper reflecting bed come in short over the area of the lens. In 
other words, our so-called high really represented an outline of the 
salt lens, while the deeper beds were actually lower at this point. If 
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we had been a little more critical of the configuration of the contours, 
we might have suspected the true nature of the apparent high, par- 
ticularly since the occurrence of salt was common in this area. 

In this connection one might wonder what sort of contour map 
would result if one were mapping a deep bed, and suddenly ran a 
traverse across some underground cavity such as the Carlsbad 
Caverns. It might be supposed that no reflections would be received 
in this exaggerated case, but there is some evidence to suggest that on 
a minor scale we have run across areas with small underground pockets 
or porous areas giving a subsurface weathering correction which has 
escaped detection, and resulted in mis-mapping of structure. 

Another type of underground velocity condition which may cause 
a mis-mapping of the lower beds has recently come to notice. In certain 
parts of northern California, a basalt flow of variable thickness is 
buried a thousand feet or so below the surface, with true basement as 
much as 8000 feet deep. The objective horizons in the area, the Eocene 
and Cretaceous, occur below the basalt. The question arises as to 
whether the seismograph can map these lower zones, even if good re- 
flections could be obtained. The basalt has a very high velocity and 
since its thickness is variable, it is debatable what value a map on 
lower beds might have. To complicate matters, it has been proved 
conclusively that multiple or echo reflections are obtained from the 
basalt, with this reflection repeating itself in a complex manner 
throughout the record. 

This matter has come to light in another connection. In the 
Southern San Joaquin Valley, where the section is normal but very 
thick, there is a certain place where a nice structure exists in the upper 
beds. The seismograph also indicates a structure in the lower beds (at 
the objective horizon), but these dips look suspiciously like multiple 
reflections from the upper beds. The question now arises as to what 
extent seismograph work is subject to these influences, and suggests 
the possibility that some seimograph structures may have been found 
barren because the later impulses which have been mapped as deep 
reflections may actually be multiples or echoes of one sort or another, 
from unconformable upper beds. There is even the possibility that 
transverse waves have been picked and computed as deeper beds in 
some cases. These problems are a source of joy to the theoretical geo- 
physicists, but are nothing but a headache to the operator who must 
decide whether to sink that $100,000 in a well or pay it in taxes. 

Regional average velocity corrections are now an important part 
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of structure mapping in most parts of California. There are now about 
150 wells surveyed for velocity in various parts of the State, and it 
can be argued with considerable success that closures have been both 
missed and erroneously introduced because of failure to consider the 
effect of lateral velocity gradient on the contours. It is certain that ties 
into wells not over 15 miles apart can be out as much as 2000 feet 
solely because of failure to consider regional velocity change, even 
when the cross-sections are computed from a velocity survey from 
either of the wells. This is an extreme case, but shows a source of error 
in our work. 

Most geophysicists in the Gulf Coast feel that such significant lat- 
eral velocity gradients do not exist in this area, even though the geo- 
_logic section is roughly the same as in California. However, it should 
be borne in mind that the density of well surveys in the Gulf Coast is 
not now sufficient to say for certain what-local effect regional correc- 
tions will have on uncorrected contours. Furthermore, since it is the 
custom in this area to work local prospects as units without tying-in to 
widely separated wells, a great many potential miss-ties have never 
come to light. Perhaps Californians are looking for trouble unneces- 
sarily. 


IV. OTHER MISINTERPRETATIONS 


Numerous oil fields have been missed in the past because insuf- 
ficient study was given the results. Let me give a prime example of how 
this can happen in dip shooting. In Fig. 3 is shown a reproduction of 
a single cross-section. This was a guiding clue to an oil field, but was 
overlooked at the time because of the apparent poor quality of the 
data. Notice that the last two stations on the left hand side show a 
distinct flattening of dip, yet the reflections are mostly graded VVP. 
Additional work proved a reversal about opposite station 9 and re- 
sulted in a prolific oil field. Unfortunately, Company B found this 
field instead of Company A, although Company A shot this line one 
year earlier. It was felt at the time that the observed flattening was 
merely a result of the poor quality of data. A more trusting soul, 
unfamiliar with the limitations of the seismograph, believed this flat- 
tening, and proceeded to check it. 

To give another example, examine the cross-section in Fig. 4. 
Would you have seen anything interesting? Well, this section is the 
clue to an oil field, when looked at with the proper amount of per- 
spective. Notice that station 1 shows a decided flattening of dip. Had 
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this line been extended two more stations to the right, Company A would 
have found an oil field years before Company B got around to tt. 

To give another example of still a different type, examine the cross- 
section shown in Fig. 5. An axis is indicated somewhere near station 5. 
Now the point is that this line is located a couple of miles regionally 
down what was thought to be the plunge from a known field. Few of 
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Fic. 4. Last station on the right shows flattening of dip. If this line had been extended 
two more stations to the right, Company A would have found an oil field. 


us suspected that this was anything other than the final expression of 
the main fold. As a matter of fact, it is a clue to a prolific secondary 
field, separate from the main feature. This single line clearly showed 
this possibility when dips and strikes were computed and plotted for 
these points. The cases of secondary fields being found near known 
fields are numerous in all part= of the country. You may even have in 
your files cases similar to the ones given, or even cases of undiscovered - 
fields. It is a sad fact that once a survey is completed in the field, it is 
apt to become lost in the files and never again consulted until some 
other company brings in a producer. Excuse me while I go consult 
my own files. 
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Another good way to miss a structure, particularly one of low re- 
lief, is to run lines at the wrong angle to the axis. Examine cross 
sections D and Q in Figs. 6 and 7. Cross-section D was shot N-S with 
stations 1000 feet apart. Cross-section Q was shot NE-SW, normal 
to the usual axis of folding in this general region, and shows a distinct 
anticline at the objective horizon of 10,000-15,000 feet. These two 
lines are on the same feature (Line D crosses line Q between station 
Q-14 and Q-15, while this point of crossing on line D is between D-12 


= + 


Fic. 5. Reversal shown on this line turned out to be a secondary oil 
field rather than final down-dip expression of the main fold. 


and D-13). Of course section ““D”’ does show a few flat or North dips, 
but some of these might be interpreted as erratics. This is particularly 
true of station 11 which apparently shows too much North dip 
throughout. The point is that running lines at unsuitable angles to 
the regional trend causes the maximum amount of distortion of the 
true cross-sectional picture. 

In one well known case, Company A discovered a long low relief 
oil field trending Northwest-Southeast. Company B wanted to outline 
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normal to the 


axis of folding. (Compare Fig. 6.) 
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Fic. 7. Running lines at odd angles to regional dip may cause distortion of picture. This section is 
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the structure with the seismograph, and on receiving permission, 
proceeded to run a grid of North-South and East-West lines along . 
roads. The resulting maps showed no closure and very little evidence 
of structure. A new grid was then run Northeast-Southwest and 
Northwest-Southeast, with the result that a structure was clearly 
shown. 

Another way we have probably missed fields is by taking the path 
of least resistance when locating lines on the ground. At this moment 
there may be undiscovered oil fields lying just off of lines which have 
been run and rerun along roads because of topographic convenience. 
When the going gets tough, it is easy to convince ourselves that we will 
get just as good a picture by running along a convenient road rather 
than across that potato patch, wheat field, or marsh land. We forget 
that companies B, C, and D also came to this conclusion, with the re- 
sult that the picture along the road (which may show nothing of inter- 
est) is now clearly defined, but that the picture across the tough area 
may not be known, in spite of all the shooting. 

Today the emphasis in seismograph work is on continuous shooting 
and detail work. A check-up reveals that the majority of companies in 
California, the Gulf Coast and elsewhere are using so-called continuous 
technique exclusively. In this connection, it is important to emphasize 
again that there is actually no such thing as continuous subsurface 
coverage, except in the case of an uninterrupted monoclinal dip, and 
it is to be hoped that we sometimes find something else. Whenever 
flexures or faulting occur, continuous subsurface coverage may cease 
to exist. The records must then be correlated by character, interval 
inspection, and any other system known to science or black magic. 
This point has been stressed before,* but bears repetition. 

As many faults have probably been missed by picking across so- 
called continuous records as were ever missed by wide spaced spot 
correlation records. Most of us tend to force picks through continuous 
records that would never be picked on separate records. The reason 
for this is partly psychological. Picking records, particularly in a tough 
area, is not an exact science, but depends to a large extent on the 
mental attitude of the picker. This was brought into sharp focus a few 
years ago when we shot a line in which the line and cross spread were 


* Joseph L. Adler, Geophysical Exploration for Stratigraphic Oil Traps. GEOPHYSICS, 
Vol. VIII, No. 4, p. 345 (1943). 

M. B. Widess, Multiple Branches in Seismic Reflection-Time Surfaces. GEOPHYSICS, 
Vol. VIII, No. 2, p. 93 (1943). 
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recorded simultaneously on separate halves of the same record. The 
party chief, when picking the records, subconsciously had in mind the 
usual type of record and tended to make his picks line up across the 
record, in spite of the fact that he knew that the two halves repre- 
sented dips at right angles. The resulting cross-section was abnormal, 
and we then re-shot the line using the same line spread but having it 
appear across the record in the usual manner. The new section showed 
entirely different dips from the original. Theoretically, this is impos- 
sible, but in practice it can happen. 

It is even suggested that one reason the seismograph has been slow 
to disclose stratigraphic traps may lie in the fact that most of us like 
to see move-outs line up on our records. We subconsciously pick dips 
which seem reasonable to the eye and which line up with picks at 
earlier and later times. We say to ourselves that that other possibility 
which shows discordant move-out is probably an erratic, and it there- 
fore never appears on a map or cross-section. Reflections from fault 
planes can easily be overlooked in this same manner, because such 
reflections frequently appear on our records with radically discordant 
move-out, often crossing bands of energy from nearer or flatter beds. 
When one considers the fact that the record move-out on a surface 
spread of only 1000 feet must be about 150 milliseconds to represent 
a dip of 70°—75° from a surface only 2500 feet deep, it is easy to see 
how such reflections can be overlooked on the records. This is par- 
ticularly true since reflections from such surfaces are frequently of 
apparent lower quality than other reflections on separate records, and 
since the normal position of the eyes, when picking, distinctly favors 
vertical move-out. 

The psychological factor in picking records is not necessarily 
limited to cases of dip shooting. I well remember being caught in a 
spot correlation boner with “marked records.” After drilling had dis- 
proved the correlation, a post-mortem revealed that the correlation, 
which looked excellent when the reflections were heavily marked with 
colored pencil on the record, looked even better when properly 
matched to fit wells. But it was necessary to consult unmarked records 
to see this better possibility. Once the records have been marked with 
heavy pencil or crayon, it is almost impossible to change the correla- 
tion, unless a fresh set of unmarked records of good quality is avail- 
able. Marking the records obviously does not improve the correlation, 
but does create a psychological impression in favor of the marked cor- 
relation which the eye cannot readily discount. 
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Events of the last few years have convinced many of us that our 
ideas of how to correlate records must be revised and improved. It has 
now been proved, for example, that major faults may exist in lower 
beds without disturbing the upper beds. The question as to whether 
faults may exist in upper beds without affecting the lower zones has 
not been settled, but exists as a possibility. Thus, the old correlation 
criterion of time-interval between reflections will no longer hold in 
many cases. We must become accustomed to anticipating this possi- 
bility. Automatic volume controls have, to a certain extent, robbed us 
of amplitude criteria for character correlations, and the increased 
number of reflections introduced by better amplifiers and filters 

greatly complicates as well as aids correlation. Things have now 
reached the point where the quiet spots on the records may be the 
clues to the proper correlations rather than the reflections themselves. 

One thing is certain. Correlations will never be made unless some- 
one has the time to study the records. The answers are usually on the 
records, but too often they never appear on maps or cross-sections, 
because the party chief—in these days of voluminous bookkeeping 
and exasperating personnel problems—never has time to really study 
the records. The reflections are picked once, a map is hastily made, 
and the data goes in the files to be lost and forgotten until some other 
company finds a field where we had that irregularity. 

This seems a strange procedure indeed when you remember the 
$10,000 a month that it costs to get the information, and the present 
great need for oil. Perhaps an oil field is buried somewhere in your 
files. How about reviewing some of that old “dope” and spending a 

little time in really studying the records? And if time permits after 

writing stuff like this, I will do the same. 


ANALYSIS FOR HYDROCARBONS IN THE PRESENCE 
OF NITROUS OXIDE* 


MONROE W. KRIEGEL{ 


ABSTRACT 


The technique for the analysis of soil air for small quantities of hydrocarbons is re- 
viewed. Various methods of expressing results and the basis for the calculation of each . 
are given. Field data on expansion ratios are presented to show that there exists in the 
soil air a gas, other than the saturated hydrocarbons, which is not removed by the 
screening agents and is retained in the condensation trap. The reasons for believing that 
this gas is nitrous oxide are given, together with plausible sources of this gas in the soil 
air. Laboratory data are presented on the analysis of pure hydrocarbons in the presence 
and in the absence of this gas. Calculations made by each of the methods show that 
when results are based on after-burning readings, a true measure of the hydrocarbon 
content is evolved, regardless of the presence of nitrous oxide. 


I. SCOPE OF THIS PRESENTATION 


In the analysis of soil air for hydrocarbons the supposition had 
long been made in this laboratory that simple reagent traps, such as 
ascarite and phosphorus pentoxide or concentrated potassium hy- 
droxide and concentrated sulfuric acid, would screen out all interfering 
gases and allow an analysis of the material condensed at the tempera- 
ture of liquid nitrogen with the confidence that an accurate measure 
could be made of both the hydrocarbon concentration and the average 
number of carbon atoms in the molecule. In'this paper results of ex- 
tensive field work on the free or unadsorbed gas in the soil air will 
show that at least one interfering gas does exist. The attempts to 
identify this gas will be given, together with laboratory data to show 
how its presence affects the hydrocarbon analysis and methods of 
calculation. 


II. ANALYSIS OF AIR FOR SMALL QUANTITIES OF HYDROCARBONS 


In order to follow the line of reasoning in this paper, it is necessary 
that the analysis procedure used here for hydrocarbons in the air be 
clearly understood. The early stages in the development of an analysis 
unit for this purpose were described by Sokolov.! Pictures and brief 


* Presented at the fourteenth annual meeting, Dallas, Texas, March 1944. 

{ Carter Oil Company, Tulsa, Oklahoma. 

1 Sokolov, V. A., “Gas Surveying,” published in Russian in 1936, translated by 
I. M. Gubkin in 1939. 
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descriptions of this method of analysis have been presented by Hor- 
vitz,? Rosaire,? and Stormont.‘ Patents on the use of this and similar 
equipment have been issued to Horvitz,® who has presented orally® de- 
tails of the analysis procedure. Probably the best published description 
of the arrangement and operation of analysis equipment for the detec- 
tion and quantitative measurement of the small concentrations of 
hydrocarbons is given by Heiland.’? The procedure followed by this 
laboratory is similar to, but not exactly the same as, that described in 
the above references. Hence, to avoid confusion, it might be well to 
describe in detail the analysis of soil air for hydrocarbons as it is done 
here. 

The principle is essentially that of screening out interfering gases 
such as carbon dioxide and water vapor, with proper reagent traps; 
condensing and separating ethane and heavier hydrocarbons from 
methane, oxygen, nitrogen, and other non-condensible gases by the use 
of low temperature and vacuum; burning the condensed gases in an 
excess of oxygen; and measuring the amount of carbon dioxide formed. 
In Fig. 1 is shown a line diagram of the essential parts of a single unit 
of the analysis apparatus. A sample pipette containing soil air with its 
small content of hydrocarbons is attached as shown. With the stop- 
cock to the pipette closed, the entire assembly, with the exception of 
the oxygen supply bulb, is pumped down to the vapor pressure of 
mercury by the mercury diffusion pump. Stopcock 2 is then closed 
and stopcock 1 is opened. The condensation trap is immersed in liquid 
nitrogen, then stopcock 2 is slowly opened and the air and hydro- 
carbons in the sample pipette are allowed to pass through the dry 
reagents, into the condensation trap, and the non-condensibles 
through the diffusion pump. When the pressure in the system has been 
reduced to a few millimeters, mercury valve 3 at the left of the con- 
densation trap is shut and the pressure in the system to the right of 
this valve is reduced to the vapor pressure of mercury. In the con- 


2 Horvitz, Leo, GEOPHYSICS, 4, 210-225 (1939). 
3 Rosaire, E. E., The Handbook of Geochemical Prospecting (1939). 

4 Stormont, D. H., Oil and Gas Journal, 52, September 14, 1939. 

5 Horvitz, Leo, U. S. Patent 2,198,619, issued April 30, 1940, and U. S. Patent 
2,261,764, issued November 4, 1941. 

6 Horvitz, Leo, Presented before Petroleum Division of the American Chemical 
Society in St. Louis, April, 1941, and before the Society of Exploration Geophysicists in 
Houston, April, 1941. 

7 Heiland, C. A., “Geophysical Exploration,” New York, Prentice-Hall, Inc. (1940). 
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densation trap will be retained all of the paraffin hydrocarbons except 
methane. Mercury valve 5 is then closed, the burning trap is immersed 
in liquid nitrogen, and the condensation trap is warmed to —80°C 
by replacing the liquid nitrogen bath with one containing dry ice 
dissolved in acetone. This causes the hydrocarbons in the light frac- 
tion, or ethane, propane, butane, pentane, and hexane, to vaporize 


MANOMETER OXYGEN SUPPLY 
9 


OIFFUSION PUME 


ORY REAGENT TRAP: 


LOW TEMPERATURE BATH MSLEOD GAUGE 


SINGLE UNIT 
OF 
HYDROCARBON ANALYSIS APPARATUS 


Fic. 1. Schematic diagram of a single unit of the hydro- 
carbon analysis apparatus used in this work. 


and to recondense in the burning trap. To get a measure of the volume 
of light hydrocarbons before combustion, the liquid nitrogen bath 
on the burning trap is replaced by a dry-ice bath, and with valves 5 
and 8 open, and 4, 6, and 7 closed, light hydrocarbons are allowed to 
vaporize into the McLeod gauge. After noting the resulting pressure, 
which is called the before-burning pressure (Pgs), the hydrocarbons 
are recondensed in the burning trap by immersing it again in liquid 
nitrogen. 

With valve 4 still shut and with valves 7 and 8 closed, 6 is opened 
and oxygen is slowly admitted through stopcock g until a pressure of 
about 2 mm., as noted on the manometer, has been added to the burn- 
ing trap. Valve 5 is closed and the trap and its contents are warmed to 


STOPCOCK CONDENSATION TRAP 
SAMPLE PIPETTE MERCURY VALVE 
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room temperature. A high frequency discharge is then passed through 
the gas, causing the combustion of the hydrocarbons. After combus- 
tion of the sample the burning trap is again immersed in liquid nitro- 
gen and the excess oxygen pumped out by means of the diffusion 
pump. Then, with valves 4, 6, and 7 closed, and valves 5 and 8 open, 
the liquid nitrogen bath on the burning trap is replaced by a dry-ice 
bath. At this temperature, the pressure of carbon dioxide may be 
measured on the McLeod gauge, while the vapor pressure of the water 
will be negligible. Knowing the volume of the burning trap, manifold, 
and McLeod system, one may, from this pressure, calculate the vol- 
ume of carbon dioxide which has been formed on combustion of the 
sample. This will be explained in the next section. : 
In order to measure the amount of heavy hydrocarbons present in 
the sample, the remaining material in the condensation trap is warmed 
to room temperature. With valve 4 open, valve 5 closed, and the burn- 
ing trap immersed in liquid nitrogen, this portion of the sample ma- 
terial is allowed to diffuse out of the condensation trap into the burn- 
ing trap. The taking of pressure readings and combustion of the sam- 
ple is accomplished in exactly the same way as was done for the light 
fraction. Often measurements of the total hydrocarbon content above 
methane are made. In such cases all of the condensed material is 
transferred in one step, by warming the condensation trap from 
— 196°C to room temperature in one operation. This procedure was 
followed in the analysis of all field samples presented in this paper. 


III. METHODS OF CALCULATION 


The general equation for the combustion of hydrocarbons is as 


follows: 


: = + (n + 1)H.0. (2) 


CrHente + 


As explained in the preceding section, measurements in a system of 
known volume of the pressure before-burning (Pgs) before the addi- 
tion of oxygen gave a means of determining the absolute amount of 
the pure hydrocarbon present, while the after-burning pressure 
(Paz) when measured after the non-condensible gases had been 
pumped off and with the burning trap at —80°C, such that water 
vapor would have a negligible vapor pressure, gave a measure of the 
amount of CO, formed. Throughout this paper the terms “‘before- 
burning pressure,” (Pg), and “after-burning pressure,” (Pas), will 
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refer to measurements made as described above and in the preceding 
section. The percentage of hydrocarbons in the gas would be expressed 
by: 

Pas‘ V, 


X 100 = % hydrocarbons (2) 
-760:V> 


where V,=volume of the system in which the pressure measurement 
is made, V»= volume of the sample in the pipette, and Pgz is the 
before-burning pressure in microns. Since the percentages are often 
very small, it is usually more desirable to express the results j in terms of 
parts per million by volume, as follows: 


Ppp: Vs 


—_——— 10° = P.P.M. hydrocarbons. (3) 
10°: 760: Vp 


Substituting & for V,-108/10*- 760-V, for a system of constant size 
and for a constant sample volume, Eq. 3 becomes 


= P.P.M. hydrocarbons. 3’) 


Results may also be calculated on the basis of the after-burning read- 
ings, as follows: 


Pap 


n 


-k = P.P.M. hydrocarbons (4) 


where »=number of carbon atoms in the hydrocarbon molecule, and 
P4g=after-burning pressure in microns. Since the light hydrocarbons 
as separated consist of ethane through hexane, it is convenient in 
practice to express the results on the basis of an average hydrocarbon, 
such as butane, in which equation (4) becomes: 


Pap 
4 


-k = P.P.M. light hydrocarbons. (4’) 


The heavy hydrocarbons may be expressed as octane, or: 


Pas 
: -k = P.P.M. heavy hydrocarbons. 


(4””) 


There is always the possibility that a small amount of COs: will pass 
the ascarite and be retained in the condensation trap. Since the volume 
of CO, does not change when exposed to the discharge in an excess of 
oxygen, the effect of CO, in the condensation trap would be to in- 
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crease both Pgz and P,z by a constant amount. Hence, to compensate 
for the possible presence of this material, calculations may be made as 


follows: 
(Pas — Ppp) 


k = P.P.M. hydrocarbon (5) 
or, if results are still expressed as butane— 
Pap — P 
(Paz Bp) -k = P.P.M. light hydrocarbon. (5’) 
3 


A consideration of the above equations shows that for pure hydro- 
carbons, the parts per million of hydrocarbons may be calculated by 
equations 3, 4, or 5. Equation (3) is based on the premise that the ma- 
terial in the condensation trap is hydrocarbon only. Equation (4) has 
as its basis that the products of combustion will give a true measure of 
the hydrocarbon concentration, while equation (5) is based on the sup- 
position that any material in the condensation trap other than hydro- 
carbon will either be non-reactive or will form materials having the 
same volume as the initial materials under the conditions of after- 
burning measurement. In the absence at this point of evidence of 
other interfering gases and in view of the possibility of getting CO 
in the condensation trap, equation (5) would seem to be the most de- 
sirable for general use in calculating the results of measurements made 
on soil air. Data will be presented later in this paper to make such a 
conclusion untenable. 

The term, “expansion ratio,” in this paper will be used to express 
the ratio of after-burning to before-burning readings, or Pas/Psz, in 
which both readings are made as pointed out in the procedure. On this 
basis, the expansion ratio of pure ethane would be 2, propane would 
be 3, or would be in the general equation (1) CnHon42=mCOs2. Since 
methane is not retained in the condensation trap, it is apparent that 
the minimum expansion ratio possible on combustion of pure hydro- 
carbons would be 2. Even if the ascarite traps were to fail to remove 
the COs, the expansion ratio of a mixture of CO, and hydrocarbons 
would approach 1 as a limit. 


IV. FIELD DATA ON EXPANSION RATIOS 


In connection with other geochemical studies, samples of soil air 
were taken from a number of areas in Oklahoma and other states. 
Sample locations were not less than .1 or more than .5 miles apart. 
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The procedure for the collection of the sample material was to dig a 
six-foot hole, insert a small pipe to within a few inches of the bottom, 
pack off the upper part of the hole, and connect evacuated pipettes to 
the top of the pipe. After withdrawing sufficient soil air to cleanse the 
sampling line of surface air, samples were retained for analysis. 

In Table I are shown the expansion ratios obtained on analysis of 
soil air samples from areas in Oklahoma, Kansas, Arkansas, and 
Louisiana. In Fig. 2 are shown the distribution curves for the data in 
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EXPANSION RATIOS —EXPANSION RATIOS EXPANSION RATIOS EXPANSION RATIOS 
AREA AREA 2 AREA 3 AREA 4 AREA S AREA 6 AREA 7 


DISTRIBUTION CHART OF THE EXPANSION RATIOS IN THE SEVEN AREAS 


Fic. 2. Schematic presentation of the expansion ratios. In each area the number of 
values lying between o-1, 1-2, etc. is plotted as ordinates. 


each of the seven areas. To illustrate how these were prepared, con- 
sider area 2. Here ten of the values have expansion ratios of greater 
than o but less than 1.00, eleven are greater than 1.01 but less than 
2.00, eight between 2.01 and 3.00, etc. Thus the various vertical lines 
represent the number of expansion ratios lying between arbitrarily 
selected limits. In order to show the variation between locations, the 
data on areas 5 and 6 are shown graphically in Figs. 3 and 4 respec- 
tively. In both of these two areas samples were taken along profiles — 
on .2 mile spacing with continuous lines in Fig. 3 along locations 1 
through 6, 7 through 27, and 28 through 50, and in Fig. 4 along 1 
through 19 and 20 through 50. 
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TABLE I. EXPANSION RATIO VALUES 
Area Number 
Number 3 5 6 7 
I 5.6 1.52 48 1.32 96 3.0 -245 
2. 1.3 .86 5.0 -65 2.0 4.1 75 
3 5.0 1.63 1.54 65 
4 8.8 2:3 255 II 1.70 -95 
5 4-33 1.72 17 1.46 10 3.95 2:0 
6 6.3 3.14 2.92 1.47 +079 66 1.06 
7 2.91 1.46 2.44 1.05 205 
8 1.45 222 252 1.40 1.9 2.14 
9 8.0 1.31 1.86 Ansa 1.68 24 2.65 
10 2.15 2.61 2.9 17.8 
II -6 1.2 2.66 .20 1.24 4-77 
12 | 3027 -30 2.28 4.3 
13 1.45 1.33 +30 1.64 1.41 
14 10.3 267 1.38 3.48 54 
15 TA 4.23 1.63 5.28 -44 4.5 
16 2599 1.07 4-55 1.6 
17 3-0 3-95 4 2.55 2.59 6.95 1.9 
18 4.0 .18 2.48 10.2 
19 2.56 4-57 1.47 1.93 235 
20 1.07 1.92 2.79 1.96 
21 1.44 -93 14 .78 1.78 ax! 1.48 
22 1.07 2.25 -43 234 3-4 2.68 
23 2.0 2.43 2.14 1.16 1.98 1.28 E72 
24 5.8 96 2.10 BG 2.2 
25 Lid 8.35 23 3.78 .69 6.9 1.78 
26 3.0 3.6 1.0 1.6 1.06 2:8 Par 
27 3-45 4-7 1.41 3-03 2.47 3-4 5.6 
28 2:5 3.47 4.3 2.52 
29 5.6 1.18 1347 4.0 3.20 4.65 .68 
30 6.45 3527 2.88 233 2.6 1.46 
31 3.0 24 4.8 3:0 4.0 
32 Bot 2.55 DotA II .40 2.4 2.3 
33 6.7 1.35 1.56 1.9 
34 10.1 .29 .034 2°. 86 
35 9.8 3-43 3.28 .10 1.87 ase 1.6 
36 3-25 4.21 2.9 2.67 Bag 5-33 1.6 
37 1.8 1.43 I.41 2.04 
38 1.8 3.30 9.3 2.4 -64 
39 223 7.0 3.96 1.45 18 4.85 2.69 
40 1.44 4.27 1225 21.5 .009 3-93 
41 163 4.4 17 5.76 .28 
42 2.0 1.24 .728 -29 36 6.73 3.08 
43 ow 4 1.54 -66 27 2.6 2.28 
44 5.43 4.1 -98 1.04 712 4.2 1.30 
45 1.66 -16 2.29 47 7.55 
46 .044 29 5.0 4.55 
47 2.7 1.1 1.48 6.1 232 
48 3.0 ASS 1.16 1.67 138 2.26 <a5 
49 1.5 3.86 3.48 14 48 4-45 1.33 
50 9.4 1.02 1.00 745 233 1.05 
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TABLE II. SuMMARY OF EXPANSION RATIO DATA 


pe Average Expansion Ratios Per Cent of 
Number All Values Below 2 Above 2 Total Below 2 
I 3-95 1.362 5-15 a 
2 2.67 975 3-89 42 
3 1.903 1.148 3.25 64 
4 © .236 1.016 3-92 54 
5 955 -677 2.99 88 
6 3.17 -992 3-97 28 
1.18 4.30 64 
Average 2.474 1.500 3-924 53-1 


It will be noted in both the figures and in the individual expansion 
ratios that there is no uniformity in the ratios in any single area. It 
will also be noted that the expansion ratios of many of the samples are 
less than 2, and that a considerable number are less than 1. 
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Fic. 3. Showing the expansion ratios obtained in a typical series of 
samples. The sample locations are .2 mile apart. 
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In Table II are shown the averages of the individual values in 
Table I. Here are presented the average expansion ratios of all the 
samples in each area, the average of the ratios above 2, and the average 
of those below 2. Included in the table is the percentage the number of 
sthtions showing values below 2 is of the total 50 stations in each area. 
The averages show that in the seven areas studied, 53.1 per cent, or 
from 28 to 88 per cent, of the samples of soil air contained a gas other 
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EXPANSION RATIOS OF SOIL AIR SAMPLES | 
IN AREA 6 


> 


EXPANSION RATIO-Rig/ Pag 


¢ s 10 5 20 25 30 35 a0 45 $o 
SAMPLE LOCATION NUMBER 


Fic. 4. Showing the expansion ratios obtained in a series of samples taken in another 
area. Most of the values are along lines with locations .2 mile apart. 


than pure hydrocarbons. This is a minimum value, since it does not 
take into account those samples in which the expansion ratio may have 
been reduced somewhat but not enough to bring its value below 2. 
Some of the values are extremely low, for example number 18 in area 
3, number 46 in area 4, number 16 in area 6, and many others. 

It is apparent from these data that the gas causing the reduction 
in expansion ratio is present in non-uniform concentrations in many 
areas. It will be remembered that the hydrocarbon content of the 
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samples, which is not discussed in this article, will also influence the 
expansion ratio. However, in general it was observed that the varia- 
tion in hydrocarbon content among samples in any one area was not 
sufficient to account for the large variation in expansion ratios. Since 
no areal correction could be applied, it was important that the gas 
be identified in order that provisions might be made for its removal 
ahead of the condensation trap or for a method of calculation in which 
its presence would not influence the results. 


V. IDENTIFICATION OF THE CONDENSED GAS 


The gas causing the decrease in expansion ratio must be one which 
has certain definite properties. It must not react with ascarite or phos- 
phorus pentoxide, must have a vapor pressure at low temperature 
similar to the light hydrocarbons (examination of many high field 
samples showed it to lie in the light fraction) which have a boiling 
point range of —110°C to +35°C, and must on combustion in excess 
oxygen form a material which is either not condensed at — 196°C or 
is not volatile at — 80°C. The latter point is important to note care- 
fully. If the product of combustion is non-condensible at the tempera- 
ture of liquid nitrogen, it will be pumped out with the excess oxygen. 
On the other hand, if the reaction product has a boiling point in the 
same range as water, it will not be volatile at —80°C, and hence will 
not be measured with the CO; in the regular after-burning reading. 

According to Scott’ a mixture of gases unadsorbed by most of the 
common reagents might contain C2He, C2Hs, N20, Ne, 
COS, and the noble gases. Hence it might be well to consider these 
initial possibilities. For convenience in following these discussions the 
boiling points’ at 760 mm. of all of the gases entering the considera- 
tions are given in Table III. Within the paraffin hydrocarbon series 
the normal boiling points give a relative measure of the vapor pressure 
characteristics at —80°C and —196°C. This analogy is not strictly 
true when applied to other gases, but is used as the basis of the division 
into the two groups in Table III. The ones which have too high boiling 
points at 760 mm., or too high vapor pressure to be retained in the 
trap at a temperature of — 196°C and a total pressure equal to the va- 


8 Scott, W. W., “Standard Methods of Chemical Analysis,” New York, D. Van 


Nostrand Co. (1939). 
® Hodgman, C. C., and Holmes, H. N., “(Handbook of Chemistry and Physics,” 
Twenty-sixth edition, Cleveland, Ohio, The Chemical Rubber Publishing Company 


(1942). 
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por pressure of mercury, are H2, Ne, and CHy, and all of the noble gases 
except Xe and Rn. Since C2H2 and C2H, would readily burn to CO2 
and cause an expansion in volume, they may be eliminated. Since 
xenon and radon are non-reactive, they would not satisfy the require- 
ments in that the volume would remain constant. Carbon oxysulfide 


TABLE III. Borminc Points oF GASES AT ONE ATMOSPHERE 


Light Hydrocarbon, Volatile at —80°C, 


Methane, or Non-Condensible Fraction Dien 
Name Formula B.P.° C Name Formula B.P.°C at 

760mm. 760 mm. 
Methane CH, —161.6 Ethane C:H¢ —88.5 
Oxygen Oz —183 Propane C3Hs —42.5 
Nitrogen N2 —196 Butane C.Hio —0.5 
Carbon Monoxide CO —190 Pentane CsHiz +36.1 
Hydrogen . —252.8 Hexane +68.8 
Nitric Oxide NO —151.8 Carbon Dioxide CO2 —78.5 
Argon A —185.7 Nitrous Oxide N.O —89.5 
Neon Ne — 245.9 Carbon Oxy-Sulfide COS —50.2 
Helium He —268.9 ° Xenon Xe —107.1 
Krypton Kr —152.9 Radon Rn —61.8 


has the proper boiling point, but according to Dennis’ this material 
will decompose with heat in the presence of oxygen according to the 
following equation: 


COS + $02 = CO2z + SOx. 


Since both of these products of combustion are condensed at — 196°C 
and volatilized at — 80°C, they will appear in the after-burning reading 
and cause a net expansion ratio of 2. Thus, all except nitrous oxide 
may be eliminated as possibilities. 

The boiling point of nitrous oxide is such as would place it in the 
light hydrocarbon fraction. Joshi"! has made a comprehensive series 
of investigations into the dissociation of nitrous oxide in the silent 
discharge at pressures above 10 cm. mercury, finding the dissociation 
to be according to the following reactions: 


2N20 = 2N2 + On 
2N20 = 2NO+ Ne 
2NO = Ne + Oz. 
10 Dennis, L. M., and Nichols, M. L., “Gas Anlaysis,”” New York, The Macmillan 


Company (1929). 
11 Joshi, S. S., Trans. Fara. Soc. 23, 227-38 (1927); 25, 137-43 (1929). 
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Stewardson” found the total pressure after combustion to be 1.5 times 
the initial pressure, indicating that decomposition is complete, i.e. 
equilibrium is reached only when all of the nitrous oxide is decom- 
posed, there being no appreciable reverse reaction. 

Although these reactions were not carried out in an excess of oxy- 
gen, it seemed reasonable to assume—subject to laboratory investiga- 
tion—that the reaction products would be the same. An examination 
of Table II will show that none of the reaction products are con- 
densible at the temperature of liquid nitrogen. Thus, nitrous oxide 
seems to satisfy the requirements for the type gas found in the soil air. 

The presence of nitrous oxide in the soil is not surprising when it is 
remembered that the element nitrogen in the form of ammonium salts, 
nitrites, or nitrates is present in fertile soils. One method" of prepara- 
tion of pure nitrous oxide is according to the equation: 


NH,NO; -+ 2H20. 


It is probable that the slow decomposition of commercial fertilizers 
might account for some of the nitrous oxide in farming areas. In con- 
nection with the studies of decomposition of vegetation under aerobic 
conditions it has been shown by this laboratory that a gas having 
properties similar to nitrous oxide forms a large portion of the con- 
densed fraction. 


VI. LABORATORY STUDIES ON NITROUS OXIDE AND HYDROCARBONS 


Having found that nitrous oxide has such properties as would 
favor a behavior in the presence of hydrocarbons similar to that of the 
unknown gas, and having established that the presence of this gas in 
the soil air is not unusual, it remained to test the effect of this gas on 
the analysis of samples of pure hydrocarbons. In Table IV are pre- 
sented the results of analysis of pure propane and pentane alone and 
in the presence of excesses of nitrous oxide. In separate five-gallon 
jugs were prepared samples of propane, pentane, and nitrous oxide in 
the concentration ranges often encountered in geochemical soil air 
samples. The initial concentration of hydrocarbon in the jug was de- 
termined in each case by withdrawing and analyzing a 335 cc. sample. 
Calculations were made on the basis of the before-burning pressure, 


12 Stewardson, E. A., Trans. Fara. Soc. 30, 1018-1027 (1934). 
13 Partington, J. R., “A Text-Book of Inorganic Chemistry,” London, Macmillan 


and Company (1939). 
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the after-burning pressure, and the difference between the before- and 
after-burning pressures, or according to Eqs. 3, 5, and 4, respectively. 
Then, an analysis was made of a sample from the same jug mixed with 
a sample from one of the jugs containing nitrous oxide. This was done 
by passing in one 335 cc. sample of air containing the hydrocarbon, 
condensing and transferring to the burning trap, making a before- 
burning reading, passing in a 335 cc. sample of air and nitrous oxide, 


TABLE IV. RESULTS OF ANALYSIS FOR PROPANE AND PENTANE IN THE 
PRESENCE OF NITROUS OXIDE 


kPAB 


k(Pap—Ppp) 
Material k Pps Pap Expansion 
Ratio n 


V,—108 PAB 


———  Microns Microns —— P.P.M. PPM. P.P.M. 
760° Vp* 108 PBB 
Propane No. 1 Alone +2345 20.1 56.7 2.82 47.1 42.9 44.3 
Propane Fraction 22345 19.2 45.0 
Nitrous Oxide Fraction 22.5 
Propane+Nitrous Oxide 41.7 1.33 97.8 43-1 
Propane No. 1 Alone +2370 20.2 57-4 2.84 47.8 44.1 45.3 
Propane Fraction 2370 19.2 45-5 
Nitrous Oxide Fraction 26.2 
Propane+Nitrous Oxide 45-4 56.8 1.25 107.5 13.5 44.8 


26.3 


Propane No. 2 Alone 2 .89 

Propane Fraction +2342 10.8 25.3 

Nitrous Oxide Fraction 12.5 

Propane+Nitrous Oxide 23-3 23.7 1.45 59-3 12.2 26.3 


Pentane No. 1 Alone .65 19.05 22.1 21.5 

Pentane Fraction +2345 8.33 19.5 

Nitrous Oxide Fraction 22.00 

Pentane+Nitrous Oxide -54 


Pentane No. 1 Alone 
Pentane Fraction +2370 9.04 
Nitrous Oxide Fraction 23.60 
Pentane+Nitrous Oxide 


Nitrous Oxide+Air +2342 12.9 4.20 -325 30.2 
Nitrous Oxide+-Air +2345 24.5 3-68 .150 57.5 —12.2* 
Nitrous Oxide+-Air 1.82* 


SE 


* Calculated as Pentane 


condensing and transferring, making a before-burning reading of the 
combined materials, and then adding oxygen and burning and measur- 
ing the products of combustion in the usual manner. To illustrate, 
consider the first unit in the table, which is propane number one and 
nitrous oxide. The initial concentration of propane in the jug, as de- 
termined on a sample withdrawn for such analysis, was found to be 
47.1, 42.9, Or 44.3, as calculated from equations 3, 5, and 4, respec- 
tively. Those results are all within the limit of experimental accuracy. 
A second propane sample from the same container showed a concentra- 


21.4 
49-4 im 77-4 9-93 23-4 
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tion of 45 P.P.M., as calculated on the basis of the before-burning 
reading, or using Eq. 3. The nitrous oxide from a 335 cc. sample taken 
out of one of the containers of that gas brought the initial pressure 
up to a value of 41.7 microns. Considering only the total before- and 
after-burning pressures, calculations were made as if this were a regu- 
lar field sample in which all of the condensed gas was hydrocarbon. 
On this basis the values of the hydrocarbon concentration were found 
to be 97.8, 11.5, and 43.1 P.P.M. by the three methods of calculation. 

In the last three experiments shown in Table IV nitrous oxide was 
added to air containing a very small amount of combustible material. 
The results indicate abnormally high values for the apparent hydro- 
carbon content when calculated on the before-burning basis, negative 
values when calculated on the difference basis, and normal values 
when calculations are based on the after-burning reading only. 


VII. CALCULATIONS OF HYDROCARBON CONTENT IN 
THE PRESENCE OF NITROUS OXIDE 


It is now apparent that from both a theoretical and experimental 
basis, no method of calculation of hydrocarbon content can be used if 
it is based on either the before-burning or on the difference between 
before- and after-burning basis when nitrous oxide is present in the 
sample material. Thus, the equation for calculation in either the pres- 
ence or absence of nitrous oxide is: 


Pap: 10% 


n- 10°: 760-V, 


= P.P.M. by volume of hydrocarbons in air 


where 


Paz=After-burning pressure in microns 
n= Number of carbon atoms in the molecule 

V,= Volume of the system in cc. 

V,= Volume of the sample in the pipette in cc. 


VIII. SUMMARY 


1. The details of the procedure for the analysis of air for small 
quantities of hydrocarbons have been reviewed. The various methods 
of calculation and means of expressing results, under the assumption 
that there were no interfering gases, have been given. 

2. Field data on expansion ratios have shown that there is often 
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present in soil air a gas which interferes with certain of the measure- 
ments and methods of calculation. 

3. A study of the literature and of properties of the disturbing gas 
indicates that it is nitrous oxide. Plausible reasons for expecting to find 
this gas in the soil air are given. 

4. Laboratory data on combustion of pure propane and pentane 
with nitrous oxide indicate that the method of analysis and one method 
of calculation is satisfactory for hydrocarbons in the presence or ab- 
sence of nitrous oxide. 
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QUANTITATIVE INTERPRETATION OF MAGNETIC AND 
GRAVITATIONAL ANOMALIES* 


ERVAND G. KOGBETLIANTZ{ 


ABSTRACT 


A new interpretation method for gravitational and magnetic anomalies based on 
the use of average values (integrals) of suitably chosen functions of the observed quan- 
tity is discussed and illustrated by application to the particular case of a symmetric 
anticline. The quantitative interpretation yields the position of the apex, the slope of 
the sides, the thickness, the depth of the base of the anticline, as well as the density con- 
trast for the gravitational case and the magnitude and direction of the magnetization 
vector for the magnetic case. 


INTRODUCTION 


Geophysical surveying became necessary only when geological 
prospecting had mapped almost all outcrops and discovered all 
mineral deposits and oil fields sufficiently shallow to be found by 
purely geological methods. Only the coming exhaustion of such de- 
posits discoverable by the study of the tectonic structure of the upper- 
most part of the earth’s crust compelled the prospectors to change 
their methods, to study and apply geophysical methods in prospecting. 

During its first twenty years, 1920-1940, geophysical surveying 
was directed mostly towards the study of relatively simple and easy 
cases. Excessively deep or excessively small deposits and gentle tec- 
tonic structures which create in general small anomalies were neg- 
lected, and exploration geophysicists looked for and studied clearly 
defined anomalies which correspond to shallow deposits and to big, 
well-defined tectonic structures. On the other hand, the regions char- 
acterized by complex anomalies resulting from superposition of many 
simple anomalies were declared as rather unsuitable for geophysical 
prospecting for the too simple reason that the interpretation of the 
maps of complex anomalies was considered almost impossible. 

This situation was not serious as long as the simple anomalies of 
fairly large magnitude were not yet exhausted. But now geophysicists 
are forced to deal with more difficult cases, and they have to interpret 
maps of small as well as of complex anomalies. The usual interpreta- 
tion methods, which are based on the systematic use of individual 
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critical values (extrema, zeros, inflection-points, etc.), and of their 
distances, helpful as they are in easy cases, are useless and cannot be 
successfully applied in the general case. In fact, they have been 
elaborated on the essential assumption that the map under considera- 
tion is that of a simple anomaly created by a single geologic phenome- 
non, i.e., an isolated mineral body or tectonic structure. But a general 
method for the resolution of a complex anomaly map into partial maps 
of component simple anomalies, whose superposition forms the ob- 
served and mapped anomaly, is lacking. Furthermore, the usual inter- 
pretation of a simple anomaly with the aid of critical values and dis- 
tances is based on the tacit assumption that these values and distances 
reflect exclusively the action of the unknown structure or ore-body 
whose study is the ultimate purpose of the interpretation. This as- 
sumption is true approximately for an anomaly of large magnitude but 
it becomes questionable if the anomaly is of average magnitude and is 
completely wrong for small anomalies. 

New interpretation methods are necessary and they must be intro- 
duced into practice if geophysical prospecting by gravimetric and 
magnetic methods is to be usefully applied in the future. These new 
methods are based, like the old ones, on the theory of potential, and 
they solve the two fundamental problems: (A) resolution of a complex 
anomaly map into maps of simple component anomalies, and (B) 
quantitative interpretation of small anomalies. 

In this paper we study the problem (B) in the case of an axial 
anomaly created by an anticline whose axis is horizontal. All results 
given here are based on the assumption that the cause of the anomaly 
is an anticline. If this fact in a concrete case cannot be proved by 
geology or by existing borings, it can be deduced from the study of 
the anomaly map itself by one of the existing general methods es- 
pecially devised for this purpose by the author and which solve the 
problem (A). Consequently, the interpretation rules formulated in this 
paper are to be considered as a method of detailed quantitative inter- 
pretation and they must be applied only after the first interpretation 
has shown that the anomaly corresponds to an anticline. 

In the two-dimensional case we are considering here there is a 
classical Poisson’s relation between the two components Z (vertical), 
X (horizontal) of the magnetic anomaly, and the gravity gradient 
U,.=G and curvature U,,=K, (Us= Uy = Uy,=0): 


+ iZ) = + 


- 
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where f, o, J, y denote the gravitational constant, density contrast, 
intensity of magnetization, and angle of inclination of the magnetiza- 
tion vector, respectively. By means of this relation we transform the 
interpretation rules formulated in this paper for the case of G- and 
K-maps into interpretation rules for Z- and X-maps of a magnetic 
anomaly created by an anticline. All considerations developed below 
with respect to gravity anomalies hold equally for magnetic anomalies. 


1. Punctual Anomaly 


Let us fyrst study the obstacles which naturally lead to an er- 
roneous interpretation of a given map if the magnitude of the anomaly 
is not very large. We do not have in mind the well known theoretical 
fact of nonuniqueness of solution.of the interpretation problem. We 
are not concerned here with this purely theoretical objection. The 
map itself, picturing the behavior of any one among gravitational 
quantities: Dg (variation of the gravity acceleration g), Uzz, Uyz, Ua, 
2U,., (second derivatives of the potential U), gradient G, curvature K, 
total anomaly A =(G?+K”)!, is supposed to be generally correct and 
to have had applied to it all the usual reductions. A map for Dg, for 
instance, must be corrected for latitude, free-air, Bouguer and terrain 
corrections. Such a map is plotted on the basis of measurements made 
at isolated station-points. All usual corrections having. been made, 
each individual value of the observed quantity obtained at a station 
is the combined effect of three gravitational actions: 1) regional anom- 
aly, 2) anomaly caused by the tectonic structure or the ore-body the 
study of which is the aim of the interpretation and 3) anomaly gener- 
ated by local irregularities of mass distribution in the immediate 
vicinity of and under the point of measurement. This local anomaly 
—we propose to call it a “(punctual anomaly”—is in general very small 
and affects only a small area around the point-station. It is this punc- 
tual anomaly which is responsible for perceptible variations in the 
station-value of the observed quantity which occur for small displace- 
ments of the apparatus used through two or three yards. In fact, the 
apparati used now in gravimetric and magnetic surveys are extremely 
sensitive and we cannot safely neglect the existence of punctual anom- 
alies. 

A map can and must be corrected for the regional anomaly, but 
there is no correction at all for the punctual anomaly affecting each 
observed value: it cannot be evaluated. After the correction for regional 
anomaly is made two circumstances are possible. If the magnitude of 
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the anomaly studied is large the punctual anomalies are negligible 
and the map can be interpreted by the usual methods. The positive 
results achieved by the old interpretation methods which use the in- 
dividual values are explained in this way. But if the magnitude of the 
anomaly is small the punctual anomalies not only modify the extremal 
values and all the ratios, but they also displace these critical values, 
altering all the distances used in the usual interpretation. This impor- 
tant fact explains the lack of success in dealing with maps of small 
anomalies. 

If we wish ta eliminate the errors of interpretation caused by 
punctual anomalies, there is only one possible way to do it. Considered 
altogether, the punctual anomalies ‘in the region mapped have a ran- 
dom distribution: they oscillate about zero and are independent one 
from another. Consequently, they must undergo an almost total com- 
pensation if we form the average value of some function of observed 
and mapped quantity, provided that the regional anomaly has already 
been removed. Using all the map at a time, combining all observed 
values in an integral, we eliminate the harmful influence of punctual 
anomalies. For any map in general the average values express much 
better the gravitational action of the phenomenon under study than 
do the individual values and their distances. Formulas based on the 
individual values hold only in very rare and exceptional cases of big 
structures such as the salt domes of Texas, Kursk, and the Kiruna- 
waara iron ores or the Great Rhodesian Dyke. 

In this paper we introduce a new interpretation method in which 
the interpretation rules use exclusively the average values and, in par- 
ticular, moment-functions and moments of the observed quantity 
and of its square. The method has been applied in practice by the 
author in France and in Iran. The case of an anticline is chosen here 
only for the sake of brevity; the author has developed his method for 
the most general cases of complex anomaly maps obtained as a result 
of a magnetic or gravimetric survey. 


2. Formulation of the Interpretation Problem 


Taking the vertical directed downward as the Z-axis, let us con- 
sider a cylindrical body of normal cross-section S, having a density 
contrast o and stretching parallel to the Y-axis, the horizontal plane 
XOY being that of the map (Fig. 1). The section S is in the plane XOZ 
and the potential U of the body does not depend on the coordinate y: 


Bi: 
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U=U(x, 2). Computing the difference U(x, 2)—U(xo, 2), where 
(x0, Yo, 20) is any fixed point, we get 
ds. f [Rot — R]dn (2) 
A 
with 
f = 66.7 X 107° cgs 


0 Line of Prorice 


Fic. 1. Illustrating the quantities used in the formulation of the interpretation 
problem. An infinite cylinder lies with its generatrices parallel to the Y axis. C(Z, 2), 
the centroid of the cross-section S of the cylinder, lies at a horizontal distance ¢ from 
the origin O; and at a depth 2 below the surface. I’ is the bounding curve of the cross- 
section. 


being the gravitational constant, and &, », ¢ being the coordinates of 
a point of the cylinder. 
Here 


f — R™|dn = 2 log (r/ro) 
if 7 
= (x — and 10? = (x9 — &)? + (20 — 
Denoting the constant 2fo by k, we get from (1): 


U(x, z) = ef f log r-dS + constant, 
A 


x 
a 
CROSS-SECTION S Vi 
0:0=X 
Curve Q:0=X 
oc=Z 
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where the double integration is extended over the area A of the cross- 


section S. 
The partial derivatives Uz, Dg=U, of the potential U considered 
at the surface, z=o, have at the point (x, 0, 0) the values 


and combining them in a complex expression U,—iU, we get the im- 
portant formula 


iDg = eff (x — = kf tog (x — (3) 


where p=&+7¢ is the radius-vector of the point (¢, ¢) which describes 

the area A of the cross-section S or the curve I’, the boundary of S. 
Differentiating the relation (3) with respect to x, we have also an 

expression involving the gravity gradient G and the.curvature K, 


K-ic=af f p)-tds = & f (x — p) dt. (4) 
A r 
Using the binomial expansion, we find for large | «| : | 


(8/2) of J = 
(RA/x) [x + (4 + | 
= + 2) + |], 
where # and 2 are coordinates of the center of gravity C of the area A. 


We see that, for large | «| , U,, Dg, Gand K for any form of the section 
S are of the orders | «| i aa | «| — and x respectively: 


DgwkAz/x?; G~—2kAz/x®; (5) 


Before applying the general formulas (3) and (4) to an anticline 
we have to choose some simple geometrical form of its cross-section S 
well adapted to express adequately its gravitational action and at the 
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same time amenable to mathematical computation. For this purpose 
we consider a symmetric anticline as a single infinite disturbing layer 
of constant density contrast o and of original thickness T (see Fig. 
2). For | «| =I this layer remains horizontal, but in the interval 
—l<x<l=EC its (originally also horizontal) part, of length 2lo, 
(lo)>1), compressed by tangential forces, forms an anticline of cross- 


LINE OF PROFILE 


x 
x 


Fic. 2. Illustrating the quantities required for application of the method of inter- 
pretation to an idealized representation of a symmetric anticline. 


section S limited by the segments of straight lines D’C’, C’A, AC, 
CD, DB, BD’. The excess-mass 2(/)>—/)-T-o and the maximum 
thickness AT=AB(A>1) are related and (1+A)-/=2l). We char- 
acterize the cross-section by the parameters a, B, y, 6 where a, B, 
go°— , 90°— 6 are the angles of the lines AC, BD, OC, OD with the 
horizontal axis OX (the line of profile). Denoting the depth OA of the 
apex A by pand the distance AE (displacement of the apex upward) 
by qg, we have 


l=q cot a=[q—(A—1)T | cot B=(p+ 9) tany= T) tané. (6) 


Therefore, /, g; \ and T can be expressed as functions of a, B,.y, 6 and 
in all there are seven parameters &, (i.e., c), a, B, y, 5, &, p. The last 
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two are coordinates of the apex of the anticline with respect to an 
arbitrary origin O; on the profile. The gravitational anomaly of an 
anticline is due not only to the excess-mass of its section C’/ACDBD’C’ 
of area Aj, but also to the suppression of the mass 2:/-0-T of the hori- 
zontal section C’/CDD’ in —/Sx<l. Therefore in what follows we 
ascribe to this cross-section C’/CDD’ the negative density-contrast —o, 
denoting its area by A». The final purpose of the quantitative interpre- 
tation in the case under consideration is the numerical determination 
of seven parameters which means location of the anticline and knowl- 
edge of its density contrast o, slope a, thickness T and AT, height q, 
width 2/ and depth # of its apex. 

We turn now to the interpretation problem of the magnetic survey. 
Poisson’s relation between magnetic and gravity elements, written 


(see Eq. 4) 
X — iZ = Ie'¥-(K — iG)/fo = ff (x—p)*dS (7) 


indicates that in this case, instead of o, the angle of inclination y and 
the magnitude J of the magnetization vector are also unknowns of the 
problem. We have the following list of eight parameters 


I, p, a, B, 4, (8) 
whose numerical values must be determined by the interpretation. 
3. First Moments and the Position of the Points O and C 


First we must find the value of #, the distance 0,0 which fixes the 
position of the unknown point O on our profile, the origin O; being 
chosen arbitrarily. This can be done with the aid of the following gen- 
eral result: 

For any cylindrical body the center of gravity of the infinite area, 
bounded by the line of profile, which is taken as the X-axis, and the 
experimental curve Dg= Dg(x), is on the same vertical as the center of 
gravity C of the cross-section S of the cylindrical body. Therefore the 
ratio M,/M, of the first two moments of the function Dg(x) is equal 


to %. 
= Mi/Mo. 


Let p=£+if be a fixed complex number whose imaginary part ¢ is 
positive and x a complex variable. Integrating (p—«)~!- dx in the plane 
of this variable along a contour consisting of the interval (—Z, L) on 
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the x-axis and of the upper half of the circumference of a circle with 
radius L and center at the origin, we get 


L 
f (p — x)-!-dx = lim (op — x)7!-dx = — in. 


Observing that kA=kff4dS, we deduce from (3) the relation 
A 


Integrating (3) and (10) with respect to x from — © to +, we have: 


The imaginary parts of these equations give the expressions for the 
first two moments Mo, M,* 


Ma = f Dg-dx = rkA, f xDg-dx = (11) 


thus proving Eq. g. If the measurements give only the curve of the 
gradient G (survey with the aid of a gradiometer), so that d(Dg) =Gdx 
is known, then, integrating by parts, we have, using for the m-th 
moment of the gradient the notation H,,: 


i= f Gdx = 0, m= f xGdx = — rkA = — Mo, 


-f = — 2rkA& = — 2M, 


whence the value of z 
t= H2/ 2H. (12) 


~ A torsion balance gives two curves: G=G(x) and K=K(z). It is 
interesting to note that for the curvature K the first two moments 
vanish and its third moment is related to the depth 2 of the center of 


* The integral M, is defined as limit for ZL of an integral over the interval 
(-L, L). 
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gravity C. Eq. 5 proves that the third moment of K can be defined 
only as the moment of the difference (K —kAx~*): This difference can 
be considered as a known quantity since the product «#7(K—kAx~”), 
whose integral is the value of the third moment, is equal to x*K+ H,/r, 
H being equal to —7kA. To compute the third moment Le of K(x) 
we first integrate Eq. 10 with respect to x over the interval (— ~, «) 
and take the real parts which gives 


f + kA/x)dx = — rk Real {i J fst - 
A 


Next integrating this integral by parts we transform it into Le, whence 
Ly [x2K + Hi/a|dx = — = 2H (13) 


and the formula for the depth 2 of the center C is 
L2/2H,. (14) 


In general a torsion balance survey gives much more information 
than a gravimeter or a gradiometer survey and the result of Eq. 14 is 
a good illustration of this. 

Let us now study the case of a magnetic survey. From Eq. 7 we 


have 
kX = 2I1(K cosy +G sin y) 


(15) 
kZ = 21G cosy — K siny). 
Using Eq. 5, we have for large | | : 
X~2IA cosy-x? and Z~ — siny- a. (16) 

Multiplying (15) by dx and x- dx and integrating over the interval 
(— ©, ©) we express the first two moments Xo, X1 of X and Zp, Z; of 
Z in terms of moments [p> of K and H)=o, of G. 
Therefore and 


x-X-dx = X; = — siny; 
(17) 
f = — 2nIA cosy. 


To determine the third moments X2, Z2 we form and integrate the 
differences 
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= x?-(X — cos y/x?) 
2k-!-Ix?|(K — kA/x*)-cosy +G sin 
x?-(Z + sin y/x?) 
2k-!-Ix?|G-cos — (K — kA/x?) sin 


with the result: 


2I(L2 cosy + He sin y)/k 


Xe -{ [x2X + 


| (18) 
Z:= — = 2I1(H2 cosy — Lz sin p)/k. 
Substituting in Eq. 18 the values from Eqs. 12 and 13 of Hz and Le 
we find 
Xo = — 4nIA(&-sin + Z-cos yp) (19) 
I 
Z2 = — 4nIA(#-cos — 2-sin y). 


The four moments X,, Xe, Z1, Zz can be used to find four among eight 
parameters listed in (8). These four parameters are #, 2, y and JA. 
Their expressions are: 


2% = +217): (X1X2 + 21:22); 
22 = (x? + = X1Z2) (20) 
= arctan (X1/Z;) and = 


Consequently, these four moments Xi, Xe, Z1, Z2 not only fix the 
position of the point C but also give the inclination y of the magnetiza- 
tion and the value of the product 7A, which will give J if A is found 
from the interpretation process to be illustrated. 


4. Reduction to a Finite Interval of Integration (Profile) 


The theoretical formulas deduced in Section 3 cannot be applied 
as they are stated for the simple reason that the experimental curve of 
the observed quantity is known only in a finite interval (—L, L), if 
we denote the length of the profile by 2Z and choose the origin O; at 
its mid-point. Integrating the experimental curve and the curves 
derived from it with the aid of a planimeter over the finite interval 
(—L, L) instead of (— ©, ©), we obtain numbers which can be con- 
sidered as first approximate values of the moments, if the profile is 
sufficiently long and 2/L is very small. This is always the case if the 
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map is a good one since then the profile was prolonged sufficiently to 
get negligibly small values of the observed quantity at both ends. 
Taking into consideration the contribution of intervals (—«, —JL), 
(ZL, ©) we compute the reduction coefficients with the aid of which 
the numerical values of the moments can be calculated, integrating 
over the finite interval (—Z, L) only. 

Taking the origin of coordinates (x’, 2’) at the center of gravity C 
of the cross-section S and the axes parallel to x- and z-axes, we intro- 
duce the characteristic constants c,, defined by 


Avena = f f dS, 
A 


Thus, Co9=1, C19 =Co1=0. They depend on the form and orientation of 
S. The radius vector p in Eqs. 3, 4 and 7 is now p=&+2'+1(2+2’) 
and, expanding (1—p/x)-! and (1—p/x)~ in series and denoting the 
constants ff4p"dS by Bn, we have 


U,—iDg=—k> K—-iG= kd (n+ 
0 0 
The coefficients B, can be expressed in terms of #, 2, A, Cmq. For in- 


stance By=A, By=A(%+iz), For 
large | | we have the following approximate expressions 


(23) 


Z~2IAx-?-sin [1-+2(#+2-cot p)a-?], (25) 


where E+61:/2, — 2?+ 3(coo— Co2) (3¢21+6¢1%— Cos) /2, 
C3= 3(#—2?+c20—Coz), Cy=62-Ci. Using these expressions, we can 
evaluate with great precision the contributions of two infinite intervals 
(— 0, —L), (L, ©) to each among the nine moments Mo, Mi, Ho, 
Le, X1, Xe, Z1, Z2 used in Eqs. 9, 12, 14, and 20. Denoting their first ap- 
proximate values, obtained by integrating in the finite interval (—L, 
L), by the same letters with horizontal bars and neglecting the terms 
of relative orders (2/L)? and (#/L)?, we have—(26): Mo=Mo- 
= H2- Xi=X1- [r+4(@-# cot y)/rL]; 


\ 
a 
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[1+4(2+4-tan y)/L]. Therefore the rules expressed by Eqs. 
g and 12 become, 


& = [1 + 2(2 + 
and &= [1 + + (27) 


We can improve the first approximate values of #; obtained by 
using in Eqs. 9 and 12, instead of Mo, Mi, Hi, Hz the numbers Mo, M1, 
Hi, Hz, only after having achieved the first step of the successive ap- 
proximation and obtained first approximate values for all parameters, 
since they are needed for the computation of the reduction factors in 
(26). Thus, beginning with #1, we calculate the first values for all 
other parameters (see Sections 5 and 6 below) and we improve the 
first value #, with the aid of Eq. 27. With the aid of the second ap- 
proximation £2 we start the computation of the second approximations 
for all other parameters and we continue this successive approxima- 
tion process until stabilization of sequences for all parameters occurs. 
The number of approximation steps required depends on the value of 
the ratio 2/L and it decreases if the length 2Z of the measured profile 
increases. 

We turn now to the rule expressed by Eq. 14 and compute the 
integral of x*K+H,/m over the intervals (— ©, —L), (L, ©). Denot- 
ing the integral of x°K+H/m over (—L, L) by Lz and observing that 
we have 


ie = 4L 2/1? + 62 — — 2) /aL. (28) 


With the aid of Eq. 28 we deduce from Eq. 14 the important 
practical rule for the depth 2 of the center of gravity of any cross- 
section 

= aLe(1 + 62/L)/H, — d#*/L, (29) 
where the coefficients a, 6B, are defined by 2a(m7?—4)=7", 
(1? —4)-A=3, 4) B= 16—T?+ C20) 

For the interpretation of a magnetic survey, given the curves 
X = X(x) and Z=Z(x) in the interval (— Z, LZ), we transform the rules 
of Eq. 20 with the aid of Eq. 26 into 


tan y= (1 — w#/L)Xi/Z,; IA=(1 + 42/aL)(X1? + (30) 


where 7X121:4=4(X1?+Z;’). We also transform Eq. 20 into a form 
which admits for #, 2 a numerical computation with the aid of inte- 
grals over (—L, L). Multiplying Eq. 7 by x and integrating, we ex- 
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press X1:—iZ, in terms of Z; and Hj. Using the fact that Z:=0 and 
Eq. 23, we find that Zi=4H,-%-/aL. Therefore 
= 2le¥- Hy, = Hye¥ (31) 


Integrating Eq. 7 multiplied by «*, we find for the reduced mo- 
ments Xe, Z, their expressions in terms of Z:, He. These moments 
being defined by 


L L 
= (x?-X + = (x?-Z 


we have the relation 
iZ>) Hy: = iH2/Hy + |}. (32) 
The relation (31) gives 2Ze¥-H,=ik(X1:—iZ,)[1—4i(#+i2)/rL]. 
Solving Eqs. 27 and 29 for the ratios Z./Hi, H2/Hi and substi- 
tuting their expressions in Eq. 32, we get in accordance with 
H,=Hy(1—42/7L) that 
— (4 + + — — + €1:/2)/L] 
= a(Xe = k/2Te¥- 3. 
Using the notation xX (X2+Z,’) -N. 2) 
= (X’+Z,’)-N., we can write the results in the final form 
Z = a-N,-(1 + B2/L) — Bx?/L; 
%=aN,: [x + 2(Bz Ne11/2)/L]; 


which is well adapted to numerical computation. 


(33) 


5. Moment-functions and Moments of Squares 


The interval of integration being (— ©, and the order s of 
moments being a continuous variable, we define the three moment- 
functions M, L, H related to Dg, K, G respectively by the integrals 


M(s) = |*-Dg(x)dx; H(s) = fi 
34 
L(s) =f | a 


The orders of Dg, K, G at infinity are | «| Therefore 
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the variable s must be limited by the conditions: —1<s<1 in M(s), 
—2<s<oin L(s) and —2<s<z1 in H(s). To simplify the calculations 
we suppose that the point O, already located by using the rules given 
in Section 4, is chosen as the origin of coordinate x on the profile. 
We begin the study of moment-functions (34) by deducing their 
general expressions which hold for any form of the cross-section S. In 
this deduction we make use of the following lemma from analytic func- 
tions theory. 
LEMMA?*: “Let F(z) bea rational function such that lim| 2: F(2)| =0 
both when z—o0 and when z>. 
(i) If F(z) has no poles on the positive part of the real axis 
~ and if s’ denotes the sum of residues of (—z)*~!- F(z) at all 
its poles, then 


x*-!.F(x)dx = ms’cosec sm, (354) 
0 
the argument of the complex variable (—2z) involved in s’ 
verifying the condition | arg(—z)| <z. 
(ii) If F(z) has no poles on the negative part of the real axis 


and if s’’ denotes the sum of residues of z*-! F(z) at all 
its poles, then 


| a |*1-F(x)da = — ms’’-cosec sm, (35.) 


the argument of z verifying the condition | arg(z)| <7.” 
From (3) we deduce 


= ik f 6) 


where p=£—i¢ is conjugate to p=£+i¢. Integrating the product 
| «| sl. Dg, we get from Eq. 36 


2M(s — 1) 


Il 


2f | 


is f J (p — p)dS 


* Whittaker & Watson. Modern Analysis (1927), Ch. VI, 6, 24, p. 117. 


4 
. 
j 


478 ERVAND G. KOGBETLIANTZ 


with F(x) If o<s<2, x*- F(x) satisfies the condi- 
tion of the Lemma. Integrating over the interval (0, ©) and using 
Eq. 35, we have 


in f [ce — p)(x — = ix(r + #)-cosec sz, 
0 


where r and # are the residues of (—z)*-!-(z—p)-!-(z—p)— at the 
poles p and 

For all points inside the cross-section S the radius vector p verifies 
the inequalities o<argp<7, while for its conjugate p we have 
—m <argp<o. Therefore, in calculating r and 7, we must choose for 
2=p, —2=—p=ep, but for z=p, —2=—p=e'*-p and, denoting 
the real part of a complex number w by R(w), we have: 


— p)(r + 7) = 


whence 


f x*-!.Dg(x)-dx = br-cosee ie ff (37) 
0 A 


In the interval (— ©, o) we apply Eq. 35, 


0 
f | a |*-1-F(x)dx = — x(r’ + #) cosec sz, 


where r’ and 7’ are the residues of 2*-!-(z—p)-!-(z—f)—! at p, p under 
the condition | arg(z)| <a. Therefore, r’+-7 = (p— (p*1— p*) and 


i(p — — | dx = -cosec sx; 


fi x = kr cosec Sf (38) 


Adding Eqs. 37 and 38, we have 


| a |*-1-De(x)dx = kr-cosec (sr/2)R ff 
A 


| 
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and, changing s into s+1: 


M(s) -f | x |*-Dg(x)-dx 


= km sec Sf (39) 


To find H(s) and L(s) we use Eq. 4, which gives: 
L(s) — i-H(s) = ef fas f | (a — p)-?-da. 
A 
Applying the Lemma to the function 2’: (z— p)—*, we have, if—2<s<o: 


f | — p)-?-da = — (s + 1)m-csc pt 
and therefore 


L(s) = | Kdx | 


= — (s+ csc (sx/2)R Sf. pris | 
(—2<s<o) 


f | |*+1-Gdx 


(s + 1) km csc (sr/2) E | 
(—2<s<1) 


+. (40) 
H(s) 


Proved for —2<s<o, the expression for H(s) holds, by analytic con- 
tinuation, if —2<s<1. It is interesting to note that if -1<s<o 


L(s) = — (s + 1) cotan (sx/2)-M(s) (41) 
that is 


fi a|*-Kdx = s-tan (sr/2) |*-1-Dg(x) (42) 


(fo <s< 2) 


| 
= 
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and this relation holds, as all relations of this section, for any form of 
the cross-section S. If the graph of the curve G=G(x) is symmetrical 
with respect to the origin O, the function H(s) vanishes. In this case 
we can use instead of H(s) the integral of the product | «| *+1- G(x) 
in the interval (— ©, o). In general, integrating by parts in the 
intervals (— ©, 0), (0, ©) we get 


0 
| |*H4Gdx -{ wHG-dx = (s+ Df | x 
0 

= (s+ 1):M(s) 


and, since the sum of the two integrals in the left hand memberis 


H(s), 
0 
of | 


‘ f = te 


*. Dg(x)dx 


H(s) + (s + 1)M(s); 
(43) 


This general result gives in the symmetric case when G(—x)+G(x) =o 
and H(s)=o 


| =— =(s+41)-M(s)/2. (44) 


In each particular case when the cross-section S is prescribed the 


integral 
ff p*-dS 
A 


can be calculated and expressed in terms of the parameters of the 
particular problem considered. Therefore each particular value of s 
gives us an equation in parameters of the problem and, choosing as 
many particular values of s as there are parameters characterizing 
our problem, we form a system of equations for them. 

We add now the expressions of the moment functions X(s) and 
Z(s) defined for —2<s<o by , 


and which are used in interpreting the magnetic anomaly map. 
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Multiplying Eq. 15 by | «| *—! and integrating, we have the expressions 
kX(s) = 21 { H(s) sin + L(s) cos v} 

kZ(s) = 21 { H(s) cosy — L(s) sin v} 

where the inclination angle y can be considered as known (see Eq. 20). 
It is important also to know how to compute the moments of 


squares Dg’, G?, X*, since they are very useful in our interpretation 
method. Defining the first four moments of Dg?(x) by 


(45) 


Dn -{ a”-Dg*(x)-dx (m = 0, 1, 2, 3) 


we use the expressions of the products «”Dg?(x) in terms of quadruple 
integrals extended twice over the area of the cross-section S. Integrat- 
ing them with respect to x in the interval (— ©, ©), observing that 


— p1)(% — p2)|-!-dx = 0, 


— p1)(% — pe) = — 


and using the fact that for e=o, 1 
f — Dg?|dx = f xte. [(U, + kA/x)? — Dg?|dx = 0, 


we get 


Al A2 


where Py= P\= pit pe, 2p1p2, P,= (p1+Be) pipe. Here the in- 

tegrals of x°Dg?(x), x3(U,+kA/x)? and xU,? over the interval (— ©, 

are defined as limits for L— of integrals over the interval (—Z, L). 
For m=o0, I, 2 we have also 


f — = f x”-K-G-dx =o 


and 


x™-G?-dx = f Pn(p1, po) (p1 po) deo, 
ri Jr2 


(m = 0,1, 2) . 
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where the double line-integrals are extended twice over the boundary 
curve (T) of the cross-section S, the function Pm(p1, p2) being the same 
as above. The relations (15) give now 


f -f x"Z2dx = 4k-?-J?- f x™-G?-dx (m = 0, I, 2). 


We add another general and very important expression of Dg as 
the Fourier transform of a function F(#) which characterizes the 
cross-section S 


Dg(x) = F(t) db, (47) 
Here . 
F(t) = ff signum t)-t. gS 
A 


with signum o=0, signum t= | t| /t for to. To prove (47) it is suffi- 
cient to substitute in (3) the following expressions for i(p—x) 


0 
i-(p — = f = it. dt 
0 


and to change the order of integration over ‘ and over A, which is 
possible, the integral being absolutely convergent since {2 min [>0. 
Applying the transform theory we get the value of the integral 
[(t)=27)?- F(t) to be 


= fe = rh ltl. dS, (48) 


The function J(#) or its real and imaginary parts can be used for the 
interpretation instead of the moment function M(s). It is easy to de- 
duce the expression of M(s) from Eq. 48 by applying to it the fractional 
integration with respect to ¢ of order 1—s, o<s<1, and observing that 
fort2o 
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Therefore, for ‘>o0 and o<s<1 we get 


f = ax f pl. dS. 
20 


Differentiating this result with respect to ¢ and observing that 
after the differentiation both members are absolutely convergent 
integrals for | s| <1 we have by analytic continuation for —1<s<1 
and ¢>o 


f x*-Dg(x)-e**!-dx = ax f p*-et-dS (| s| <3, 
A 


For to we obtain the result 


= (| s| <1) (50) 


from which it is easy to deduce the expression (39) for M(s), observing 
that in the range (— ©, 0), x= | «| e*, 

The Parseval theorem for the transforms gives us also an important 
relation holding for (¢>o0) 


dx 
(51) 


from which we can easily deduce another proof of the results expressed 
by Eq. 46. 
For G we have an expression similar to Eq. 47, 


Gla) = g(a, (52) 


where g(t) is expressed by a line integral extended over the boundary 
curve (IT) of the cross-section S, 


g(t) = 
r 
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Finally we point out the expression 


K(x) = 


I(t) = g(t) -signum ¢, 


The formulae (52), (53) give corresponding expressions for X(x), Z(x), 
and the transform theory applied to all these expressions gives the 
general rules for computing integrals which are very useful in the in- 
terpretation. 


6. Interpretation of Geophysical Measurements Made 
Above a Symmetric Anticline 


Gravity map. Among the seven parameters %, p, k, a, B, y, 6 the 
first, %, is given by Eq. 27, Section 4, which in the symmetric case 
(¢11=0) reduces to 


x= (z + 28/mL)-Mi/Mo. (54) 


Using M,/M, as its first value we can plot the first position O, of the 
point O. The value 7kA of Mp is now 1-k(A1—Az), A; and Ag being 
the areas of the sections C’/ACDBD’C’ and C’CDD’C’ (Fig. 2) re- 
spectively, since to the last one we have ascribed the negative density 
contrast —o. Therefore My=2mfo(A—1)T and for o we have the rule 


o = k/2f = + — 1) IT. (55) 


Expressing p, T, g, \ as functions of 1, a, 8, y, 6 we obtain T=I-sin 
(y—4)-csc y-csc 6, p=1-seca-csc y-cos (a+), g=l-tan a, \=1+sec 
B-sin y-sin 6-sin (a—£)-csc (y—6). We can consider the inter- 
pretation problem as solved if the values of five parameters /, a, B, y, 6 
are found. Calculating the explicit expression of the function M(s) 
for —1<s<1, we will use the abbreviations y/(s), f(s). These functions 
of four parameters a, 8, y, 6 and of the number s are defined by 
¥(s)=(s; a, 

= sin a-csc*ty- [cos*t* (a + y)-sect* a — cos (a + 2y + sy) | 
and f(s) =f(s; a, B, 7, 6) 
= on sec (sx/2):[(x + + 3) [W(s; — Ws; B, 


| 
where 
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With these notations we have the result: 


= f(s), (56) 


In particular, ¥(0; a, y)=tan a, f(o)=7(tan a—tan y), and 
y(1; a, y)=tan a-(3 ctn y—tan a). Therefore 


My = M(o) = rki?-sin (a — B)-sec a-sec B. (57) 


The depth Z of the center of gravity is a function of J, a, B, y, 6 and, 
defining the function ¢=¢(a, B, 5) by 


¢=¢(a, B, 7, 6) =371-cos a-cos B-csc (a—8) a, B, 6) ], 
we have for 2 the rule Z=/-¢, that is 


Now every numerical value of s~o between —1 and 1 yields an 
equation in five unknowns /, a, B, y, 6, since the numerical value of ° 
M(s) can be computed with the aid of a planimeter, using the experi- 
mental curve Dg= Dg(x). This curve being known only over the inter- 
val (—L, L), a correction is necessary. The difference M(s)—M(s), 
where M(s) is the moment function reduced to the finite integration 
interval (—ZL, L), is proportional to (//L)!~*. Therefore the negative 
values of s are better suited for interpretation, the correction needed 
being much smaller. Once and for all we fix the five values of s defined 
by s=—m/8, m=3, 4, 5, 6, 7. We then have five equations of the 
general type M(s)/Mo=I*-f(s)/f(o), and eliminating /, we can re- 
duce them to a system of four equations in four unknowns a, 8, 7, 6. 
Each of them is of the form E(a, 8, y, 5) =”, where the known number 
n denotes M’(s)-M-*(r)-M,*"=n, the function E being defined by 
E=f"(s)-f-*(r) -f*-*(0). For j=1, 2, 3, 4. we can define E; by fixing the 
_ values of s, r and letting s=s,;, r=r,;. For instance, we choose once and 
for all so=ri=—.75; Ss=re=—.625; 53=rs=—.5 and 
14= —.375. The tabulation of the functions £; is facilitated by the 
fact that in general 


sin (a — — 8, 
[v(r; a, — ¥(r; B, 8) ]* 


(59) 


E = 


cos a cos B 


485 

Z= p+ (A+ 2)T/3 — (4 — A)g/3(4 — 1) =1-6(a, B, 7, 8). (58) 
= | 
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the constant C,, being defined by 
(sr/2) cos (r/2) |*. 


Therefore it is advisable to plot the graphs of five families ¥(s; x, y) 
=constant, for s= —m/8, m=3, 4, 5, 6, 7. Solving the system of four 
equations in four unknowns 


E;(a, B, 6) = Nj, (j = I, 2, 3; 4); (60) 


we find the values of a, 8, y, 6. Any moment-function M(s) now gives 
1 since in general 


= a, B, 5) “M(s)/Mo- f(s; Qa, B, Y 6), 


and o is given by Eq. 55. 

Now, the curve Dg=Dg(x) being known only over the range 
(—L, L), the numbers in Eq. 60 defined by n= M"(s)- M-*(r) - 
must be replaced in the first approximation by #, defined as i= M"(s) - 
M-*(r)-M,*. Denoting the solution of the system E;=%; by ai, 1, 
1, 51 we deduce with their aid /, and the first corresponding value 2;, 
_ using the rule expressed by Eq. 58. To improve these first values by 
means of the process of successive approximations we must evaluate 
the difference M(s)—M(s). For | «| =L, the profile being sufficiently 
long on both sides of the axis of the anomaly (the anticlinal axis), for 
the gravitational action of the anticline can be substituted that of a 
horizontal line source at the depth Z, so that for | «| = Z we can use the 
approximate formula 


Dg(x) ~ + 2) (| «| = ZL). (61) 
Therefore, neglecting the terms of relative order (2/L)?, 
M(s) — M(s) = 2Mo-2L*"/x(x — s) 
and for M(—s) the reduction formula is 
M(— s) = M(— s) + + (s> 0), (62) 


where My = M)(1—22/rL)—. With the aid of a1, 81, 1, 61 we can 
compute, using (Eq. 62), the improved values of M(—s), and deduce 
the improved second approximations for the second members 1; of 
Eqs. 60, as well as improve the position of the point O. Solving anew 
the system (60) with new values of second members, we deduce the 
set of second approximations for our five parameters /, a, B, y, 6 which, 
in turn, used in Eq. 62, gives us better values for the second mem- 
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bers m; in (60). This process is to be continued until stabilization of 
the values of the four parameters a, 8, y, 6 occurs. Substituting the 
final values into the expression (56) for the moment function M(s) for 
any fixed value of s other than the five values 8s= —m, m=3, 4, +--+ 7, 
and comparing the result thus obtained with the experimental value 
of the same M(s) given by Eq. 62 where M(s) is computed by means 
of a planimeter, we have a control on our results. 

The solution of the system (60) can be facilitated, if a set of charts 
of four families E,(a, B; x, y)=constant, E2(a, B; x, y)=constant, 
E;(a, B; x, y)=constant, E4(a, 8; x, y) =constant is plotted for differ- 
ent values of a and B<a. To four experimental data m1, m2, m3, 4 Cor- 
respond four curves (60) on each chart for a=5°, 10°, 15°: - + upto. 
a=85°; B=0°, 5°, - - --<a. These four curves form on every chart a 
curvilinear quadrilateral, whose area considered as a function of a, B 
must vanish for the values a=&, B=8, roots of the equations (60); 
the coordinates #, j of the intersection point (£, 7) common to these 
four curves being two other roots 7 =, 5=7 of Eqs. 60. Consequently, 
in practice, we have two such charts 6’) and (a@”’, 6”) with 
a!’ =a'+5°, =6'+5°, such that the above quadrilaterals have the 
smallest areas on these two charts. It is plain that a’<a@<a”’, 
6’ <B<B" and 7, 5 are given approximately by the coordinates of the 
center of gravity of the quadrilateral. More precise results are there- 
after obtained, plotting the set of charts for a=a’+1°, a’ +2°, a’+3°, 
a’+4° and B=6’+1°, B’+2°, B’+3°, B’+4° and studying the same 
quadrilaterals on these detailed charts. 

Another interpretation can be based on the use of the function J(?) 
(Eq. 48) for different values of the positive constant ¢ which has the 
dimensions of (/ength)—}. Its explicit expression in our symmetric case 
is 


862) cos xt-dx 
= I(t) = [@(a, It) — It)] (63) 


~ 


where 


u) = sin a: [cos — cos (a — u)]. (64) 


To find the four unknowns a, 6, 1, T we form the system of four 
equations 


Q;(a, B, L, T) = Mi, (j = I, 2, 3, 4) (65) 
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where m= 2'-*-]*(t)-I-*(s)-Mo** are known numbers, the functions 
Q; being defined by 


Q(a, B, 1, T) = [¥(e)]*- [W(s) (66) 


where 
(lt)?-sin (a — B)- V(t) = [®(a, lt) — It)] cos a-cos B. 


Choosing once and for all the eight constants ¢;, s; (7=1, 2, 3, 4) 
and using Eqs. 64 and 66 we can tabulate the four functions Q;. On 
the other hand the numbers m; are computed with the aid of the ex- 
perimental curve Dg=Dg(x). Solving the system (65) by successive 
approximation we achieve the interpretation. This solution proceeds 
exactly as that of Eq. 60 where instead of Eq. 62 we use the reduction 
formula 


I(t) = I(t) — 2Ma[sin (tL) — 2 cos (tL) /tL | (67) 


where as always J(t) denotes the integral of Dg(x) cos xt- dt over the in- 
terval (—L, L). Knowing the values of a, B, 1, T, we deduce g=/- tan a, 
\=1+(q—/ tan B)/T and the value of is given by Eq. 63 using any 
I(¢;). The parameters Z and a are found with the aid of Eqs. 58 and 55. 

Thus, each of two independent methods, based on M(s) and I(#) 
respectively, yields the numerical values of all parameters involved, 
locating the apex A, determining the slope a, thickness T and AT, 
depth of base p+, width 2] of the symmetric anticline, and giving the 
density contrast ¢, deduced from the geophysical measurements. 


Gradient Map 


A gradiometer (torsion-balance with three masses) is a much more 
precise apparatus than the ordinary torsion-balance. While it gives 
only U,, and Uyz, being insensitive to the curvature components U4, 
and 2U,,, this is compensated for by the fact that the gradient compo- 
nents are given with the remarkable precision of a half-Eétvés unit 
(.5X107® cgs.). 

With this apparatus a map of the gradient is obtained and we pro- 
ceed to study its interpretation in our symmetric case. Since H(s) =o 
we use instead of the moment-function M(s) the integral h(s) defined 
by 

h(s) = f | | G(x) | (67) 
(—-3<s<1) 


The moment replaces My and —Hi1=Mp. 
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Nothing is changed in the general procedure of interpretation ex- 
cept for the tables of functions used, charts and some formulas. Thus, 
for instance, Eqs. 60 and 62 are replaced respectively by 


Ej‘(a, B, = (j = I, 2, 3, 4) (60) 
h(— s) = h(— s) + + (62’) 


where h(s) denoting the 
integral of Eq. 67 but taken over the finite interval (—Z, L) rather 
than over (— ©, «). The successive approximation process is applied 
to six parameters a, B, y, 6,1, #. For the last one the formula of Eq. 28, 
Section 4, is to be used. The moment Hz involved in it is always de- 
fined as an integral over (—L, L) instead of over (— ©, ©). We ob- 
serve that s can now take all values between —3 and 1. For s=—2 
we must use the fact that 


O(a, 7) = lim [(2 + s)-1-¥(s)] 


s—2 


= sin a cos a[y tan a + log cos (a + y) — log cos a] 
which gives for 4(— 2) the special formula 
h(— 2) = — k lim f(s) = 2rk[0(B, 6) — O(a, )]. 


The arguments cos (a++), cos (8+4) of logarithms are always posi- 
tive since a+7<go0°, B+6<g0°. To diminish the importance of the 
correction for the range of the infinite integrals (— «©, —L), (LZ, ~) 
it is advisable to use, as values s;, numbers in the interval —3<s;S 
— 2. The function 


H(t) = fe@-sin xt-dx (68) 


can also be used for the interpretation of the gradient map since in- 
tegrating by parts we have 


H(t) = tf“ Dg(2)-cos at-dx = — t.I(t). 


The system (65) is to be solved with the only difference that its 
second members m; must be computed by the rule 


m = s‘-t-*-H*(t)-H-*(s), 
‘the reduction formula for H() being 
H(t) = H(t) + 2hoa[sin (tL) — 4 cos (tL) /tL]/nL*. 


| 
‘ 
/ 
AS 
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Torsion Balance Map 


If we have to interpret the two curves G=G(x) and K= K(x), the 
first approximate values for # and 2 are given by the results of Eqs. 28 
and 30 of Section 4. It is interesting to observe that a control of the 
value of # is obtained by using the moments of the product GK. De- 
fining them by 


= x" -G(x)-K(x)-dx 
we find easily that the first three moments IIo, II, IIz are zero, but 
= — rk?A?, Il, = — 
which gives the relation 


The reduction to the finite interval (—Z, L) in our symmetric case is 


L L 
Tl; -{ x°GKdx & — 82/mL), x*KGdx = Il, 


L L 


which yields the important formula 
& = (1 + (69) 


The five parameters /, a, 8, y, 6 can be found with the aid of the ex- 
perimental curve K= K(x) and this controls the set of values of the 
same parameters obtained by interpreting the curve G=G(z). 

For —1<s<o we have in general 


L(s) = fi K(x)dax = — (s+ 1) cot (sr/2)-M(s). (70) 
(s ¥ — 1) 


In our symmetric case substituting for M(s) its expression (56), we 
obtain for L(s) the expression which holds by analytical continuation 
for —2<s<o 


L(s) = + s)!-k-esc 8, 8) — a, (71) 
(—2<s<o). 
Eqs. 60 are now replaced by 


Ej'"(a, B, 5) = n;"", (j 4) 
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where n”’ = L’(s)- L~*(r) -o** and the function E” can be defined as E 
by Eq. 59 with the only difference being that the constant C,, in 
(59) is to be replaced by C,,”’, this coefficient being defined by 
2*-(2+5s)"- | sin (sm/2)| = 2"-(2+1r)*-| sin (rm/2) | The reduction 
formula for L(—s) is 


L(— s) = L(— s) + 


We add two control relations obtained comparing the moments 
Hi, He, Tz and Ily. Reducing them to the interval (—Z, L), we have 


It is also possible to use for the interpretation the function K(é) 
defined by the integral 


K(t) = cos (xt) -dx. 


In our symmetric case the integral of the product K(x) -sin (x#) is 
zero and therefore the transform theory shows (see Eq. 53) that 
K(t) =(27)?-1(t). For t20 we have /(t)=g(#) and also K(t)=A(2), 
since H(t) = (2m). g(#). The relation 


K(2)-cos (xt)-dx = (xt):dx (#20) 


which holds only in a symmetric case for t20, proves that all that has 
been said on the use of H(t) = —¢-I(¢) can be repeated for the K-map, 
using K(t) instead of H(t). The reduction formula for K(é) is 


K(t) = K(#) — 2ho- [sin — 2 cos (tL) /tL /atL?. 


Magnetic Maps 


Eqs. 30 and 33 give the rules for computing the parameters £, 2, y 
and JA. Their first values are obtained by omitting in Eq. 33 in the 
expression for 8 as defined following Eq. 29, coz and C29 (¢11=0 in the 
symmetric case). To find the other five parameters we will use instead 
of M(s), k(s) or L(s) the moment functions X(s), Z(s) defined in Sec- 
tion 5. In our case H(s) =o and Eq. 45 becomes, using the value of Mo 
given by Eq. 57: 


X(s) = cosp-L(s)/Mo, Z(s) = — 2nIA-sin y-L(s)/Mo. 


492 ERVAND G. KOGBETLIANTZ 


Since L(s) is given by (71) and 27JA and y are known, it is plain that 
each of these two functions X(s) and Z(s) can serve for the interpreta- 
tion in the same way that M(s) serves for a gravity map. According to 
Eq. 15 their reduction formulae are 

X(—s) = X(—s) + 4IA cos p/sL*, 

Z(— s) = Z(— s) — 4IA sin y/sL*. 

The interpretation can be based also on the use of any one of four 

functions X,(é), X.(¢), Z.(¢), Z.(¢) defined by the integrals 


X.(t) + iX,(1) = f “X(a)edx, Z(t) + iZ,() = "2 


In fact, with the aid of Eq. 15 it is easy to prove that in our sym- 
metric case, where the integrals over the interval (—~, ©) of 
G(x) cos (x#) and K(x) sin (xt) are zero, all these four functions reduce 


to H(t): 
k[Z.(t) — iZ.(t)] = &[X.(t) + iX,()] = H(t) (#2 0). 
Therefore any one of them can be used in the interpretation of a 


magnetic map in exactly the same way that H(t) is used for a gravita- 
tional map. We note the relation 


tan = X,(t)/X-(t) = — Z.(t)/Z,(t) (¢ 2 o). 


which holds for any positive value of ¢. The reduction formulae for our 
four functions are 
— (X,— X.)/cos = (Z, — Z,)/sin y 
= 4IA [sin (tL) — 2 cos (tL) /tL |/tL? 
(X, — X,)/sin = (Z, — Z,)/cos y 
= 421A [sin (tL) — 4 cos (tL) /tL |/L?. 
When the interpretation is achieved the values of all parameters 
are known and the areas A; and Az can be computed. Therefore the 
parameter A =A,—Az is given by the interpretation of a magnetic 


map and the intensity of magnetization J is obtained dividing the 
value of the parameter JA by that of the area A. 


CONCLUSION 


In the foregoing an outline is given of a new mathematical method 
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for the interpretation of the maps of magnetic and gravitational 
anomalies. These methods are general, and they are applicable to all 
possible cases of two-dimensional or centrally symmetric anomalies 
observed in any of the various types of gravitational and magnetic 
measurements. Based on the use of average values only, they yield 
quantitative interpretation, the determination of density contrast 
and magnitude and direction of the magnetization-vector included. 
They do not depend on the existence of well-pronounced maxima and 
they are applicable equally well to small anomalies. We confined our 
discussion to the particular case of a symmetric anticline in order to 
carry it out with enough detail to make the concrete application of our 
methods in this particular case as clear as possible. It is important to 
emphasize that, in general, the application of such quantitative 
methods presupposes an accurate correction for the regional anomaly, 
but cannot be expected to yield good results if this correction is made, 
as is customarily done, simply by arbitrarily smoothing the experi- 
mental curves. 

These methods, including the solution of the problem (A) (see In- 
troduction) which is not treated in this paper, have been applied by 
the author to the practical interpretation of field-problems of various 
kinds. 


\ 
a 
‘ ‘ 


PRECISE MEASUREMENT OF DEEP 
ELECTRICAL ANOMALIES* 


THOMAS S. WESTt AND CLARENCE C. BEACHAM}{ 


ABSTRACT 


A method of electrical prospecting has been developed in which the deleterious ef- 
fects of superficial inhomogeneities in resistivity are eliminated. An electrode configura- 
tion is used consisting of a current electrode at “infinity,” a pair of current electrodes 
with a constant separation of a few hundred feet which are used alternately, and a pair 
of potential electrodes with a constant separation of a few hundred feet, collinear with 
the pair of current electrodes, and at a distance from them which is varied to secure 
depth resolution. 

It is demonstrated that superficial influences are eliminated by forming the resistiv- 
ity increment, the ratio between apparent resistivity determined by use of the closer 
current electrode of the alternate pair and that determined by use of the farther current 
electrode of the alternate pair. It is also demonstrated that the area explored is below a 
point approximately half-way between the alternate pair of current electrodes. 

Resistivity increment curves, called Resistologs by the authors, can be correlated 
from station to station, and after they are correlated the lateral variation in the electri- 
cal properties of selected subsurface zones can be studied. As a result of the latter char- 
acteristic of the measurements, a possibility exists for the direct location of oil and gas, 
observational evidence for which is presented. 

The results of surveys made in Hart County, Kentucky, and in the Sam Fordyce, 
Seven Sisters, Pettus, Branyon, Dunlap, Oakville, and Rhode Ranch fields of Southwest 
Texas are displayed. 


INTRODUCTION 


This paper is the first publication of results of a research on elec- 
trical prospecting extending over a period of ten years. It became ap- 
parent early in this work that the effects of shallow inhomogeneities 
were the fundamental barrier to the practical extension of electrical 
methods to greater depths and smaller anomalies. This work may 
therefore be characterized as primarily a study of the nature of such 
inhomogeneities and the invention and experimental investigation of 
various schemes for eliminating the influence of such undesired 
anomalies from the final results. Quite a number of such procedures 
were investigated during this time. This discussion is confined to the 
procedure showing the greatest promise of practical application. 


* Presented at the twelfth annual meeting, Denver, Colorado, April 1942. 
1 Simmons Oil Company, Three Rivers, Texas. 
¢ Harbough-Potts Coal Company, Somerset, Ohio. 
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The method of investigation has been primarily experimental. 
The procedure was first investigated on a small scale using a water- 
filled test pond over which the electrode system was suspended by a 
framework. This permitted the insertion and removal of inhomoge- 
neities of various types without other change of conditions. After com- 
pleting small scale experiments this procedure was extensively investi- 
gated under field conditions on an experimental and semiexperimental 
basis. In field investigation it was endeavored to go beyond the usual 
procedure of accepting apparent agreement between experimental 
data and known or subsequently determined subsurface data as final 
proof of validity. Auxiliary tests were applied in each of such cases 
which were designed to determine whether results were explained by 
shallow inhomogenieties or other common sources of error. 

The necessity of relying almost entirely upon experimental investi- 
gation of electrical methods in general is well indicated by the com- 
plexity and diversified nature of the various phenomena known to be 
involved. The conduction of electric current in the earth is electrolytic 
in nature. Current flow results because of actual movement of the 
positive and negative ions into which formation electrolytes dissociate. 
Subsurface anomalies result because of the existence of the following 
electro-chemical phenomena, although the relative importance of 
each has never been ascertained.! 

(1) Polarization at the interface of formations containing different 
quantities of electrolyte per unit volume. 

(2) Polarization at the interface of formations containing electro- 
lytes of different concentrations. 

(3) Polarization at the interface of formations containing different 
electrolytes. 

(4) The presence of a medium which is permeable to some species 
' of ions and not to others. The relations existing when two solutions 
are separated by a membrane which is non-permeable to an ion, or 
ions, in one solution have been investigat@ extensively in connection 
with biological processes. Equilibrium occurs with the diffusible ions 
distributed in accord with Donnan’s Equilibria.? Because of the 
decrease in porosity of shales with depth*, anomalies from this source 
should increase in importance with depth. 


1 Creighton and Fink. Electrochemistry, Vol. 1. 
2 Bodansky. Physiological Chemistry, p. 22. 
3 Athy, L. F. A.A.P.G. Bulletin, Vol. 14, No. 8, 1930. 
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The presence of non-diffusible ions or charged particles appears to 
be a possible explanation of most of the potential difference variation 
of the ‘‘porosity” curve of electric well logs. For example, the anoma- 
lous flow of current resulting when two electrolytes of different con- 
centration are separated by a shale section so as to simulate the con- 
dition existing between a sand, shale, and a penetrating well, may be 
explained theoretically if the shale is impermeable to ions which are 
more concentrated in one electrolyte than the other, or to non-diffusi- 
ble colloidal material which is dissolved from the shale to a greater 
extent by one electrolyte than the other. Anomalous current flow 
under the above conditions and its relation to electrical logging has 
been studied experimentally by Mounce and Rust.* 


SMALL SCALE EXPERIMENTAL DATA 


Fig. 1 illustrates by small scale experiment the influence of shallow 
inhomogeneities on apparent resistivity data. The electrode system 
consists of two fixed current electrodes and a series of fixed potential 
electrodes along the extension of the line determined by the current 
electrodes. Connection was made to one pair of potential electrodes at 
a time, the separation between current and potential electrodes having 
been varied by making connection between different pairs of potential 
electrodes. A means was thus provided by which a given series of 
potential electrode positions could be accurately duplicated as to sepa- 
ration between potential and current electrodes. This electrode system 
was suspended over a water pond and in contact with the water. The 
water level was maintained constant and the pond shielded from air 
currents. A commutator by which potential and current leads were 
simultaneously reversed was used for overcoming polarization effects 
at electrodes. Potential electrodes were 6 cm apart in all cases and were 
moved by 6 cm increments. This arrangement provides a means for 
insertion and removal of inhomogeneities without other change of con- 
dition. Since the apparent Tesistivity curve for the water pond was an 
approximately straight horizontal line with all inhomogeneities re- 
moved, the influence of the inhomogeneity was readily apparent by 
inspection. 

It will be noted that the shallow inhomogeneity produces a resistiv- 
ity maximum which might be interpreted as resulting from a relatively 
deep high resistive subsurface layer, if all maxima are interpreted as 
resulting from horizontal subsurface layers. 


4 Mounce and Rust. A.I.M.E. Technical Publication No. 1626. 
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It will also be noted that the apparent resistivity curve obtained 
over the deeper inhomogeneity indicates that the resistivity maximum 
will be at an electrode separation somewhat greater than twice the 
depth. This curve has an inflection point at about 25 cm corresponding 
to a depth to the inhomogeneity of 24 cm. Such inflection points have 
been utilized as a basis for correlation and interpretation of field re- 
sults which will be subsequently described. 


PROCEDURE EMPLOYED FOR ELIMINATING THE 
INFLUENCE OF SHALLOW INHOMOGENEITIES 


The method for eliminating the influence of the shallow inhomo- 
geneities along the potential electrode traverse is illustrated by Fig. 2. 
A shows the type of apparent resistivity curve obtained if current 
electrodes remain stationary and potential electrodes are moved across 
a shallow inhomogeneity for otherwise homogeneous conditions. If the 
potential electrodes are now placed on or near the shallow inhomo- 
geneity and electrode separation varied by moving current electrode 
1, the potential electrodes and distant current electrode 3 remaining 
stationary, a curve such as shown at B, Fig. 2 is obtained. This curve 
is an approximately horizontal straight line except for relatively small 
electrode separation. It is evident that for electrode separations great 
enough to be on the flat portion of this curve the influence of the shal- 
low inhomogeneity is approximately independent of electrode separa- 
tion. This fact is utilized for eliminating the influence of shallow in- 
homogeneities from the data by the procedure shown at C, Fig. 2.5 
An additional current electrode 2 is placed in the vicinity of 1, prefer- 
ably along the line determined by 7 and the potential electrodes. For 
each position of the potential electrodes two measurements of poten- 
tial difference are made, one with current flow between 1 and 3 and 
the other with current flow between 2 and 3. These two potential dif- 
ference values along with corresponding current values and electrode 
separations provide data from which two apparent resistivity curves 
R; and Rz may be plotted. For electrode separations a of interest in 
deep prospecting switching from 1 to 2 results in varying the electrode 
separation by the distance between these electrodes along the flat 
portion of the curve of B of Fig. 1. Obviously no appreciable change in 
the influence of any inhomogeneity which might be present in the 
vicinity of the potential electrodes results from this operation. A ratio 


5 U.S. Patents $2,231,048 and #2,217,780. 
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is then taken between the two resistivity values for each position of 
the potential electrodes. This provides data for a resistivity increment 
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Fic. 2. Illustration of the incremental resistivity (Resistolog) method 
employed to eliminate the effect of superficial inhomogeneities. 


curve Ri/Re, which is approximately independent of shallow inhomo- 
geneities along the potential electrode traverse. This curve will not, 
however, be independent of subsurface inhomogeneities below elec- 
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trodes 1 and 2 since such inhomogeneities will produce an influence on 
the apparent resistivity curve which varies with electrode separation. 
Resistivity increment curve Ri/Re thus indicates the rate of variation 
of apparent resistivity with electrode separation and will therefore 
have maxima or minima corresponding to apparent resistivity curve 
inflection points. Inhomogeneities in the immediate vicinity of 1 
which differ from those in the vicinity of 2 will alter the resistivity 
increment curve, but the resulting deformation will be that of intro- 
ducing a smooth curve similar to that of B in Fig. 2, which is approxi- 
mately a straight horizontal line for large electrode separation. Ob- 
viously, adding such a curve to the resistivity increment curve will 
not alter the electrode separation at which maxima or minima occur, 
nor the relative magnitude of such features. For uniform conditions 
in the vicinity of z and 2 the resistivity increment curve Ri/Rg will be- 
come unity at the point the apparent resistivity curve has a maximum 
or minimum, will be less than unity if apparent resistivity values are 
increasing with electrode separation, and greater than unity if these 
values are decreasing. Interpretation of resistivity increment curves 
is based entirely on resistivity increment curve maxima and minima, 
as described above. These extrema are independent of inhomogeneities 
in the vicinity of both current and potential electrodes. For conven- 
ience in description apparent resistivity increment curves obtained in 
the above manner have been named Resistologs. 


SHALLOW RESISTIVITY DATA 


Application of this procedure to relatively shallow prospecting is 
illustrated by Fig. 3. These data were obtained in the Horse Cave Area 
of Hart County, Kentucky. Except for a surface layer of soil which 
varied from o’ to 10’ in thickness, the underlying beds were predomi- 
nantly high resistive limestones to a depth of at least two thousand 
feet. Electrode 3 was placed on the same side of z and 2 as the potential 
electrodes but was about twice as far from the potential electrodes as 
1 and 2. Curve Ri is the apparent resistivity curve obtained with cur- 
rent flow between electrodes z and 3 and Rz is a similar curve for flow 
of current between 2 and 3, except that the apparent resistivity values 
for Re are plotted against the electrode separation @ rather than the 
distance from potential to current electrode. Because curves R; and 
Re are practically coincident, curve Re has been offset by 20 ohms 
to facilitate comparison as to form. The two resistivity values 
plotted-against each value of a are therefore the two resistivity values 
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APPARENT RESISTIVITY DATA 
HORSE CAVE AREA 
HART COUNTY, KENTUCKY. 
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Fic. 3. Example from the Horse Cave area of Hart County, Kentucky, showing 
the elimination of the surface effects in the apparent resistivity by use of the resistivity 
increment. The apparent resistivity curve Re has been shifted down 20 ohms to facilitate 
comparison. 
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for each position of the potential electrodes. It is obvious that the 
various maxima and minima of these two individual apparent resistivity 
curves occur at a given position of the potential electrodes rather than 
at a given separation between potential and current electrodes and are 
therefore a result of shallow inhomogeneities along the potential elec- 
trode traverse. This conclusion is further supported by the resistivity 
increment curve (Resistolog) Ri/Rz, which was obtained by plotting 
the ratio between each of the two apparent resistivity values, as 
functions of a. Since each point on this curve is the ratio between two 
resistivity values which differ in separation between potential and 
current electrodes by 60 feet, any of the sharp maxima or minima of 
R,; and Rg occuring at a given separation between current and poten- 
tial electrodes (i.e. due to subsurface beds) would have produced a 
similarly sharp change in Ri/Re. 

Apparent resistivity curve maxima and minima such as those of Ri 
and Re have been frequently attributed to horizontal subsurface beds 
by previous investigators. Use of a greater separation between poten- 
tial electrodes would, of course, tend to decrease the relative magni- 
tude of the effect of these inhomogeneities. However, variation in ap- 
parent resistivity values by 20% along a potential electrode traverse 
600 feet long because of shallow inhomogeneities is quite common 
even under favorable conditions. 


MEASUREMENT PROCEDURE AND ELECTRODE 
ARRANGEMENT FOR DEEP PROSPECTING 


In the application of the above procedure to deeper prospecting the 
electrode arrangement shown in Fig. 4 has usually been employed. As 
illustrated the separation between potential electrodes and between 
electrodes 7 and 2 is 600 feet. The distance between electrodes 2 and 3 
is usually equal to or greater than the maximum depth to which pros- 
pecting is desired. The position of 3 is not critical. It is usually placed 
approximately collinear with z and 2 but various other positions, such 
as at right angles to the line of 1 and 2 or on the opposite side of the 
potential electrodes from r and 2, such as in Fig. 3, have been em- 
ployed. A practical difficulty as to the position of 3 is that it is difficult 
to maintain the degree of insulation necessary for precise work if the 
connecting current carrying line passes within 1000’ to 1500’ of the 
potential electrodes. The separation between potential electrodes is 
maintained constant for a given station. As previously described, 
electrode separation is varied by moving the potential electrodes, all 
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current electrodes remaining stationary for a given station. On initial 
attempts to apply this procedure to deep prospecting the distance 
between potential electrodes and electrodes 1 and 2 was varied by in- 
crements of 300feet to 600 feet on the assumption that deep subsurface 
formations would produce rather broad maxima and minima. Dis- 
covery that features actually obtained in the field were larger and 
sharper than originally assumed to be likely led to variation of 
electrode separation by smaller and smaller increments, the various 
increments being a result of successively dividing 600 by 2. The pres- 


3 2 


Approx. 6 Max. (3000. 6000')-» 600" 


Fic. 4. Schematic diagram of the circuit employed, showing also the 
electrode separations used for relatively deep prospecting. 


ent practice is to use an increment of 373 feet, although it is possible 
to obtain additional detail by use of an increment of 183 feet as in- 
dicated by the Resistologs of Figs. 14 and 18. 

It is apparent that the resistolog procedure minimizes errors in- 
troduced by line leakage and instrument calibration since most errors 
of this type will similarly effect the two resistivities between which a 
ratio is taken. To minimize the effect of variation of instrument cali- 
bration while observations are in progress the various potential dif- 
ference values which are taken for each position of the potential 
electrodes are alternated. The resulting two groups of potential values 
are then separately averaged for obtaining the two potential difference 
values which are used for calculation of the resistivity increment 
value. 

A rather detailed description of the procedure employed for poten- 
tial difference measurement will be given. This procedure is apparently 
critical for successful prospecting. Although the operation is quite 
tedious and requires considerable training in this particular type of 
measurement, skillful operation results in a precision of measurement 
which appears to be somewhat greater than otherwise attainable in 
the presence of ordinary earth current fluctuations. 
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Relatively small errors in measurement of potential difference may 
result in resistivity increment curves which do not positively correlate. 
A considerably higher degree of accuracy must therefore be main- 
tained than ordinarily justified, or secured, in electrical prospecting. 

Care must be exercised in laying out traverses in the presence of 
artificial conductors, particularly those extending from the vicinity 
of potential electrodes to the vicinity of the current electrodes. Wire 
fences, even where no connection with the earth exists, except through 
wooden posts, constitute the most troublesome source of error of this 
type, since errors of the same order of magnitude as anomalies due to 
subsurface features may result. In all cases in which wire fences 
cross or parallel the potential electrode traverse insulators are in- 
serted in each wire on both sides of this traverse unless gates or other 
breaks in continuity of the wire exist nearby. Where such insula- 
tion is impossible traverses are laid out in another direction. As pre- 
viously discussed, however, the magnitude of errors resulting from 
artificial conductors is greatly reduced with the Resistolog procedure 
relative to that obtained with usual electrical prospecting methods. 

At least ten separate measurements of potential difference are made 
for each position of the potential electrodes, that is, five for current 
flow between electrodes 1 and 3 and five for current flow between 2 
and 3. At least two independent sets of data are taken for each sta- 
tion, sets of data referring to two independent series of potential 
electrode positions. The usual practice is to first take data for obtain- 
ing an increment curve covering the desired range of movement of 
potential electrodes, and then to repeat these measurements one or 
more times depending upon the amount of discrepancy between the 
various sets of resistivity increment values. When sufficient data for 
at least two independent Resistologs, which check generally, are ob- 
tained, all data are averaged for the final curve. 

Potential difference measurements are made by the null method, 
utilizing a potentiometer having a slide wire by which potential differ- 
ence is varied by steps of one microvolt, and a total potentiometer 
range of 16 millivolts. A high sensitivity D’Arsonval type galvenom- 
eter having a voltage sensitivity of about 0.05 microvolt per milli- 
meter deflection on a scale at one meter is employed for balancing the 
potentiometer. The usual type of porous pot non-polarizable electrode 
is used, the resistance between these electrodes usually being about 
500 or 600 ohms. The galvanometer is therefore highly damped and 
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not responsive to rapid earth current fluctuations, since its critical 
damping resistance is 2400 ohms. A relay arrangement is provided by 
means of which the observer measuring potential difference can switch 
current to electrodes 1 and 3 (or 2 and 3) on and off by remote control, 
thereby enabling the truck containing potential difference equipment 
to be stationed near the potential electrodes and the truck containing 
current equipment near electrodes r and 2. This separation of current 
‘and potential difference equipment eliminates errors due to leakage 
between them or between current-carrying lines and potential equip- 
ment. 

. Observation of potential difference is made by simultaneously 
switching, on and off, potentiometer working current and the current 
to the earth, and adjusting the potentiometer until the deflection of 
the galvanometer is not altered by this switching. Deflection of the 
galvanometer during this switching operation is arrested by short-cir- 
cuiting the coil. It is, of course, necessary to provide an auxiliary 
source of potential difference in the galvanometer-potential electrode 
circuit to compensate for potential differences existing between these 
electrodes when there is no flow of current between electrodes r and 3 
or 2 and 3. This procedure permits very precise measurement of po- 
tential difference in the presence of ordinary earth current fluctuations 
without the use of commutating devices, although the measurement 
is actually made during short quiet periods which usually occur. 
Fig. 4 shows a schematic diagram of the circuit employed. 

Observation of the galvanometer is approximately during the 
period of time between five and ten seconds after the circuit to the 
current electrodes is closed. The current flow is manually adjusted toa 
given value by means of a rheostat during the entire period of flow. 
The same current value is usually maintained for flow to each current 
electrode and for all positions of the potential electrodes throughout 
a given set of data. The source of current is storage batteries, the 
resultant non-fluctuating applied E.M.F. permitting a precise adjust- 
ment of current flow. Magnitude of current flow is usually between 4 
and 10 amperes for relatively deep prospecting in low resistive beds. 

By means of the above described procedure, five potential differ- 
ence values can practically always be obtained at a given potential 
electrode position which have a difference between the highest and 
lowest value not exceeding 5 microvolts. The overall accuracy of 
measurement of resistivity increment values attained in routine 
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prospecting has been such that an increment curve at a given station 
could be duplicated as to form within one quarter of one per cent 
(.25%) by an independent check. 

An average of about thirty minutes time is required for an expe- 
rienced operator to obtain one set of data for a given potential elec- 
trode position. Since at least two such sets of data are taken at least 
one hour is required to take data for each point on Resistologs. If 
potential electrodes are moved by the usual increment of 373 feet 
about 15 hours are required to cover a depth interval of 500 feet, ex- 
clusive of the time required for surveying, laying lines, etc. 

The distance between electrodes z and 2 and distances to potential 
electrodes are measured along a surveyed line by means of a sur- 
veyor’s chain, measurement being made along the horizontal. An 
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Fic. 5. System of overlapping potential electrode traverses which demonstrate 
that Resistolog features are not due to superficial inhomogeneities at the potential 
electrodes. Data from stations 36-40, Fig. 14. 


independent check was usually made of these measurements. Elec- 
trode 3 is usually located approximately by compass and automobile 
speedometer. 

The average of the elevation of electrodes z and 2 is used in all 
cases for designating the elevation of the station. The location of the 
station is designated as the midpoint of the line joining electrodes 7 
and 2, the location of the subsurface area prospected having been.ex- 
perimentally established as lying below this point by tests which will 
be subsequently described. 

A number of rather simple but positive tests may be applied under 
field conditions for proving whether features on Resistologs are of sub- 
surface origin. 

The most positive test of this type is that illustrated by Fig. 5. 
In this case potential electrodes were moved by increments of 150 
feet. By taking a series of stations along a given line at 150 foot inter- 
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vals (i.e., electrodes z and 2 both moved 150 feet) a series of Resisto- 
logs were obtained for which two of the potential electrode positions 
were identical with two of the potential electrode positions of the 
preceding station. However, the separation between potential and 
current electrodes for each of these two potential electrode positions 
differed by 150 feet. Any residual influence from shallow inhomoge- 
neities which might be present would result in features which occurred 
at a given potential electrode position, while any feature resulting from 
horizontal subsurface beds would consistently occur at a given elec- 
trode separation. It ‘will be noted that the Resistolog feature shown 
consistently occurred at the same separation between current and 
potential electrodes. The Resistolog feature shown in this case appears 
very pronounced because of the large scale to which the increment 
value is plotted. Such features would probably not be recognized even 
in the absence of shallow inhomogeneities with the accuracy of meas- 
urement ordinarily attained in electrical prospecting. 

Another rather positive test of validity of data under field condi- 
tions can be obtained by rotating the electrode system about the 
midpoint of the line joining electrodes z and 2 so as to obtain two 
or more independent Resistologs at the same station. As will be sub- 
sequently demonstrated the point midway between 1 and 2 overlies 
the subsurface area prospected. By turning the electrode system 
about this point any desired number of Resistologs can be obtained at 
a given station, in which all potential electrode traverses differ. 

Fig. 6 shows results of rotation of the electrode system sufficiently 
to place the potential electrode traverses 330 feet apart. Both resistiv- 
ity and resistivity increment curves are shown in this case, this station 
being number 38 of Fig. 14. Although the corresponding apparent 
resistivity curves obviously do not have definite similarity, the Resist- 
ologs are quite similar. 

Fig. 7 shows a case in which results were verified by rotating the 
electrode system through go degrees, the position of electrode 3 remain- 
ing stationary. The various sets of independent data which were aver- 
aged for obtaining the average, or final, Resistolog for each position of 
the electrodes, are shown. These two Resistologs are obviously similar 
as to features. The difference in slope probably results principally 
because electrode 3 was assumed to be at an infinite distance in mak- 
ing resistivity calculations. These curves are shown as computed and 
without rotation to the parallel position. The care necessary for ob- 
taining data of the required accuracy under field conditions is indi- 
cated by the fact that over 350 independent measurements of poten- 
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tial difference were employed in calculating the above two Resistologs, 
approximately 36 hours having been required for making these meas- 
urements. 

Fig. 8 shows experimental verification of the subsurface origin of 
features by rotating the electrode system through 180 degrees. This 
illustration shows the electrode arrangement employed, the various 


Potential Electrodes Collinear with !¢2-e-e— Potential Electrodes 330 off Line of 162 


RESISTIVITY Resistivity RESISTIVITY 
Ohms/Cm3 Increment 3 


(a) 


330 


é - 


3300 


Station *38 


Fic. 6. Repetition of Resistolog features when potential electrodes are offset 330 
feet, also illustrating the effect of surface inhomogeneities on apparent resistivity obser- 
vations. Data at station #38 of Fig. 14. 


sets of data which were averaged for obtaining the required accuracy 
of measurement, and curves indicating the accuracy of measurement. 

The average position of the potential electrodes with these 
electrodes southeast of the station of Fig. 8 was approximately 10,000 
feet from the average position of these electrodes when northwest 
of the station. The Pettus sand is estimated from a subsurface map 
of the area to be approximately 300 feet deeper under the potential 
electrodes when to the southeast of the station than when to the 
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northwest, and 150 feet deeper below the point halfway between cur- 
rent electrodes r and 2 than for the former position of potential 
electrodes. Since the two Resistologs obtained with the potential 
electrodes in these positions definitely correlate at the same electrode 


O oO 


Potential Electrodes North of Station —»4—Potential Electrodes West of Station 


Electrode Separation (a) 


° 
° 


Fic. 7. Repetition of Resistolog features when the electrode system (except the 
electrode at “infinity,’’ No. 3) is rotated 90° about a point midway between current 
electrodes 1 and 2. Data at station 41, Fig. 14. 


separation, it is evident that the subsurface area prospected in this 
case at least approximately underlies electrodes z and 2. More positive 
tests for determining the subsurface area prospected will be discussed 
later. The definite similarity of the above two independent Resistologs 
at station 2 is also rather positive proof that these features are not a re- 
sult of inhomogeneities along the potential electrode traverse. So accu- 
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type of prospecting shows the various characteristics usually re- 
garded as typical of resistivity prospecting. The features shown on the 
Resistologs are very large in comparison to those obtained in deep pros- 
pecting. It will be noted that the correlations check the variation in 
depth to the base of the Oakville sand at the three wells in the area 
within the limits of detail of Resistologs, or 182 feet. Results in other 
areas, however, indicate that cutting out of high or low resistive mem- 
bers produce a variation in the relation between depth of prospecting 
and electrode separation in such shallow prospecting. It is nevertheless 
possible to reliably detect abrupt changes in depth such as that pro- 
duced by a fault with considerable accuracy, if stations are sufficiently 
close together. Faults result in abrupt changes in indicated depth from 
one value to another, whereas stratigraphic changes usually result in 
a progressive variation of indicated depth, several intermediate 
values being obtained between the two extremes. As will be subse- 
quently indicated by field results, variation in relation between 
electrode separation and depth of prospecting is much less prevalent 
in deep prospecting and is probably of a different nature than that 
encountered as shallow depths. 


STRUCTURAL SURVEYS 
Branyon Field, Caldwell County, Texas 


Fig. 10 shows a traverse across a known fault in the Branyon Field 
of Caldwell County, Texas. The purpose of this survey was the deter- 
mination of the accuracy with which the Austin Chalk trace of a fault 
could be located. After locating the point along the traverse at which 
this fault was indicated by Resistolog correlations the electrode 
system was rotated through 180° about the point halfway between 
electrodes 1 and 2, and Resistologs at stations 6 and 7 verified, these 
being the stations between which the subsurface trace of the fault was 
indicated. Resistologs obtained with potential electrodes both East 
and West of the station (point half way between 1 and 2) are shown in 
Fig. 10 with those obtained with the potential electrodes to the West 
of the station plotted deeper to avoid the confusion of overlapping; 
the features shown actually occurred at the same electrode separa- 
tion as those from the potential electrodes East of the station. Since 
electric logs and straight hole surveys were not made on wells from 
which the subsurface trace of this fault was determined, the location 
indicated by Resistolog data appears to be within the limits of ac- 
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*rode system through 180° about the point midway between current electrodes 1 and 2, that the area prospected is below this point. 
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curacy with which it was determined by drilling. Resistolog corre- 
lations check the displacement of this fault in direction and to within 
1o feet in magnitude. Verification as to the subsurface origin of 
the features in this area, in addition to that obtained by rotation of 
the electrode system, is supplied by Resistologs at stations 6 and 7 
with potential electrodes to the East. The deeper five points of the 
Resistolog at station 6 were obtained with the potential electrodes 
in identically the same series of positions as were the five shallowest 


‘ points of station 7, the only change in electrode position being that 


electrodes 1 and 2 were 75 feet further West for station 6 than for 
station 7. It will be noted that the feature at the lower end of the Re- 
sistolog at station 6 corresponds to a smooth curve at station 7. Further- 
more, the validity of Resistologs at both stations 6 and 7 was proved 
by duplication of these Resistologs with the electrode system rotated 
through 180° about the point halfway between electrodes z and 2. 

A fault located by a similar traverse about one half mile west of 
the above was established by subsequent drilling to have an Austin 
chalk trace at least within 150 feet of that indicated by Resistolog 
correlations, while the magnitude of displacement was within 10 feet 
of the indicated value. 

The separation employed between electrodes z and 2 and between 
potential electrodes in this area was 300 feet. 


Dunlap Area, Guadalupe County, Texas 


Fig. 11 shows results along a traverse which was made for more 
accurately locating a fault, the position of which was approximately 
known from subsurface data. Oil production is obtained in this area 
from the Austin chalk, when this formation has sufficient permeability 
and porosity. Development in the area indicated that these required 
conditions were practically always present in or immediately adjacent 
to fault planes. Locations for D. G. Herring #1 and #2 Frank Appling 
and #1 C. F. Appling were made on the basis of the above survey. 
These wells encountered permeability and porosity in the Austin 
Chalk and were completed as oil wells. #2 Frank Appling encountered 
the top of the Austin Chalk 29 feet and 40 feet lower than did #1 
Frank Appling and #1 C. F. Appling, respectively, apparently verify- 
ing the 35 foot fault indicated by Resistolog correlations. The exist- 
ence of this fault had not been established by subsurface data prior to 
the Resistolog survey. The position of the second fault, having a dis- 
placement of 108 feet according to Resistolog correlations, has not 
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Fic. 11. Location of the Austin Chalk trace of a fault in the Dunlap area of Cioadahape 
County, Texas, by correlation of Resistologs. 


been accurately located by drilling, although, as previously stated, its 
existence in the approximate vicinity is indicated by subsurface data. 

It will be noted by reference to the Resistolog scale in Fig. 11 
that features on which correlations are based are relatively small. 
Correlations appear to be positive, however. It is also of interest to 
note that Resistologs at stations 4 and 5 which cover the subsurface 
interval in which the fault plane occurs, according to correlation of 
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Resistologs at stations 2 and 3, are not similar to other Resistologs in 
the area. The separation employed between current electrodes z and 2 
and between potential electrodes in this area was 300’. 


Pettus Area, Bee County, Texas 


Fig. 12 shows Resistologs along a traverse extending down dip 
from the Pettus field along the indicated line. Resistologs have been 
plotted against depth by use of a relation between depth and electrode 
separation to be shown in Fig. 18. The relation between depth of 
prospecting and electrode separation was not actually determined for 
this area, so it is not definitely known that the point on Resistologs 
correlated with the Pettus sand is actually that shown. Resistolog 
correlations conform to subsurface data, however. It is therefore ob- 
vious that Resistolog markers parallel and probably are near the 
Pettus sand. 

Correlations employed in making this section are shown at the 
right of Fig. 12. The depth to the Pettus sand along this traverse 
ranged from 3900 feet to 4100 feet. The range through which the 
electrode separation a was varied was between 4000 feet and 4900 
feet. The electrode arrangement was similar to that of Fig. 4. 

Because of the unexpectedly large features obtained on these 
Resistologs unusual effort was made to prove positively that they 
were not a result of superficial effects or errors in measurement. From 
three to five days were devoted to taking data at each station. More 
than 750 independent measurements of potential difference were used 
in calculating Resistologs at each of stations 1 and 3, a total of more 
than 1500 values for, the two stations. In case of each potential differ- 
ence measurement the apparently correct potentiometer setting was 
checked some five or ten times for definitely determining that no dis- 
cernible galvanometer deflection resulted, before recording the poten- 
tial difference value. 

The large Resistolog features and positive correlations obtained 
in this area have also been obtained on other surveys of similar depth, 
particularly along this trend. It is evident from reference to the elec- 
trical log of Humble Oil and Refining Co. #9 McKinney, on the section 
of Fig. 19, that the relatively large magnitude of these features is not 
explained by resistivity contrast as measured by electrical well logs. 
Experimental verification as to the subsurface origin of these rela- 
tively sharp features by rotating the electrode system about the point 
halfway between electrodes 1 and 2 has been shown by Fig. 8. 
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The obtaining of positively correlating Resistologs such as those 
of Fig. 12 also constitutes rather definite proof that features are of 
subsurface origin, this being particularily true if resulting correlations 
also indicate a variation in depth of subsurface beds conforming to 
that indicated by well data. Shallow inhomogeneities and errors of 
measurement, the only other likely source of such features, are random 
in nature, and would therefore practically never result in definite 
similarity, while both similarity and correlations conforming to varia- 
tion in depth of subsurface beds would be still less likely to occur acci- 
dentally. 


POSSIBILITY OF DIRECT DETECTION OF OIL AND 
GAS SATURATION—DISCUSSION 


Incremental resistivity observations which evidently reflect varia- 
tions of electrical properties of subsurface beds at interesting depths 
suggest the possibility of direct detection of oil and gas saturation. 
Due to the effects of shallow inhomogeneities it has heretofore been 
impossible to determine accurately by direct measurement whether 
oil or gas saturation at depths of practical importance could be de- 
tected by surface measurements. Williston and Nichols,’ Lee and 
Swartz’ and Karcher and McDermott® have described possibly suc- 
cessful efforts at the direct detection of oil and gas saturation by stand- 
ard four electrode DC surface resistivity measurements. However, no 
means of correcting data for shallow inhomogeneities was employed in 
any of these cases. It follows that successful detection was possible 
only in those rare cases in which the anomaly resulting from oil or gas 
saturation was large in proportion to the influence of shallow in- 
homogeneities. 

As has been frequently recognized, shallow inhomogeneities such as 
variation in moisture content of shallow beds, rock outcrops, gravel 
deposits, etc., introduce anomalies into the data of all basic electrical 
prospecting procedures. Such shallow inhomogeneities produce a 
wide variety of anomalies of various magnitudes and areal extents, 
which sometime appear to coincide with oil saturated areas either 
by chance or because of some influence of associated structural con- 


6 Williston and Nichols, U. S. Patents $1,842,361, $1,841,376, $1,842,362, #1,841,- 
975, #1,841,976, and personal communication. 

7U. S. Bureau of Mines Technical Publications $488 and #521 and U. S. Patent 
#1,951,760 to F. W. Lee. 

8 Karcher and McDermott. Bulletin of the A.A.P.G. Vol. 19, p. 64 (January 1935). 
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ditions on shallow beds. It is probable that at least some of these sup- 
posedly successful cases of direct detection were a result of coincidence 
of shallow inhomogeneities with ‘oil ar gas saturated areas. It is like- 
wise probable that failures in wildcat territory reported by Williston 
and Nichols® and Karcher and McDermott?’ were also a result of 
shallow inhomogeneities rather than the postulated variation of elec- 
trical properties of deep-lying beds. The latter conclusion is supported 
by electrical well logs. Although variations in thickness and resistivity 
of subsurface formations do occur, beds lying at the depths from 
which oil and gas are usually produced are ordinarily shown to be 
remarkably uniform laterally by electric logs. It is probable that 
electrical properties of such formations are at least no more erratic 
than are the physical characteristics upom which seismic and other 
prospecting procedures are based. 

It may be proved both experimentally and theoretically that re- 
placing all or a portion of the water which normally saturates a given 
specimen of sand with oil or gas always results in an increase in resis- 
tivity." It is also established by electrical well logs, however, that for- 
mations which are not oil or gas saturated sometimes have a lateral 
variation of resistivity of the same order as that resulting from oil and 
gas saturation. This is no doubt true of any other physical character- 
istic which might be employed for the direct location of oil or gas 
saturation. The problem of direct location of oil or gas is therefore re- 
solved into two steps: the measurement of the anomaly; and the dis- 
tinction between anomalies resulting from oil or gas saturation and 
those due to other causes. However, field data which have previously 
been described, such as those of Fig. 3, in which the relatively large 
influence of shallow inhomogeneities is demonstrated, very strongly 
suggest that the failure to successfully apply electrical methods to 
direct location of oil has thus far been entirely because of failure to 
separately meacure the anomaly. This results because of the predom- 
inating influence of shallow inhomogeneities rather than because of 
erratic variation in the resistivity of deeper subsurface beds. 

Electrical well logs provide some information as to the accuracy 
with which resistivity may be employed for determination of oil or 
gas saturation. Resistivity as measured by electric logs is not, however, 
the true in-place resistivity, since invasion of water from the drilling 


Loc. cit. 
10 Loc. cit. 
11 Martin, Murray, and Gillingham, Geopuysics, Vol. III, No. 3 (July 1938). 
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mud into permeable formations and the ‘“‘shortcircuiting”’ effect of 
the drilling mud tend to decrease the apparent resistivity of oil sands 
more than that of other formations. This results in a reduction of 
resistivity contrast of oil sands as compared to adjacent beds. The ac- 
curacy of interpretation of electrical well logs is also a function of the 
number of logs available in the immediate vicinity for comparison. It 
is believed that if electric logs were available on a sufficient nunrber 
of wells distributed over an oil productive area, including wells just 
off the saturated area in all directions, that a very high degree of ac- 
curacy of interpretation as to oil or gas saturation could be obtained 
solely on the basis of electric log data. With this quantity of data 
available it can be determined that: (1) an area of higher than normal 
resistivity exists in a particular horizon; (2) that this area has a down 
dip edge following a given structural contour; (3) and that geological 
conditions favorable for oil accumulation exist. With the possible 
exception of very shaly sands and of areas in which oil accumulation 
is controlled entirely by variation of porosity these data should be 
adequate to approximate positive determination of oil or gas satura- 
tion. The quantity of data necessary for such a scheme of inter- 
pretation is practically never available in the interpretation of electric 
logs. However, in the case of a surface geophysical method data at a 
large number of stations over the area under investigation can be 
employed for interpretation. It is thus evident that the accuracy of 
interpretation of such electrical prospecting data is primarily deter- 
mined by the accuracy with which the variation of resistivity with 
depth can be measured. If an accuracy of measurement of variation 
of depth to a given bed and the lateral variation of resistivity within 
that bed approximating that of electric logs is attained a high degree 
of accuracy of interpretation in terms of oil and gas saturation and 
structure should be possible. So high a degree of accuracy of measure- 
ment is hardly to be expected in practice but it seems evident that no 
limitation of importance exists other than the difficulty of attaining a 
measurement which accurately indicates the variation of resistivity 
of subsurface beds with depth. 


POSSIBILITY OF DIRECT DETECTION OF OIL AND 
GAS SATURATION—FIELD EVIDENCE 


Seven Sisters Field, Duval County, Texas 


Fig. 13 shows the results of an experimental survey made in an 
effort to make direct measurement of the anomaly resulting from oil 
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DETERMINATION OF RESISTOLOG ANOMALY RESULTING FROM OIL 
SATURATION AND RELATION BETWEEN DEPTH AND ELECTRODE 
SEPARATION 
SEVEN SISTERS FIELD, DUVAL COUNTY, TEXAS. 
ELECTRIC LOG 
Reliance Oil @ Rayalty Co. Atlantic Oil Prod. Co. RESISTIVITY RESISTOLOG N24-) 
PL Telford N23 J F Welder N9A-6 Welder N°A-6 = | RESISTOLOG N25 
Sec 386 Sec. 385 Telford N23 : 
2000 
9 
2 
} 
2500 2500 
\ \ 
waren 2500 / 
8 
Humble Westheimer 6 Daube LEGEND 
R. Serna 384 224 Z. Campos 
3651386, @ OIL WELL IN GOV'T. 
Atlantic Reliance ORY 
J. F. Welder PL. Telford > “HOLE. 
0 1000° 


Fic. 13. Illustrating the change in electrical well log (next to last panel) and Resisto- 
log (last panel) in passing from a well on the water side of the oil-water contact (Reli- 
ance Telford #3) to a well on the oil side (Atlantic Welder #A-6) in the Seven Sisters 
Field of Duval County, Texas. 


saturation. This traverse was across the East edge of the Seven Sisters 
field of Duval County, Texas. In this case the water-oil contact in both 
the Government Wells and Chernosky sands occurs between the 
stations 4 and 5 for which Resistologs are shown. Both the Chernosky 
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sand at a depth of 2200 feet and the Government Wells sand at a depth 
of 2500 feet are oil saturated in Atlantic #A-6 Welder and water 
saturated in Reliance #3 Telford. The resistivity increment curves 
obtained in the vicinity of these wells are shown plotted against actual 
depth by use of the depth-electrode separation curve of Fig. 18. The 
difference between these two Resistologs is shown on the extreme 
right of Fig. 13 alongside the difference in electric log resistivity of the 
above wells. In this case it was not only necessary that two Resistolog 
maxima be present at the station over the oil saturated area which 
were absent off this area, but also for these features to occur at a dif- 
ference in electrode separation corresponding to a difference in depth 
of 300 feet. 

The apparent relation between Resistolog features and the occur- 
rence of oil saturation was verified by other traverses in the Seven 
Sisters and the nearby Loma Alto fields. 

To determine whether results were influenced by well casings and 
pipe lines the north edge of oil saturation in the Government Wells 
sand in the Seven Sisters field was located in advance of development. 
Subsequent drilling established the edge of saturation to be between 
the two limiting stations at this point, these stations being 600 feet 
apart. 

As indicated by the curve of Fig. 18 an oil sand at 2200 feet pro- 
duced a Resistolog maximum at an electrode separation of 2275 feet, 
while an oil sand at a depth of 2500 feet produced a Resistolog Maxi- 
mum at an electrode separation of 2575 feet. Electrode separation is, 
of course, the distance a of Fig. 4. 


SamFordyce Field, Starr County, Texas 


The application of the Resistolog procedure to the determination 
of the extent of an already discovered oil field is shown by Figs. 5, 6, 7, 
14, 15, and 16. The area covered by this survey was the northwest end 
of the SamFordyce field of Starr and Hildago Counties, Texas, this 
survey having been confined to the Starr County portion as shown 
by the map of Fig. 14. The primary purpose was to determine the 
extent of the oil saturated area, in advance of development. In the 
area surveyed all oil and gas production is from the SamFordyce sand 
of the Frio formation, at a depth of 2800 feet. Deeper sands later 
proved productive on other portions of the SamFordyce structure but 
were not productive in the area under consideration. 
At the time of the initiation of this survey Phillips and Barnsdall 
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#1 Ytturia was the northwesternmost oil or gas well in the Sam- 
Fordyce field. Stations 1, 2, and 3 were taken over the known oil - 
saturated area for comparison with similar data off this area, such as 
at stations 5 and 32. The various resistivity increment curves, or 
Resistologs, obtained in the area are shown in the upper right hand 
corner of Fig. 14. The scheme of plotting is such that the inflection 
point associated with a resistivity maximum produces a resistivity 
increment curve peak having an apex which points to the right. The 
upper group of curves are those on which electrode separation was 
varied by moving the potential electrodes by 150 foot increments. 
Potential electrodes were moved by 75 foot increments in the case of 
the lower group. The influence of oil or gas saturation is obvious from 
comparison of the various Resistologs over the oil saturated area with 
those off this area. It will be noted that these Resistologs may be di- 
vided into two well defined groups of definitely similar curves, those 
on the oil saturated area, and those off this area. These Resistologs are 
shown exactly as computed, except that each has been rotated to ap- 
proximately the parallel position. Each Resistolog is the average of 
two or more independent Resistologs which were similar as to features. 
No adjustment of values or elimination of ‘‘bad”’ data has been em- 
ployed. 

A more quantitative determination of the net anomaly apparently 
resulting from oil saturation is shown on the right-hand side of Fig. 
16. In this case a Resistolog over the oil saturated area (Station 17 on 
Highway 12) is subtracted from a Resistolog off this area (Station 15, 
600 feet to the northwest.) It will be observed from these data that 
a Resistolog maximum is produced at an electrode separation a of 
2950 feet by an oil sand at a depth of 2800 feet. 

Resistologs covering sufficient range of electrode spacing show cor- 
relatable features which persist regardless of the presence or absence 
of the features attributed to oil or gas saturation. This fact is also 
illustrated by the Resistologs of Fig. 17. Except for short range curves 
for filling in details between control stations, the procedure in the 
Sam Fordyce field was first to correlate Resistologs with those at 
known points and then to note whether a maximum occurred at the 
point corresponding to the Sam Fordyce sand as determined by trav- 
erses across known edges of oil saturation in this formation. With 
this scheme of interpretation the location of a subsurface boundary in 
a given horizon is not necessarily based on the disappearance from 
the Resistolog of the maximum (or minimum) resulting from this hori- 
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zon. It is only necessary that the position of the Resistolog of this 
feature change relative to that of maxima or minima on which the cor- 
relation is based. It is believed that the very sharp determination of 
the edges of oil saturation in the Sam Fordyce area with this procedure 
is a result of such a shift in position on the Resistolog of the corre- 
sponding feature rather than to any abrupt decrease in magnitude of 
this feature. The depth to the SamFordyce sand was approximately 
the same throughout the area covered by this survey. This fact, 
along with the relatively large increment by which potential electrodes 
were moved, resulted in all Resistologs correlating at the same elec- 
trode separation. 

Because of the apparent fidelity with which Resistologs delineated 
the known oil saturated area it was decided to take so great an amount 
of data of the maximum accuracy obtainable as to definitely preclude 
accidental agreement with subsurface conditions. With this in view 
over three months were devoted to taking date in the relatively small 
area covered by the map of Fig. 14. These data were also subjected to 
the various auxiliary tests of validity. 

It will be noted by reference to Fig. 14 that interpretation as to 
oil or gas saturation on the above basis coincides with the known or 
subsequently determined oil productive area in considerable detail 
both on traverses across the oil-water contact and across edges of the 
saturated area in which the oil sand was replaced by shale. 

The isolated oil productive area contained within the shaded 
area of Fig. 14 was discovered as a direct result of this survey. H. J. 
Porter et al. (Porter and Grimes) #1-C Garcia, the discovery well for 
this area, was an offset to two dry holes, Porter #1-B and #2-B Garcia, 
and was beyond the limits of the field as indicated by the then existing 
subsurface information. 

The barren area in the vicinity of Porter #1-B and #2-B Garcia 
was indicated in advance of drilling by the Resistolog at station 7, and 
later detailed by additional stations. The barren area in the vicinity of 
Sun #3 and #4 Chapotal and #1 Garcia was also indicated in advance 
of drilling by stations 35, 18, and 44. 

Sections A and B of Figs. 15 and 16 respectively are taken along 
the lines indicated on the map of Fig. 14. Sand thickness and satura- 
tion shown are as determined by drilling. In this case Resistologs are 
plotted against actual depth by the relation between depth of pros- 
pecting and electrode separation shown by Fig. 18. The latter curve 
was determined by traverses across known subsurface boundaries in 
areas of generally similar subsurface conditions. 
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Fic. 18. Empirical relation between the electrode separation and the depth of prospecting, 
compiled from traverses across the oil-water contacts and other known subsurface discontinui- 
ties along the SamFordyce, Pettus, and Jackson trends of South Texas. 


As previously discussed a number of auxiliary tests are applicable 
to electrical prospecting data of this type by which the origin of any 
maximum or mimimum observed may be rather conclusively estab- 
lished. One such test is that of similarity of results at nearby stations. 
Subsurface beds are generally quite uniform laterally. Resistologs at 
stations 100 and 200 feet apart should therefore ordinarily be practi- 
cally identical both as to magnitude of features and the electrode sepa- 
ration at which they occur. The duplication throughout an area of 
Resistologs having numerous features such as those of Figs. 10, 11, 
13, 14 and 17 is practically positive proof as to the subsurface origin 
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of these features. The two most likely sources of error, shallow in- 
homogeneities and errors of measurement, are random in nature, and 
would practically never result in Resistologs of definite similarity. 

Because of the length of the electrode spread employed in deep 
prospecting it is necessary to know the location of the subsurface area 
prospected with reference to the electrode system for accurate lateral 
location of subsurface boundaries such as the edge of saturation in a 
given sand or the trace of a fault at a desired level. The experimental 
location of the subsurface area prospected with reference to the elec- 
trode system is based on reversed traverses across a subsurface bound- 
ary. If the point in the electrode system which overlies the subsurface 
area prospected is chosen as the location of the station, the location 
of the boundary must be determined to be at the same point, regard- 
less of the orientation of the electrode system. Such reversed traverses 
at stations 15 and 16 of Fig. 14 (on Highway No. 12) disclosed the 
point overlying the subsurface area prospected to be the point half 
way between electrodes 1 and 2. 

Fig. 10 showed the results of such a reversed traverse in the Bran- 
yon field of Caldwell County, Texas, the subsurface boundary being 
a fault trace in this case. 


“2100 Foot Sand” Area of the Oakville Field, Live Oak County, Texas 


Fig. 17 shows a detailed structural and direct location survey of 
the portion of the Oakville field which is productive in the ‘‘2100 foot” 
sand. The position of this sand in the geologic column is indicated on 
the electric log of Simmons Oil Co. #7 McGriff on the cross section of 
Fig. 19. The upper 7 feet of the 2100 foot sand is oil saturated in this 
well as compared with about 17 feet of oil saturation in Simmons Oil 
Co. #3-B McGriff. The purpose of this survey was the determination 
of the extent of the oil saturated area in this sand and the associated 
structural conditions. This survey was also used for determining the 
value of employing the relation between the edge of the indicated 
saturated area and structural contours, as determined from Resistolog 
correlations, as a means of distinguishing oil saturation from other 
subsurface lateral variations in resistivity. 

Resistologs at stations 1 (on the NW-SE road; at McGriff #1-B) 
and g (at McGriff #3-B) were used as standards of comparison for 
other Resistologs in the area. Simmons Oil Company #1-B McGriff, 
where station 1 was located, was just below the water level in the 
2100 foot sand. Simmons Oil Company #3-B McGriff was the discov- 
ery well in this horizon. 
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All Resistologs covering the 2100 foot sand zone in detail are shown 
in the correlated position in the upper right-hand corner of Fig. 14. 
The dotted line designated as 2050’ represents an electrode separa- 
tion @ of this amount. The increment by which potential electrodes 
were moved was 18.75 feet. Resistologs are plotted so that the feature 
resulting from an inflection point on the apparent resistivity curve due 
to a high resistive subsurface layer produces a peak on the Resistolog 
having an apex which points to the right. Electrode separation in- 
creases in the downward direction in this Fig., the electrode separation 
being 2050 feet or greater for all Resistologs. 

All Resistologs obtained in the area are shown as computed. It will 
be noted that these curves may be divided into two well defined 
groups, those over the oil saturated area and those off this area. These 
Resistologs do not appear to have quite so high a degree of similarity 
as those of Fig. 14. This is due, at least in part, to a difference in the 
scale of plotting and to the increment by which potential electrodes 
were moved. The anomaly resulting from oil saturation, as well as 
general features, are, however, smaller in this case than in case of 
deeper prospecting. This apparent increase in size of features with 
depth will be discussed further in connection with Fig. 19. 

It will be noted, by reference to Fig. 17, that all stations except 
Station 12 which were interpreted as oil or gas saturated are within 
the oil saturated area as determined by drill stem tests, this water 
level being — 1914 feet, and that all stations indicated as oil saturated, 
including 12, are within the oil saturated area as determined by elec- 
trical logs, this water level being—1922 feet. As is evident in Fig. 17, the 
upper map, which was determined from Resistolog data, is quite simi- 
lar to the subsurface map as determined from well data after complete 
development, the latter being the lower map of Fig. 17. 

The increment of 18.75 feet by which potential electrodes were 
moved makes possible an error in determination of depth equal to the 
corresponding variation in depth, or 17 feet, even if the relation be- 
tween depth of prospecting and electrode separation remains con- 
stant. Reference to the subsurface map of Fig. 17 reveals that only the 
Resistolog at station 7 failed to check actual subsurface data by the 
amount determined by the increment with which potential electrodes 
were moved. The error in this case was 29 feet which is 12 feet greater 
than could be accounted for by lack of detail. This error possibly re- 
sulted from refraction of lines of current flow by a hard and high re- 
sistive sandstone member which was disclosed to be present in this vi- 
cinity by Simmons #6-B McGriff. This sandstone, which occurred 
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within 50 feet of the surface in this well, was either absent or softer at 
other locations in the area. 

A special test was made to investigate whether the apparent addi- 
tional detail obtained by moving potential electrodes by an increment 
of 18.75 feet rather than the usual 37.5 foot or 75 feet increment was 
not because of a small residual influence of shallow inhomogeneities. 
For this purpose a second pair of current electrodes r’ and 2’ were 
placed 18.75 feet from r and 2 respectively, and in the direction of pro- 
file. For each position of the potential electrodes separate potential 
difference measurements were made for current flow between elec- 
trodes 1, 1’, 2, and 2’ and the distant electrode 3. This provided data 
for calculation of two resistivity increment curves or Resistologs at 
stations (the point halfway between 1 and 2 and 1’ and 2’) which were 
only 18.75 feet apart. Furthermore, a given potential electrode posi- 
tion represented a separation between current and potential electrodes 
which was 18.75 feet greater for one Resistolog than the other. Any 
features occurring at a given potential electrode position would there- 
fore be apparent upon inspection. This procedure was applied to all 
stations shown in Fig. 17 having higher number than 16. Stations 16A 
and 16B were obtained in this way, these stations being only 18.75 feet 
apart. These Resistologs obviously correlate at the same electrode 
separation and therefore any influence from shallow inhomogeneities 
was negligible in comparison to the effect of subsurface features. Simi- 
lar data were taken at all stations having higher number than 16 but in 
all other cases these two entirely independent Resistologs have been 
averaged for the final curve after data of sufficient accuracy was ob- 
tained to result in two Resistologs which had similar features. Since 
these two Resistologs were averaged at the same electrode separation 
rather than at the same potential electrode positions, any influence of 
shallow inhomogeneities in the vicinity of potential electrodes will be 
decreased in magnitude by 50% and distributed between adjacent 
points. This procedure may be applied more than two times with cor- 
respondingly greater distribution and decrease in magnitude of errors. 
The Resistolog procedure may also be successively applied to such 
data by taking a ratio between the two resistivity increment values 
obtained for a given position of the potential electrodes, etc. However, 
any errors present appear to be principally errors of observation and 
are more satisfactorily reduced to a negligible quantity by averaging. 

It will be noted that the position of the “2100 foot sand” trace of 
the fault which limits the Oakville field on the Northwest as deter- 
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mined by the Resistolog survey is very close to the position indicated 
by subsurface data. The indicated displacement of this fault also 
checks that determined by drilling to within 5 feet, the actual dis- 
placement being about 38 feet. 

This fault has a displacement of only 16 feet at a déeth of 1800 feet 
and cannot be located by shortening of section or correlation of elec- 
tric logs at shallower depths. The successful location of a fault having 
negligible displacement at depths less than 1800 feet provides addi- 
tional proof as to the depth of prospecting in this area. Neither core 
drilling to shallow markers nor a seismograph survey gave any evi- 
dence of this fault, although its presence was known and a special 
effort made to find it in both cases. 

The various field surveys which have been exhibited and discussed 
are not in all cases in chronological order, but rather in the order which 
seemed best fitted to logical exposition of the subject. The order in 
which the surveys were made is as follows: Fig. 3 (Hart County, Ky.); 
Fig. 14 (SamFordyce); Fig. 13 (Seven Sisters); Fig. 12 (Pettus); Fig. 
10 (Branyon); Fig. 11 (Dunlap); Fig. 17 (Oakville); Fig. 9 (Rhode 
Ranch). 


GENERAL CHARACTERISTICS OF RESISTOLOGS 


Fig. 19 is a cross section compiled from electrical well logs which is 
approximately 38 miles in length extending from the White Creek 
field of Live Oak County, Texas, to an area immediately down dip 
from the Pettus field of Bee County, Texas. Resistologs along this 
traverse are shown in the correlated position and approximately cor- 
related with an electric log in their vicinity. 

Resistologs are arranged in order of similarity rather than geo- 
graphical order along the traverse. A group of Resistologs taken in the 
Fitzsimmons field of Duval County, Texas and a representative elec- 
tric log from this field are shown attached to this section. The Fitz- 
simmons field is about 50 miles Southwest and slightly down dip from 
the Pettus field, oil production being principally from the Pettus sand. 

The area in which each group of Resistologs was taken was in the 
vicinity of the wells at which the associated electric logs were ob- 
tained. Reference lines are shown for each group of Resistologs from 
which the electrode separation a for any point on each Resistolog may 
be determined. 

Curve A shown above the cross section of Fig. 19 shows the aver- 
age electrode separation for each group of Resistologs. The electrode 
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separation given by this curve is the distance from the point halfway 
between electrodes 1 and 2 to the nearest potential electrode, that is 
a+300’, 

The middle curve, B, shows the average potential difference value, 
for each group of Resistologs, of the point correlated with the Pettus 
sand. The electrode separation of the potential difference value used 
in obtaining this average is the distance from electrode 2 to the nearest 
potential electrode, that is, a+600’, this being the smaller of the two 
potential difference values used in calculating a resistivity increment 
value. 

The lower curve, C, shows the average potential difference value 
of the point correlated with the Pettus sand for a current flow of 5 
amperes. The values of curve C were calculated from those of curve 
B by dividing the average potential difference value of B by the aver- 

E av. 
Tav. * 5)" 

It will be noted from Fig. 19 that if the interval prospected is cen- 
tered around the Pettus sand, in accordance with the relation between 
electrode separation and depth shown by Fig. 18, the magnitude of 
Resistolog features and the similarity of Resistologs in a given area 
definitely increase with increasing depth to this sand. This phenome- 
non has been observed consistently in areas along strike from various 
points of the cross section of Fig. 19. Similar increase in magnitude of 
features with depth of prospecting have also been noted in other parts 
of the geologic column. For example, such results were obtained in 
prospecting centered around the top of the Jackson formation in an 
area in which this formation occurred at 4800 feet, and in another area 
in which prospecting was carried approximately 1000 feet below the 
top of the Pettus sand, the depth of prospecting being approximately 
4500 feet. 

Where Resistologs have extended from the surface downward they 
appear to be roughly divisible into three portions. The first of these 
portions is that extending from the surface to about 1000 feet. In this 
range very large features obviously associated with shallow beds are 
obtained. These features exhibit all of the orthodox characteristics of 
resistivity prospecting predicted by mathematical analysis on the as- 
sumption of ohmic conductivity.. Resistolog features become increas- 
ingly sharp with decrease in depth, very broad features being obtained 
at the greater depths. Resistivity anomalies of the ohmic type do not 
appear to terminate in the above range (to 1000 feet) but simply result 
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in such broad features as to render correlation increasingly difficult. 
In areas in which high resistive limestones underlie formations having 
lower resistivity such as the section of Figs. 10 and 11, very broad 
Resistolog features are obtained with large electrode separation if 
Resistologs cover sufficient range. For example, the apparent resistiv- 
ity values in the latter area decrease from the surface down to an elec- 
trode separation of roughly 4000 feet after which they gradually begin 
to increase. The deeper features employed for correlation in this as 
well as other areas appear to be superimposed on the smooth curves 
resulting from the above type of variations in subsurface resistivity. 
Resistologs on the o—1000 foot range also show marked changes in rela- 
tion between electrode separation and depth of prospecting which are 
not obtained in deeper prospecting. Several such shallow traverses on 
which electric log control was available definitely establish that change 
in resistivity of the section such as cutting out of high or low resistive 
members results in a variation in relation between depth of prospect- 
ing and electrode separation. Notwithstanding this limitation such 
surveys may yield information of considerable value, this being par- 
ticularly true if the problem is the mapping of abrupt changes in depth, 
or character, of beds such as occurs at the surface trace of a fault. 
Shallow structural mapping is also possible in areas of sufficiently uni- 
form surface formations, or if well log control is available at several 
points. Fig. 9 showed an application of the Resistolog procedure to 
shallow prospecting. Faults, such as that indicated, may be more posi- 
tively determined by placing stations closer together in the vicinity of 
the suggested fault. Faults result in an abrupt change in indicated 
depth from one value to another whereas stratigraphic changes result 
in a progressive variation of indicated depth, several intermediate 
values being obtained between the two extreme: 

The second of the three portions into which Resistologs may be 
divided is the range between 1000 feet and about 2000 feet. In this in- 
terval Resistolog features are consistently small and are either less 
consistent than at greater depth or require a higher degree of accuracy 
of measurement that it was possible to obtain consistently. Several 
successful surveys have, however, been conducted in this interval but 
the time required is greater and the results less satisfactory than for 
deeper prospecting. One case of successful direct detection of oil or gas 
saturation was obtained in this interval, this being the traverses across 
the water-oil contact of the Ezzell field of Live Oak and McMullen 
Counties which are indicated on the curve of Fig. 18. The depth to the 
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oil sand was only 1500 feet, and the survey was made under known or 
approximately known subsurface conditions. No cases of definitely 
successful direct detection of oil or gas saturation at shallower depths 
have been obtained, although on most such attempts the detail of 
Resistologs as to increment with which potential electrodes were 
moved and the accuracy of measurement were not so advantageous as 
in the Ezzell area. Field results suggest that the difficulty in obtaining 
satisfactory Resistolog features of the type used for deeper prospecting 
at depths less than 2000 feet is a result of these features covering an 
increasingly smaller depth interval with decrease in electrode separa- 
tion as well as decreasing in magnitude. 

The third of the three portions into which Resistologs may be 
divided is that deeper than 2000 feet. For a given subsurface section 
Resistolog features appear to increase in magnitude and depth interval 
covered to at least 6000 feet, the greatest depth to which prospecting 
has been attempted. As previously stated, Resistolog features are also 
usually more persistent and definitely correlatable at depths greater 
than 2000 feet. Of course, it does not follow that large features are al- 
ways obtained at depths in excess of 2000 feet, the above conclusion 
being true only for prospecting covering a given geologic column at in- 
creasing depths. In some sections such as that shown by the cross 
sections of Figs. ro and 11 no large features are obtained on Resisto- 
logs to a depth of at least 4000 feet. Since in these cases the beds deeper 
than 2000 feet are known to be mostly high resistive limes and chalks 
the suggestion is very strong that the features on which Resistolog cor- 
relations are principally based are not contrasts in ohmic resistivity. 
The resistivity contrast in some of these areas was quite similar, as 
to the first 2000 feet of section, to that along the Oakville-Pettus 
portion of the cross section of Fig. 19. In this case the first few hun- 
dred feet of section was occupied by high resistive fresh water sands 
of the Wilcox series resulting in shallow resistivity variation quite 
similar to that due to fresh water sands in the Oakville series of 
Fig. 19. The conclusion that Resistolog features are not due to ohmic 
resistivity is further strengthened by considerable evidence that areas 
in which the most positive results are obtained both as to structural 
correlation and direct location of oil are those in which subsurface 
beds are of relatively low resistivity and show little resistivity con- 
trast. This was true in the Sam Fordyce area as well as in the Seven 
Sisters, Oakville and Pettus areas for which electric logs are shown. 

The reason for the apparent in¢rease in magnitude of Resistolog 
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features with increase in depth is not established. It is also not deter- 
mined whether this variation is fundamentally stratigraphic or due to 
physical changes associated with depth changes. At least two physical 
changes in the character of geologic sections of the type considered 
are known to occur with increase in depth—(1) the porosity of the 
shales decreases and (2) the mineralization of the water with which 
formations are saturated varies, usually increasing. 

Athy has shown” that the porosity of shales decreases with in- 
crease in depth according to the following formula, 


f = fe™ 


where f is porosity, fo is the average porosity of surface clays, b is a 
constant, and z the depth below the surface. According to Athy the 
average porosity of surface clays (not fresh deposits) is 45 to 50 per- 
cent, while the porosity at 6000 feet is roughly 5 percent. 

Carpenter and Spencer have shown experimentally that consoli- 
dated sands are substantially incompressible. Results on cores from 
the Woodbine sand in the East Texas field, as well as on cores from 
other sands in other areas, indicate a maximum reduction in porosity 
of 3.6 percent at a pressure of 8000 pounds per square inch. This 
pressure corresponds to the difference between hydrostatic pressure 
and pressure due to the weight of overlying formations at a depth of 
approximately 14,000 feet. The relative incompressibility of sands is 
also indicated by failure of core analyses of sands to indicate any con- 
sistent decrease in porosity with depth which may not be accounted 
for by change in silt content, grain sorting, cementation, or similar 
change in character. 

On the basis of the above conclusions, the porosity of the shales will 
progressively decrease while that of associated sands will remain con- 
stant as these formations become deeper. Since the magnitude of Re- 
sistolog features also increases with depth the accompanying increase 
in porosity contrast between sands and shales is suggested as a possible 
explanation for this increase in the magnitude of Resistolog features. 

No extensive measurements of porosity of shales are available 
along the section of Fig. 19. However, several measurements of total 
porosity of the Pettus sand and the shale immediately above it have 


12 Athy, L. F. Amer. Assoc. Petrol. Geol. Bull. 14, 8 (1930). 
13 Carpenter, C. B., and Spencer, G. B. Oil Weekly. Dec. 23, 1940. 
Carpenter, C. B., and Spencer, G. B. U. S. Bureau of Mines Report of Investigations 
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been made in the Oakville area. The average depth of shale cores was 
2750 feet and the average total porosity 36.4%. The average depth 
of sand cores was 2773 feet and the average total porosity 37.2%. 
The method of measuring porosity was by determining the loss in vol- 
ume when the sample was pulverized to the unconsolidated state. 
Since the porosities of shales and sands are about equal in the Oakville 
area it follows from the relation between depth and porosity of shales, 
that the porosity of shales is probably greater than that of sands of the 
portion of the section of Fig. 19 which is up-dip from Oakville. Con- 
versely, on the down-dip portion of this section, the shales should have 
lower porosity than the associated sands. This is based on the assump- 
tion that no material change in the porosity of formations occurs other 
than that resulting from the increase in depta. 
Variation in porosity within shale sections also undoubtedly occurs 
since some shales are probably more compressible than others. Varia- 
tion in sorting of grains of which the shale is composed would also vary 
porosity. Relative porosity of sands and shales and variation of 
porosity within shale sections would not, of course, be reflected by the 
S. P. or “porosity” diagram of electric well logs, since this curve indi- 
cates permeability and difference in electrolytic content of formations 
relative to that of the drilling mud, rather than porosity. 

Available data fail to establish a definite relation between Resisto- 
log features and mineralization of water contained in formations at the 
level prospected. Since, along the section of Fig. 19, the Pettus sand 
was at least approximately centrally located in the subsurface interval 
prospected mineralization of the water produced from this sand should 
be indicative of the mineralization of formation waters in this interval. 

The usual mineral constituents of water from the Pettus sand is in- 
dicated by the following analysis of water from Union Producing Co. 
#18 Ray. This well is located in the Pettus field of Bee County, 
Texas (Fig. 12). Depth to the Pettus sand was 3927 feet. 


Hypothetical Combination Paris per Million by wt. 

Silica 2. 

Iron and Alumina 2. 

Magnesium Carbonate 74.2 
Calcium Carbonate 139.2 
Sodium Carbonate 217.2 
Sodium Sulfate 8.5 
Sodium Chloride 13,273.6 


As usual for oil field waters sodium chloride is the principal mineral 
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TABLE I 
Depth to . NaCl Part sper 
Pettus Sand Field County Million by wt. 

goo’ Jacobs McMullen 8,400 
2511’ So. White Creek area Live Oak 17,100 
2750’ Oakville | Live Oak 30,075 
3900’ Dirks Bee 17,480 
3930’ Pettus Bee 13,700 
‘42007 Fitzsimmons Duval 26,650 
5082’ Southland Duval 42,000 


constituent. Table I, which gives the chloride content, calculated as 
sodium chloride, of water produced from the Pettus sand in fields on, 
or along strike from, the section of Fig. 19, should therefore be indica- 
tive of mineralization of this water. 

It will be noted that the salt content of-water from the Pettus sand 
increases with depth along the section of Fig. 19 to a point in the gen- 
eral vicinity of the Oakville field, and then decreases to the Pettus 
field. The salt content of water in the Fitzsimmons and Southland 
fields suggest a general increase in salt content with increase in depth. 
It is probable that interruption of continuity of beds by faults has re- 
sulted in greater variation in mineralization of subsurface waters along 
this section than would otherwise be the case. 

It is possible that greater than normal contrast between the salt 
content of water in the Pettus sand with respect to that of adjacent 
shales exists in the Fitzsimmons and Pettus fields. Although the salt 
content of water from the Pettus sand appears to be relatively low con- 
sidering the depth, electric log resistivity values suggest that adjacent 
shales are saturated with water which is not correspondingly low in 
salt content. Table II gives average electric log resistivity of a uniform 
shale section, about 100 feet thick, which consistently occurs immedi- 
ately above the Pettus sand. This shale includes the Textularia Dibol- 


TABLE II 

No. of Av. Depth Av. Abd. Res. Est. Est. Res. 

Area Well To Pettus Res.Ohm-  Diboll Temp. Ohm-M 

“ Sand M Shale Deg. F at120°F 
White Creek 5 1897’ 3.8 1.20 104 1.06 
Oakville 9 2769’ 2.3 1.05 115 1.01 
Pettus 5 4018' 2.7 1.09 133 1.19 
Fitzsimmons 4276’ 2.5 .89 137 1.02 
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lensis zone which is employed as a micro-paleontological marker in 
this area, for which reason it is commonly designated as the Diboll 
shale. Temperature correction is approximate in Table II because 
temperature data were not available on all wells. The indicated 
resistivity values are, however, in close agreement with those of other 
wells along strike from this section, on which temperature data were 
available. 

Comparison of the electric log resistivity of the Diboll section in 
the Oakville field with that in the Pettus area is of particular interest. 
The ratio of the resistivity of water in the Pettus sand in the Pettus 
field to that of the Pettus sand water in the Oakville field is indicated 
by salt content to be .49/.24 or 2.04. The ratio of resistivity of the Di- 
boll shale section in the Pettus area relative to that in the Oakville 
field is 1.19/1.01 or 1.179. This suggests that the salt content of water 
in the Diboll shale section is greater in the Pettus area relative to the 
salt content of Pettus sand water than at Oakville. This conclusion is 
based on the assumption that the porosity of this shale section is at 
least not greater in the Pettus area than at Oakville. 

Consideration that the porosity of shales decreases with depth in 
accord with Athy’s formula* suggests that the salt content of the water 
in the Diboll shale section increases with depth. Tables IT and III sug- 
gest that resistivity of this shale section remains approximately con- 
stant as the depth increases. If the porosity decreases with depth, an 
increase in mineralization of the water with which the shale is satu- 
rated is necessary if the resistivity is to remain constant. 

Additional, and more complete, data are being obtained along the 
above trend as rapidly as made available by drilling operations. The 
above study of subsurface data in an effort to obtain a theoretical 
explanation of results established experimentally is preliminary in 
nature. 

All relatively deep prospecting in areas in which large features were 
obtained involved the measurement of quite small potential difference 
values, usually between 100 and 200 microvolts. Since approximately 
the same current flow was employed in all relatively deep prospecting 
in low resistive areas, the current density at the level prospected was 
considerably less for large electrode separation. As previously ex- 
plained, curve B of Fig. 19 shows the average potential difference value 
for the point correlated with the top of the Pettus sand. It is possible 
that the smaller current density in the case of deep prospecting had 
some influence on the magnitude of Resistolog features. However, 
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‘ several experiments performed in an effort to establish such a re- 
lation, although not sufficiently extensive to be conclusive, failed 
to indicate any variation of the magnitude or character of features 
with change in magnitude of current. One experiment of this nature 
in the Oakville area in which current was varied from 0.5 to 7.0 
amperes with potential ele¢trodes 600 feet apart and the nearest 
potential electrode 2800 feet from electrode r resulted in a ratio be- 
tween potential difference and current which varied less than five- 
tenths of one percent through the range of current variation. A num- 
ber of Resistologs in the Ezzell area which were taken with current 
flow of both 3 amperes and 6 amperes failed to show any definite in- 
crease in the magnitude of features. Some increase in magnitude was 
noted with the smaller current flow in one or two cases but the amount 
of increase was not definitely greater than could be explained by error 
of observation. 

The group of Resistologs shown in Fig. 19 with the electric log of 
Holland and Smith #1 Bonham are typical of Resistologs in areas along 
this trend in which depth to the Pettus sand is less than 2000 feet. Al- 
though features are quite small in this area it was possible to make 
fairly reliable structural correlations. The maximum error in determi- 
nation of depth disclosed by subsequent development in the area was 
30 feet. 

The group of Resistologs shown with Grimes #1 Friedle represent 
an apparently positive case in which a high resistive Resistolog feature 
resulted because of higher than normal resistivity of the Pettus sand 
due to abnormal cementation and resulting hardness of this forma- 
tion. Resistologs having the higher electrode separation are on the low 
side of a down-thrown-to-west fault in this area and near the above 
well. The existence and location of this fault has been recently verified 
by core drilling. No attempt was made to study the down dip edge 
of the high resistive area in relation to structural contours in this 
case. 

The group of Resistologs shown for the Fitzsimmons area were ob- 
tained in connection with a survey for mapping the fault which limits 
oil accumulation on the up-dip or Northwest edge. The Resistologs 
having the higher electrode separation are those on the low side of this 
fault, the remainder being near the fault on the high side. The direc- 
tion and displacement of this fault indicated by Resistolog correlation 
checked that determined by drilling at all known or subsequently de- 
termined points. An apparent shift in the location of the subsurface 
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area prospected with reference to the electrode system was noted in - 
the vicinity of this fault, the location of this area being shifted towards 
the potential electrodes by a distance of about 1320 feet or about one 
fourth the distance between the point half way between electrodes 
rand 2 and the midpoint of potential electrodes. No reversed traverses 
were made across this fault but the same relation appeared to exist 
on all traverses in the area. Well casing, pipe lines, and lease bounda- 
ries prevented reversed traverses in this case. No attempt was made 
to determine the relation between depth of prospecting and electrode 
separation or whether it was possible to detect oil saturation directly 
in this area because of limitation of time available for this survey. 
Only one other case has been encountered in which the subsurface area 
prospected was established as being displaced from the usual position 
below the point halfway between rz and 2, this also being a traverse 
across a fault. 

The group of Resistologs shown in the Oakville and Pettus areas 
by Fig. 19 have already been discussed. 

As suggested by the results at station No. 7 of Fig. 17, variations in 
the relation between depth and electrode separation have been ob- 
served in deeper prospecting although such changes are very much less 
prevalent than in shallow prospecting. These variations appear to be 
associated with the cutting out of prominent high resistive members 
near the surface. This phenomenon has been noted on several tra- 
verses across the outcrop of the Oakville sand in which 150 feet to 250 
feet of fresh water-saturated sand having resistivity about thirty 
times that of underlying beds goes out of the section within a mile or 
less. In all such cases an increase in resistivity of the surface layer re- 
sulted in the depth indicated by Resistolog correlations being shal- 
lower than the actual. This relation seems to be explained by the 
bending of current lines away from the normal in passing from the sur- 
face layer to subsurface beds of lower resistivity. The intersection of 
the current lines flowing in subsurface beds, and therefore the appar- 
ent current electrode position, would in this case be below the surface, 
resulting in an apparent decrease in depth. All variations in relation 
between depth of prospecting and electrode separation observed in 
deep prospecting appear to have been of this type. As suggested by ac- 
companying data material errors of depth determination are relatively 
rare in deep prospecting. It should be possible to make correction for 
errors of this type by means of auxiliary data for determining the 
thickness and resistivity of the surface layer, or by placing current 
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electrodes in core drill holes which have been carried through the high 
resistive surface zone. 
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INDEX OF WELLS SHOT FOR VELOCITY“ 
B. G. SWANt 


ABSTRACT 


This index lists 811 wells in which seismic velocity surveys have been made. In so 
far as possible to obtain, information has been listed concerning the location, survey 
depth, by whom and when made, and the company or co-operative association respon- 
sible for the survey. It isintended that the index be brought up to date at approximately 
six-month intervals in future issues of GEOPHYsICs. 


INTRODUCTION 


For many years those engaged in seismograph exploration have 
been measuring seismic velocities by lowering a geophone in a drill 
hole and recording shots set off at the surface. In most instances, 
many companies are quite willing to share the data from a survey with 
others who are willing to share the expense either at the time the sur- 
vey is made or at some later date. Such cooperation has been difficult 
in the past because information regarding the existence and location of 
old surveys has not been readily available. The chief purpose of this 
index is to overcome this difficulty. It is to be hoped that it will 
eventually lead to a wider distribution of velocity information and 
consequent improvement in the general level of seismic interpretation. 
It should be emphasized, however, that the listing of a survey in the 
index does not in any way imply that the data are available for sale or 
trade. 

Without the almost unanimous response of all the seismograph 
contracting companies, co-operative well velocity survey associations, 
and oil companies which conduct their own geophysical work, the 
preparation of this index would not have been possible. This co-opera- 
tion has been greatly appreciated by the author and the editorial staff 
of GEOPHYSICS. 

Information from which the index is compiled was obtained by. 
direct solicitation from the sources mentioned above. A preliminary 
index was furnished, and each source was requested to make correc- 


* Society members who are particularly concerned with well velocity determinations 
can secure a reprint of this index by addressing a request to the Business Manager. 

+ Continental Oil Company, Ponca City, Oklahoma. Serving as Special Editor to 
assemble this index. 


540 


INDEX OF WELLS SHOT FOR VELOCITY 541 


tions and additions from information in their files. From these 
amended and revised lists, the comprehensive index was assembled. 

It is to be expected that some errors in the initial index will be 
discovered by those who use it to any extent. It is requested that all 
such errors which withstand careful checking be reported to the editor 
of Gropuysics and likewise information that is at present lacking on 
listed wells. Since it is planned to bring the index up to date at ap- 
proximately six-month intervals, it is also requested that information 
regarding new surveys made in the future be reported voluntarily by 
all those having knowledge of the survey. This is important if the in- 
dex is to be kept current. 

The Society plans to prepare enough reprints of the index to pro- 
vide on request a copy, gratis, to any man in any company who is par- 
ticularly concerned with well-velocity determinations. 

It is intended that the “Depth” heading refer. to the maximum 
depth of the survey. This definition was not made clear when informa- 
tion was first solicited from the various companies, and it is suspected 
that many may have reported instead the total depth of the well. Cor- 
rections reported for errors in this figure will be greatly appreciated. 

Little difficulty should be encountered in interpretation of the ab- 
breviations in the “Shot By” column. Under ‘Sponsored By,” how- 
ever, the meanings of the abbreviations are as follows: 

OK WSA= Oklahoma-Kansas Well Shooting Association 

CWVSG= Co-operative Well Velocity Survey Group (California) 

SSWSA = Southern Seismic Well Survey Association 
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DISCUSSION AND COMMUNICATIONS 


[In response to a request by the Editor, Mr. M. B. Dobrin, who acted as instructor 
in an ESMWT course in petroleum geophysics which was presented in Washington this 
Spring and Summer, has prepared the following account of the course for the assistance 
of those who may be interested in organizing similar courses in the geophysical centers.] 


DESCRIPTION OF ESMWT COURSE IN PETROLEUM GEO- 
PHYSICS CONDUCTED IN WASHINGTON, D.C. 


Although 700 miles away from the nearest area where geophysical field parties or- 
dinarily prospect for oil, thirty professional men and women in the vicinity of Washing- 
ton, D. C., have evinced a sufficiently strong interest in petroleum geophysics to register 
for a fifteen week course in this subject when it was offered during the past spring and 
summer by George Washington University. The course was a part of the Engineering, 
Science, and Management War Training Program sponsored by the U. S. Office of 
Education to provide specialized training for work in essential war activities where 
there is a critical shortage of technically trained personnel. The government finances all 
instruction offered under this program, but leaves its administration to established insti- 
tutions of collegiate standing. The program is in operation at several colleges and uni- 
versities located near centers of geophysical activity and similar courses could probably 
be established at such places if there is sufficient demand. A brief description of the 
course taught at Washington might be useful to anyone interested in starting a class of 
the same kind elsewhere. 

Despite an acute shortage of trained technical men for petroleum prospecting, ad- 
mittedly essential war work, no geophysics courses had been offered under the auspices 
of the ESMWT program during the first three years of its existence. In January, 
1944, the writer, believing that there was considerable potential interest around 
Washington in a Petroleum Geophysics course, approached Prof. Frank A. Hitch- 
cock, in charge of the war training program at George Washington University, of- 
fering to teach such a course if the University would undertake its sponsorship. Prof. 
Hitchcock agreed to apply to the U. S. Office of Education for its approval. This was 
granted after a preliminary canvass, which Mr. J. F. Gallie made at the Interior De- 
partment, showed that seventeen persons there would be interested in enrolling. Classes 
were scheduled in two hour sessions two evenings a week for fifteen weeks beginning 
May 1. 

Of the thirty people that enrolled, eleven were employed by the Petroleum Admin- 
istration for War, three by the Geophysical Instrument Company, three by the Naval 
Ordnance Laboratory, two by the U. S. Geological Survey, and the remainder by vari- 
ous regular and wartime government agencies. Most of the PAW representatives were 
geologists, economic analysts, or refinery technologists. Four members of the class were 
in the Army, three of these being officers on detail to the Refining Division of the PAW. 
Three were women, two of them geologists, the third an attorney who dealt with oil 
lands and leases for the Interior Department. Almost all of those in the class had college 
degrees, three holding doctorates. Only two had been engaged directly in petroleum 
prospecting by geophysical methods, but one had had twelve years’ experience on five 
continents. 

For a group so heterogeneous in background and interests, it was difficult to decide 
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on a curriculum. The vast majority of the enrollees were too well established in other 
professions to look upon the course as vocational training for employment in geophysical 
work. These were more interested in learning the fundamental principles and applica- 
tions of the various geophysical methods from the viewpoint of the oil company execu- 
tive, geologist or technical expert in auxiliary or related fields. The textbook chosen for 
the course, “Geophysical Prospecting for Oil” by L. L. Nettleton, was especially appli- 
cable for a class of this composition. 
The subject matter covered during each of the fifteen weeks was as follows: 


1. Introductory material; fundamentals of petroleum geology 
2. General survey of geophysical methods 
3. Gravitational methods—fundamental principles 
4-5. Gravity reductions and interpretation 
6. Gravity instruments and field operations 
7. Magnetic methods—fundamental principles 
8. Magnetic interpretation 
g. Magnetic instruments and field operations 
Seismic prospecting—fundamental principles 
11. Refraction methods—theory and interpretation 
12. Reflection methods—theory and interpretation 
13. Reduction of sample reflection records (problem) 
14. Seismic instruments and field operations 
15. Electrical prospecting and well logging; geochemical methods. 


Class work consisted of lectures, discussions, and the working of problems requiring 
the reduction, plotting, contouring, and interpretation of geophysical field data. The 
problems were kindly furnished by Dr. Nettleton, who had devised them for courses he 
had taught in the subject at the University of Pittsburgh. The Gulf Research and De- 
velopment Co. furnished prints of six seismograph records representing a mile of reflec- 
tion shooting. These were worked up as a project for the class. Often various members of 
the class would bring in such valuable demonstration material as actual seismic reflec- 
tion maps, electrical well logs, and even geophysical instruments. 

Prof. Hitchcock obtained authorization for the instructor to invite a limited num- 
ber of guest lecturers to address the class. Under this arrangement, the following speak- 
ers were heard, all being well known for their attainments in various geophysical fields: 


1. Malvin G. Hoffman, Chief, Reserve Division, PAW: “Isostasy and its Relation 
to Gravity” 

2. George B. Somers, Senior Physicist, Bureau of Ordnance, Navy Department: 
“Gravity Field Operations” 

3. L. L. Nettleton, Staff Geophysicist, Gulf Research and Development Co.: 
“Gravity Effects around Salt Domes” 

4. A. G. McNish, Magnetician, Department of Terrestrial Magnetism, Carnegie 
Institution of Washington: “The Earth’s Magnetic Field” 

5. L. C. Paslay, Partner, National Geophysical Co., and Physicist, Naval Ordnance 
Laboratory, “Seismic Instrumentation” 


Although several members of the class remarked on occasion about the usefulness 
of the course in connection with their work, it must be acknowledged that its value in 
helping relieve the shortage of trained technical manpower in petroleum prospecting was 
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limited by two factors. One is the geographical separation of Washington from the areas 
of geophysical activity; the other is the fact that all those who enrolled in the course 
were already engaged in essential war work. However, as pointed out earlier, the same 
machinery is available for establishing similar courses of instruction in cities where they 
might be put to more immediate use. If the course given in Washington can set an exam- 
ple leading to further activity in the same direction it will have served its major pur- 


pose. 
Mitton B. Dosrin 


Naval Ordnance Laboratory Instructor 
Washington, D. C. 
August 27, 1944 


PATENTS* 


ELECTRICAL PROSPECTING 
U.S. No. 2,343,140. H. M. Evjen. Iss. 2/29/44. App. 1/8/42. Assign. Nordel Corp. 
Apparatus for Making Geophysical Explorations. An apparatus using commutated 
d-c for making four electrode earth resistivity measurements and made direct reading 
by passing a part of the ground current through a separate commutator and transformer 
whose secondary is synchronously rectified to produce d-c for feeding the potentiometer 
against which the earth potential is balanced. 


U.S. No. 2,344,672. H. Blasier. Iss. 3/21/44. App. 3/23/40. 


Method of Measuring Earth Potentials. A method of measuring potential differences 
between buried conductors and the ground or between ground points by using an elec- 
trode insulated from the ground and measuring the charging current of the condenser so 


formed. 


U.S. No. 2,345,608. F. W. Lee. Iss. 4/4/44. App. 5/7/40. 


Geophysical Prospecting. A method of making four electrode resistivity or im- 
pedisivity measurements to determine dip and strike by comparing potentials trans- 
verse to and in line with the current electrodes, and rotating the entire configuration 
until an azimuth of symmetry is found. 


GEOCHEMICAL PROSPECTING 
U.S. No. 2,345,219. R. F. Sanderson. Iss. 3/28/44. App. 12/30/39. Assign. Stanolind Oil ~ 
and Gas Co. 


Geochemical Prospecting. A method of collecting soil gas samples in which a con- 
tainer with hydrocarbon sorbent material is suspended in a shallow hole, allowed to re- 
main a definite length of time, then closed and removed, the container having material 
to remove acid, alkali and water from the soil gas before adsorption of hydrobarcons. 


U.S. No. 2,346,735. P. F. Daugherty and P. D. Barton. Iss. 4/18/44. App. 8/22/41. 
Assign. Sun Oil Company. 


Method of Detecting Petroleum Deposits. A method of soil analysis prospecting in 
* Abstracts by O. F. Ritzmann, Gulf Research & Development Company. . 
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which soil samples or drill cores are heated to determine the amounts of water and of 
low boiling hydrocarbons and the ratio used as an index to locate petroleum. 


U.S. No. 2,348,103. A. O. Beckman. Iss. 5/2/44. App. 1/31/40. Assign. American Geo- 
chemical Corp. 


Method of Soil Analysis for Location of Oil Deposits. A method of soil analysis pros- 
pecting in which samples from shallow holes are analyzed for sulphide content. 


U: S. No. 2,349,250. R. L. Doan. Iss. 5/23/44. App. 10/9/39. Assign. Phillips Petroleum 
Co. 
Gas Detection. A method of analyzing soil gas for various hydrocarbons by placing it 
in a closed container having a platinum filament whose temperature change is observed 
in a bridge circuit and tested at successively higher filament temperatures. 


U.S. No. 2,349,472. M.S. Taggart, Jr. Iss. 5/23/44. App. 6/7/39. Assign. Standard Oil 
Development Co. 
Oil Prospecting Method. A method of soil analysis in which samples are placed in a 
sealed container with hydrocarbon gas and the rate of consumption of the gas observed 
to determine the amount of hydrocarbon consuming bacteria present. 


GRAVIMETRIC PROSPECTING 


U.S. No. 2,346,593. A. R. Lindblad and J. D. Malmqvist. Iss. 4/11/44. App. 1/5/40 and 
1/2/41. 
Damping Device in Gravimeters. A gravimeter having two suspended condenser 
plated and filled with a liquid to provide damping which is symmetrical about the null 
position. 


U.S. No. 2,349,404. D. W. Blair and F. G. Boucher. Iss. 5/23/44. App. 7/22/41. Assign. 
Standard Oil Development Co. 


Clamping Means for Force-Responsive Elements. A non-sticking suspended system 
clamp which has auxiliary fine wire springs to free the system from the main clamp, the 
wires being removed from contact by further motion of the clamp screw. 


SEISMOGRAPH PROSPECTING 


U.S. No. 2,345,288. W. E. Pugh. Iss. 3/28/44. App. 9/22/42. Assign. Seismograph Serv- 
ice Corp. 
Method of Seismic Prospecting. A method of obtaining continuous subsurface seis- 
mograph control by shooting short and long shots into the same detector set-up from 
both ends. 


U.S. No. 2,346,369. J. D. Eisler. Iss. 4/11/44. App. 2/28/42. Assign. Stanolind Oil and 

Gas Co. 

Seismic Surveying. A distortionless seismograph amplifier with a band pass filter 
having a linear phase shift curve and whose cut-off frequencies are designed so that an 
extrapolation of the phase shift curve to zero frequency will give an apparent phase shift 
which is an integral multiple of 7 radians. 
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U.S. No. 2,348,225. O. S. Petty. Iss. 5/9/44. App. 2/13/40. 

Magnetic Seismometer. A variable air gap reluctance type seismometer whose sus- 
pended system consists of an annular magnet and internal coil, the magnetic return cir- 
cuit and air gap being inside the magnet. 


U.S. No. 2,348,245. P. H. Dijksterhuis. Iss. 5/9/44. App. 4/9/41. Vested in Alien Prop- 
erty Custodian. 

Device for Converting Variations of a Mechanical Quantity into Variations of an 
Electric Voltage. A modulated carrier system of measuring displacement in which the 
demodulated output controls the grid bias of an amplifier tube, so as to compensate for 
non-linearity in the transducer itself. 


U.S. No. 2,348,401. L. Manzanera. Iss. 5/9/44. App. 7/28/42. Assign. Olive S. Petty. 

Apparatus for Recording Wave Form Signals. A toothed wheel for putting timing 
lines on seismograph records and so designed to mark small intervals only on the edges 
of the record and large intervals all the way across. 


U.S. No. 2,348,409. J. O. Parr, Jr. Iss. 5/9/44. App. 6/14/41. Assign. Olive S. Petty. 

Seismic Surveying. A method of seismograph recording in which mixing is auto- 
matically varied with time by means of a motor driven cam, first arrivals being re- 
corded unmixed and the degree and extent of mixing subsequently increased. 


U.S. No. 2,348,411. O. S. Petty. Iss. 5/9/44. App. 2/8/41. 

Method for Correlating Seismographic Curves. A method of redrawing seismograph 
records taken at different or variable tape speeds by tracing the record onto a trans- 
parent template also having timing marks and sliding it into proper position as the 
tracing progresses. 


U.S. No. 2,349,186. E. Merten. Iss. 5/16/44. App. 8/31/42. Assign. Shell Development 
Co. 
Sensitivity Control for Seismic Recording. A seismograph expander in which resistors 
shunting the detectors are initially immersed in mercury, the removal of the mercury 
being initiated by the first arrival and its rate controlled by the recorded signals. 


U.S. No. 2,350,803. R. D. Newcomb. Iss. 6/6/44. App. 6/18/41. 

Automatic Volume Control. An avc circuit for audio amplifiers in which a portion of 
the output is rectified and applied to the grids of two tubes whose plate circuits control 
the impedance of a center tapped shunting transformer. 


U.S. No. 2,351,456. N. H. Ricker. Iss. 6/ 13/44. App. 8/29/41. Assign. Standard Oil 
Development Co. 
Seismic Prospecting. A method of seismic recording in which the amplifier has a fre- 
quency characteristic which is the reciprocal of the ground attenuation-frequency 
characteristic and also a linear phase shift. 


U.S. No. 2,351,524. R. L. Lay, B. H. Treybig, Jr. and L. M. Hubby. Iss. 6/13/44. App. 
5/23/42. Assign. Texaco Development Corp. 


Method of Seismic Prospecting. A method of avoiding secondary waves in water 
shooting by suspending the charge from a float at such depth to allow the water to blow 
out above the charge. 
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WELL LOGGING AND SURVEYING 


U.S. No. 2,342,827. A. L. Ackers. Iss. 2/29/44. App. 5/28/40. Assign. Stanolind Oil & 

Gas Co. 

Apparatus for Well Logging. A device for locating fluid seneteinil in a well consisting 
of an electric heating coil with an adjacent thermostat which is lowered through the 
fluid and recording the operation of the thermostat so as to log the power required to 
maintain constant temperature. 


U.S. No. 2,344,598. W. L. Church. Iss. 3/21/44. App. 1/6/42. 


Wall Scraper and Well Logging Tool. A fluid pressure operated device for the end of 
the drill stem which scrapes the sides of the hole and collects a sample in a removable 
barrel and also equipped with a-c electrical logging electrodes. 


U.S. No. 2,345,119. G. D. C. Hare. Iss. 3/28/44. App. 6/26/40. Assign. Texaco Devel- 
opment Corp. 
Subsurface Prospecting. A method of well logging through the casing by lowering a 
source of gamma rays and a source of neutrons and simultaneously detecting with ap- 
propriate ionization chambers the scattered gamma radiation and scattered neutrons. 


U.S. No. 2,345,770. S. Park. Iss. 4/4/44. App. 10/5/40. Assign. Eastman Oil Well Sur- 
vey Corp. 

Means for Surveying Well Bores. A well surveying instrument which is stalled 
through the drill stem just before the drill pipe is pulled and which records photo- 
graphically at intervals during removal of the pipe, a record of the mane orientation 
of the drill pipe being hapt at the surface. 


U.S. No. 2,345,829. W. J. Opocensky. Iss. 4/4/44. App. 1/20/40. Assign. Eastman Oil 
Well Survey Corp. 
Means for Marking Record Elements. A device for marking accurately concentric 
reading lines on well surveying instrument disks and having a marking block which cuts 
lines into the disk when the block is rotated. 


U.S. No. 2,345,864. J. C. Arnold. Iss. 4/4/44. App. 6/6/42. Assign. Lane-Wells Co. 

Conductor Core for Wire Ropes and Method of Manufacture. A welllogging cable con- 
sisting of a rubber insulated conductor reinforced with spiral cords which are covered 
with synthetic rubber and wound in opposite directions. 


U.S. No. 2,346,203. W. M. Zaikowsky. Iss. 4/11/44. App. 12/7/40. Assign. Consoli- 
dated Engineering Corp. 
Well Logging Method. A method of mud logging in which circulation is saenaed fora 
short time to permit formation fluids to enter the mud and when circulation is resumed, 
analyzing samples taken from the mud corresponding to various well depths. 


U.S. No. 2,346,481. A. D. Garrison. Iss. 4/11/44. App. 11/27/40. Assign. Texaco De- 
velopment Corp. 

Determination of Underground Strata. A method of fluorescence logging of uncased 
hole by transmitting ultra-violet light from a lamp to the wall through a quartz rod, 
and with another quartz rod transmitting the fluorescence to a filter and to a moving 
photographic film. 
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U.S. No. 2,346,789. E. E. Roper. Iss. 4/18/44. App. 6/15/42. Assign. Stanolind Oil and 
Gas Co. 
Well Logging. A method of eliminating the effect of potassium on radioactivity logs 


by taking samples of the formations, spectrographically determining their potassium 
content and its effect on the instrument and subtracting this from the log. 


U.S. No. 2,347,794. R. G. Piety. Iss. 5/2/44. App. 3/15/41. Assign. Phillips Petroleum 
Co. 
Well Surveying Device. A method of electrical well logging using a surface electrode 
and a long cylindrical electrode having in the middle an insulated section whose current 
or potential may be measured separately as the assembly is moved in the hole. 


U.S. No. 2,349,225. S. A. Scherbatskoy, G. Swift, R. E. Fearon and J. Neufeld. Iss. 

5/16/44. App. 9/15/41. Assign. Well Surveys, Inc. 

Well Logging Instrument. A gamma ray well logging apparatus having an ionization 
chamber and two resistors in series, one of which is fed by current supplied and recorded 
at the surface, the condition of balance of the two resistor potentials being automati- 
cally signaled to the surface and automatically adjusted. 


U.S. No. 2,349,366. C. A. Moon. Iss. 5/23/44. App. 3/12/42. Assign. Socony-Vacuum 
Oil Co., Inc. 
Method of Geophysical Prospecting. A method of well logging in which a radioactiv- 
ity logging instruent is dropped into the drill stem when it is raised off bottom to add a 
new length of pipe, the log being taken while the pipe is lowered to bottom, and the in- 
strument then fished out on a wire line. 


U.S. No. 2,349,712. R. E. Fearon. Iss. 5/23/44. App. 11/10/38 and 1/28/43. Assign. 
Well Surveys, Inc. me 


Well Logging Method and Device. A well logging device having a neutron emitting 
source and a properly shielded ionization chamber to detect scattered gamma radiation 
produced in the formations by the neutrons. 


U.S. No. 2,349,753. B. Pontecorvo. Iss. 5/23/44. App. 2/5/42. Assign. Well Surveys, 
Inc. 


Method and A pparatus for Geophysical Exploration. A method of gamma ray well 
logging in which the hardness of radiation is determined by using a divided Geiger 
counter and measuring the ionization in each side and also the coincidences. 


U.S. No. 2,350,371. G. A. Smith. Iss. 6/6/44. App. 7/23/40. Assign. Sperry-Sun Well 
Surveying Co. 


Borehole Logging Apparatus. An acoustic logging device having a sound transmitter 
and receiver with electrical feedback to cause oscillation whose frequency depends on 
the rock velocity. 


U.S. No. 2,350,832. F. F. Segesman. Iss. 6/6/44. App. 2/21/41. Assign. Schlumberger 
Well Surveying Corp. 


Electrical Depth Marker. An electrically detectable marker for wells consisting of a 
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ring of an electrically dissimilar metal or an insulating ring on the inside surface of the 
casing, or in open hole a gun perforator bullet made of two metals constituting an elec- 
trolytic couple. 


U.S. No. 2,351,028. R. E. Fearon. Iss. 6/13/44. App. 9/17/41. Assign. Well Surveys, 
Inc. 


Well Surveying Method and Apparatus. A method of gamma ray well logging in 
which the ion pulse produced in the ionization chamber is proportional to the gamma 
ray energy and having a recording circuit which will record each separate pulse so that 
a gamma ray frequency spectrum may be obtained. 


MISCELLANEOUS 


U.S. No. 2,342,685. A. S. Norcross. Iss. 2/29/44. App. 5/8/41. 

Measuring System. A follow-up device for continuously measuring and indicating a 
condition manifested as a high electrical resistance and using a thermionic tube to main- 
tain constant current through a resistance in series with the one to be measured.. 


U.S. No. 2,343,063. E. L. Kent. Iss. 2/29/44. App. 10/31/42. Assign. C. G. Conn, Ltd. 


Method of and Apparatus for Measuring Mechanical Vibrations. A mechanical vi- 
brometer having a mass on an arm which is supported by damped springs whose angle 
with the arm may be changed to vary the natural frequency of the system. 


U.S. No. 2,344,296. F. W. Frink. Iss. 3/14/44. App. 1/25/38 and 11/29/40. 


Means and Method for Cathode Ray Oscilloscope Observation and Recordation. A 
method of observing several cathode ray tube screens simultaneously by looking at one 
with one eye and at a second with the other eye or by optically superimposing images of 
the screens through the use of partially reflecting mirrors, or by observing the screens in 
rapid succession by using a rotating mirror. 


U.S. No. 2,344,344. C. Dorsman, M. F. Reynst and R. Veldhuizen. Iss. 3/14/44. App. 
2/26/41. Vested in Alien Property Custodian. 


Device for Oscillographing Electrical or Mechanical Processes. A device for balancing 
a remotely located a-c bridge circuit whose unbalance is observed with a cathode ray 
oscilloscope, by introducing into the observed bridge voltage a component derived from 
the a-c exciter and adjusting its phase and amplitude. 


U.S. No. 2,344,642. A. C. Ruge. Iss. 3/21/44. App. 9/16/39 and 8/5/42. 

Temperature Compensated Strain Gauge. A filament type electric strain gauge having 
a strain-free temperature compensating element suspended as a cantilever slightly above 
the strain element and the two elements connected in adjacent arms of a Wheatstone 
bridge circuit. 


U.S. No. 2,344,647. E. E. Simmons, Jr. Iss. 3/21/44. App. 2/23/40and 1/1/42. 


Method and Apparatus for Making Strain Gauges. A production jig for winding 
filament type electric strain gauges which insures tautness of the wires and uniformity in 
characteristics. 
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U.S. No. 2,344,648. E. E. Simmons, Jr. Iss. 3/21/44. App. 6/20/42. 

Moisture Proof Strain Gauge. A filament type electric strain gauge which is bonded 
to the inside of a sealed tube also containing a strain free dummy for temperature com- 
pensation, the tube being arranged to engage the test body. 


U.S. No. 2,345,654. C. H. Bowman. Iss. 4/4/44. App. 6/6/38. Assign. Trojan Powder 
Co. 
Explosive Cartridge. A waterproof explosive wrapper having longitudinal scores on 
the cylindrical surface so that when end pressure is applied in tamping the explosive ex- 
pands in diameter to completely fill the hole. 


U.S. No. 2,345,887. J. E. Rothrock. Iss. 4/4/44. App. 4/5/40. Assign. Hercules Powder 
Co. 
Method of Sealing Containers. A method of sealing plastic explosive containers by 
having an internal recess filled with mastic material into which a pre-cooled internal 
plug is thermally expanded. 


U.S. No. 2,345,935. G. L. Hassler. Iss. 4/4/44. App. 4/6/42. Assign. Shell Development 
Co. 
Method and Apparatus for Permeability Measurements. A method of measuring per- 
meability of rock samples under polyphase fluid flow by forcing fluid around a closed 
system and measuring the pressure drop across the sample in a separate closed system. 


U.S. No. 2,346,043. K. J. Mysels. Iss. 4/4/44. App. 3/2,//42. Assign. Shell Development 
Co. 


Pipe-Line Leak Detecting Method. A method of detecting leaks in pipe lines by dis- 
solving a small amount of radioactive substance in the material carried by the pipe and 
subsequently measuring the radioactive intensity along the pipe line. 


U.S. No. 2,346,093. W. A. Tolson. Iss. 4/4/44. App. 10/22/41. Assign. Radio Corp. of 
America. 
Sonic Depth Indicator. A sonic depth indicator having a slowly charging condenser 
discharged by the initial pulse and again discharged by the amplified reflected pulse, the 
depth being indicated by the current on discharge by the amplified reflected pulse. 


U.S. No. 2,346,486. D. G. C. Hare. Iss. 4/11/44. App. 10/7/41. Assign. The Texas Com- 
pany. 

Method and Apparatus for Measuring Thickness. A method of measuring thickness 
of a plate by placing a source of penetrating radiation on one side and observing the 
transmitted radiation at three random points on the other side and geometrically de- 
termining the transmitted intensity which would be observed at a point directly oppo- 
site the source. 


U.S. No. 2,346,655. H. Benioff. Iss. 4/18/44. App. 4/18/41. Assign. Submarine Signal 
Company. 


Electrodynamic Vibrator. A supersonic submarine signalling vibrator having a num- 
ber of electro-dynamic driving elements in each of which the moving coil is connected to 
the diaphragm through a half-wave length resonant vibrating system. 
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U.S. No. 2,347,702. H. B. Maris. Iss. 5/2/44. App. 4/11/41. 


Device for Measuring Extremely Small Angles. An optical device for measuring small 
angular motion between plane mirrors by measuring the displacement of a high order 
multiple reflection image. 


U.S. No. 2,347,759. T. C. Wherry. Iss. 5/2/44. App. 3/29/41. Assign. Phillips Petro- 
leum Co. 
Fishing Tool. A fishing tool which is supported on an electric conductor cable, the 
tool being equipped to indicate at the surface the proximity or contact with the fish 
either magnetically, electrically or by radioactivity. 


U.S. No. 2,348,520. W. T. Cardwell, Jr. Iss. 5/9/44. App. 5/5/42. Assign. Standard Oil 

Co. of California. 

Suspendometer. A device for comparing the relative capacities of fluids to support 
suspended material and consisting of an immersed frame with a shelf holding balls of 
graduated diameters, the frame being moved aside for a definite time and then returned 
to catch the balls on lower shelves. 


U.S. No. 2,348,810. D. G. C. Hare. Iss. 5/16/44. App. 5/29/40. Assign. The Texas Co. 


Method and Apparatus for Determining Liquid Level. A method of locating the liquid 
level through a closed container by moving vertically on the outside a source of gamma 
rays together with a detector of scattered radiation. 


U.S. No. 2,348,985. J. A. Lewis. Iss. 5/16/44. App. 4/19/40. Assign. Core Laboratories, 
Inc. 
Method for Determining Permeability. Measuring the permeability of rock samples 
by forcing air into a small circular area on a large face of the sample and observing the 
pressure and rate of flow obtained. 


U.S. No. 2,349,429. G. Herzog and J. H. Stein. Iss. 5/23/44. App. 4/29/42. 


Method and A pparatus for Measuring Thickness. A method of measuring thickness of 
pipe by passing penetrating radiation through the wall in a tangential direction and 
measuring the absorption. 


U.S. No. 2,350,072. E. E. Simmons, Jr. Iss. 5/30/44. App. 2/23/40 and 1/1/42. 


Torque Measuring Apparatus. A torsion dynamometer in which the torque trans- 
mitting member has tangential arms carrying bonded filament type electric strain 
gauges. 


U.S. No. 2,350,073. E. E. Simmons, Jr. Iss. 5/30/44. App. 5/8/42. 


Strain Gauge and Method for Making Same. A filament type electric strain gauge in 
which an insulated wire filament is embedded in a thermoplastic material which may be 
bonded to the test member by application of heat. 


U.S. No. 2,350,566. G. R. Mustaparta. Iss. 6/6/44. App. 3/25/41. Assign. Hercules 
Powder Co. 


Primer Tube for Dynamite Cartridges. A tube having a narrow inside pocket into 
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which a blasting cap may be placed so that the cap is not disturbed when the dynamite 
is inserted to make up the primer. 


U.S. No. 2,350,712. W. Barsties. Iss. 6/6/44. App. 3/25/41. Vested in Alien Property 

Custodian. 

Fluid Level Gauge. An immersion type tank fluid level gauge consisting of a long 
right angle prism having an extended light source along one base and an extended pho- 
tocell along the other base, the light being totally reflected at the hypotenuse to the 
photocell only where not immersed in liquid. 


U.S. No. 2,350,972. A. C. Ruge. Iss. 6/6/44. App. 9/16/39. 

Strain Gauge. A wire filament type strain gauge having two wires of different and 
opposite thermal coefficients of resistance so that temperature compensation may be ob- 
tained. 


U.S. No. 2,351,003. M. Camras and W. Korzon. Iss. 6/13/44. App. 9/11/40. Assign. 
Armour Research Foundation. 
Recording and Reproducing of Vibrations. A recording and play-back head for longi- 
tudinal magnetic recording on a steel wire which enters and leaves the recording gap 
through the pole pieces. 


U.S. No. 2,351,004 to 2,351,011. M. Carmas. Iss. 6/13/44. App. 12/22/41 to 2/25/43. 
‘Assign. Armour Research Foundation. 
Method and Meansof Magnetic Recording. A method of magnetic recording on a steel 
wire, elimination of disturbances, recording head designs, and control devices for mag- 
netic recorders. 


U.S. No. 2,351,201. R. Gillis. Iss. 6/13/44. App. 7/16/42. Assign. Western Electric Co. 


Method and Means for Measuring. A direct reading a-c four electrode method of 
measuring the thickness of a metallic sheet by passing a current between two points far 
apart and measuring the potential between two intermediate points close together. 


U.S. No. 2,351,353. M. D. McCarty. Iss. 6/13/44. App. 12/19/41. Assign. Socony- 
Vacuum Oil Co., Inc. 
Galvanometer. A galvanometer element having a double coil, one of which is used 
for recording the signal and the other either for exciting an ave circuit, mixing, inverse 
feedback or damping. 


U.S. No. 2,351,457. W. M. Rust, Jr. and M. A. Arthur. Iss. 6/13/44. App. sn/aulam 

Assign. Standard Oil Development Co. 

Grating Spectrograph. A spectrograph control unit which opens a shutter a predeter- 
mined time after the light source is started, closes the shutter a predetermined interval 
later and gives an audible signal which may be used to gauge photographic developing 
time. 


U.S. No. 2,351,508. H. G. Hamilton. Iss. 6/13/44. App. 5/8/43. Assign. The Singer 
Manufacturing Co. 


Electric Motor Control System. A speed control for small series motors consisting of 
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tandem rheostats, one of which is in series with the motor and the other of which con- 
trols the rectified current energizing a magnetic brake so that the brake is released as the 
motor is speeded up. 


Erratum 


In the July 1944 number, under Geochemical Prospecting, patent U.S. No. 2,239,- 
651 should read U.S. No. 2,339,651. 


PUBLICATIONS RECEIVED 


(The publications listed below were received since issuance of Groruysics for 
July, 1944, and are available for loan to the membership. Articles cited are those of 
particular interest to geophysicists.) 


American Journal of Science, Vol. 242, No. 6 (“Thorium-Uranium Ratios in Rocks and 
Minerals,” by N. B. Keevil, and “Radioactive Substances, III, The Analysis and 
Age of a North Carolina Monazite,”’ by Allen Douglass Bliss), and 7 (“The Distri- 
bution of Helium and Radioactivity in Rocks, VI, The Ayer Granite-Migmatite at 
Chelmsford, Mass.,” by N. B. Keevil, E. S. Larsen and F. J. Wank), and 8 (June, 
July, August, 1944), New Haven. 


Basic Principles of Weather Forecasting, by Victor P. Starr, Harper and Brothers, Pub- 
lishers, New York. 


Bulletin of the Academy of Sciences of the Union of Soviet Socialist Republics. No. 4 
(1943), and 6 (“The Study of Electric Current on Models,” by J. P. Bulashevitsh), 


1943, Moscow. 


Bulletin of the American Association of Petroleum Geologists, Vol. 28, Nos. 7, 8 (July, 
August, 1944), Tulsa. 


Economic Geology, Vol. XXXIX, No. 4 (June-July, 1944), Lancaster. 


Independent Petroleum Association of America Monthly, The, Vol. XV, Nos. 2, 3 and 4 
(June, July, August, 1944), Tulsa. 


Journal of Applied Physics, Vol. 15, Nos. 5, 6, 7 and 8 (May, June, July, August, 1944), 
Lancaster. 


Journal of the Institute of Petroleum, Vol. 30, Nos. 244, 245 (April, May, 1944), London. 
Petroleum, Vol. VII, Nos. 5, 6 and 7 (May, June, July, 1944), London. 


Review of Scientific Instruments, The, Vol. 15, Nos. 5, 6 (May, June, 1944), and 7 
(““Magnetometer for the Determination of the Vertical Component of the Earth’s 
Magnetic Field,” by S. L. Ting and S. T. Lin), July, 1944, Lancaster. 


Terrestrial Magnetism and Atmospheric Electricity, Vol. 49, No. 2 (June, 1944), Balti- 
more. 


World Petroleum, Vol. 15, No. 6 (“Oil Search in the Eastern Gulf Region,” by Joseph A. 
Kornfeld), 7 and 8 (June, July, August, 1944), New York. 
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gree in 1934, M.S. in 1936, and Ph.D. degree in 
Chemical Engineering in 1939, all from the Uni- 
versity of Texas. While doing graduate work at 
that school he was an instructor in chemistry 
for three years, and was a research associate in 
the Bureau of Industrial Chemistry for two 
years, working on problems associated with the 
decomposition of methane by high voltage dis- 
charge. 
After leaving school, he worked as field 
testing engineer for the Oil and Gas division of 
the Texas Railroad Commission. For a year and 
a half he served as professor of chemical engi- 
neering at Texas College of Arts and Industries. 
In December, 1940, he joined The Carter Oil 
Company as a research geophysicist in geo- 
chemical work. Mownro™ W. KRIEGEL 
He is a member of the American Chemical 
Society, The American Institute of Mining and Metallurgical Engineers, Phi Lambda 
Upsilon, Tau Beta Pi, anil associate of Sigma Xi. He holds a Texas professional engi- 
neer’s license in chemical engineering. 
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CONTRIBUTORS 


ERVAND GEORGE KOGBETLIANTz was born 
in southern Russia of Armenian parents. He 
graduated from the University of Moscow, and 
received the degree of Docteur es Sciences 
Mathematique at the Sorbonne in Paris. In 
1927 he founded, with the financial participa- 
tion of the Banque de Paris et de Pays Bas and 
the Compagnie Francaise des Petroles, the firm 
Société de Prospection Géophysique, which en- 
gaged in gravity and magnetic explorations. In 
1933 Dr. Kogbetliantz was sent by the govern- 
ment of France on a scientific mission to Iran, 
where he was appointed by the Iranian govern- 
ment as Professor of Applied Geophysics at the 
University of Teheran, a post he filled for five 
years. In 1938 he returned to France as a staff 
member of the National Center of Scientific 
Research, assigned to study and improve the E. G. KoGBETLIANTZ 
interpretation of geophysical maps. Dr. Kog- 
betliantz came to the United States in June of 1942 at the invitation of Lehigh Univer- 
sity, where he has been teaching mathematics and continuing with his researches. 


Tuomas S. WEsT attended the University 
of Oklahoma for four years during the period 
1922-28, and returned in 1930 for additional 
work in physics. During 1928-29 he was em- 
ployed as an engineer by the Tidal Refining 
Company. During the years 1931-41 he was 
associated with the McLaren, Grimes and Pat- 
terson interests of Dayton, Ohio, in electrical 
prospecting research and petroleum exploration 
and development. Mr. West has been in charge 
of geological and engineering work for the Hen- 
derson-Coquat interests of San Antonio, Texas, 
since 1941. 


Tuomas S. WEST 
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CONTRIBUTORS 


CLARENCE C. BEACHAM received the de- 
gree of B.S. in mining engineering from Ohio 
State University in 1928. He was employed by 
the Gypsy Oil Company as an engineer in the 
Seminole area and as district engineer in the 
Drumright district during 1928-32. During 
1931-41 he was associated with the McLaren, 
Grimes and Patterson interests of Dayton, 
Ohio, in electrical prospecting research and 
petroleum exploration and development. Mr. 
Beacham has been part owner and manager of 
the Harbough-Potts Coal Company of Somer- 
set, Ohio, since 1942. 


CLARENCE C. BEACHAM 


B. G. Swan’s photograph and biography appeared in Vol. VII, No. 4, p. 425 (Oc- 
tober, 1942). 
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THE SOCIETY ROUND TABLE | 
MEMBERSHIP APPLICATIONS APPROVED FOR PUBLICATION 


The Executive Committee has approved for publication the names of the following 
candidates for membership in the Society. This publication does not constitute an elec- 
tion, but places the names before the membership at large. If any member has informa- 
tion bearing on the qualifications of these nominees, he should send it to the Secretary 
within thirty days. (Names of sponsors are placed beneath the name of each nominee.) 


ACTIVE 
Charles Vidor Aderman 
C. C. Zimmerman, H. W. McDonnold, J. H. Deming 
Joe Bill Barbisch 


W. Harlan Taylor, P. E. Narvarte, W. Lee Moore 


Carl V. Benz 
W. Harlan Taylor, Leo Horvitz, L. W. Blau 


Victor Joseph Blum, S.J. 


James B. Macelwane, S.J., Stanley W. Wilcox, Cornelius G. Dahm 


Weldon L. Crawford 


W. Lee Moore, Jr., P. E. Narvarte, W. Harlan Taylor 


Lawrence Andrew Davis 

W. Harlan Taylor, R. B. Pratt, P. E. Narvarte 
Fred Joseph Di Giulio 

Henry Salvatori, Booth B. Strange, Robert Dyk 
Archibald Gilbert Eberts 

L. A. Scholl, Jr., Roy L. Lay, Geo. D. Mitchell, Jr. 
Jack Benjamin Ferguson 

A. L. Ladner, R. H. Tucker, W. Harlan Taylor 
Leon Fischer 

F. Goldstone, L. K. Mower, Robert L. Geyer 
Raymond Ellsworth Halsey 

W. Harlan Taylor, P. E. Narvarte, C, D. Whitsitt 
Carl Wilburn Hemmi 

O. S. Petty, W. Harlan Taylor, P. E. Narvarte 
Theron E. Hobbs 

Willis R. Dortch, H. R. Moorman, P. P. Conrad 
Judson George Jackson 

T. A. Manhart, A. J. Barthelmes, G. W. Westby 
Cornelius Dourve Keen 

Wm. M. Barret, L. I. Freeman, Ashton Crain 
Harold James Kidder : 

W. M. Tottenham, Paul E. Nash, H. C. Cortes 
Frank Russell Kittredge 

W. Harlan Taylor, O. S. Petty, R. B. Pratt 


Adelbert Kohler 
W. Harlan Taylor, P. E. Narvarte, R. B. Pratt 
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Ralph H. Lang 

D. L. Scott, J. B. Lovejoy, W. J. Osterhoudt 
Clair Edmond McClure 

Stanley W. Wilcox, Albert J. Barthelmes, Hugh M. Thralls 
Jack Voorhies McManus 

Charles M. Moore, Jr., Cecil H. Green, H. B. Peacock 
James Roy Maxey 


A. B. Bryan, G. E. Wagoner, R. W. Gemmer 
Roderick Mathew Nugent 

Henry C. Cortes, Paul E. Nash, J. C. Menefee 
James Earl Payne 

Henry C. Cortes, Paul E. Nash, J. C. Menefee 
Theodore Rozsa 

F. W. Oudt, D. S. Hughes, F. Goldstone 
Paul Truman Rumsey 

A. L. Ladner, R. H. Tucker, W. Harlan Taylor 
Donald Wilson St. Clair 

F. Goldstone, F. E. Deacon, L. K. Mower 
James LeRoy Sauls, Jr. 

L. A. Scholl, Jr., George D. Mitchell, Jr., Roy L. Lay 
Garland Newman Shell 

W. W. Newton, Hal R. Adams, Thos. R. Shugart 
E. Joe Shimek 

H. C. Cortes, P. E. Nash, J. P. Minton 
William Gordon Smith 

H. C. Maliphant, F. Goldstone, Donald R. Brown 
Roy Jackson Tribbey, Jr. ; 

H. R. Moorman, P. P. Conrad, D. R. Dobyns 
Merrill Douglas Tucker 

Joseph A. Sharpe, Sidon Harris, L. H. Wise 
John Averett Waller, Jr. 

W. W. Newton, T. R. Shugart, B. D. Farrow 
Donald Bean Winfrey 

W. Harlan Taylor, W. Lee Moore, R. B. Pratt 


ASSOCIATE 
Dorothy Earle Albertson 
L. A. Scholl, Jr., Roy L. Lay, C. R. Wallace 
Harry Raymond Conley 


L. A. Scholl, Jr., Roy L. Lay, C. R. Wallace 
Patrick Eugene Haggerty 

Cecil H. Green, Eugene McDermott, J. E. Jonsson 
John Blythe Halton Henderson 

Paul E. Nash, J. C. Menefee, W. C. Woolley 
John M. Kendall 

L. A. Scholl, Jr., Roy L. Lay, B. H. Treybig, Jr. 
William Priestley Leverich 

Tom D. Mayes, T. J. Bevan, Jack C. Pollard 
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James Harold McLay 
M. C. Kelsey, H. B. Peacock, Cecil H. Green 
William Edward Burckhardt Pappert 
Raymond A. Peterson, Chas. L. Heald, Barthold W. Sorge 


James Wheeler Vaikus 
M. C. Kelsey, E. J. Stulken, Earl Thomas 
Cal Yates 
R. S. Jackson, T. I. Harkins, D. F. Broussard 
STUDENT 
Paul Clarence Wuenschel 
R. Maurice Tripp 


PRELIMINARY ANNOUNCEMENT OF FIFTEENTH ANNUAL MEETING 
OF THE SOCIETY 


The S.E.G. will meet concurrently with the A.A.P.G. and S.E.P.M. in Tulsa, 
March 20-22. Mr. Henry C. Cortes, Vice-President of the S.E.G., is General Chairman 
of the Program and Arrangements Committee. He will be assisted by a host city ar- 
‘rangements committee, a host city program committee, and program committee mem- 
bers from various regions of the country. 

The chairman of the host city arrangements committee is T. A. Manhart, Seismo- 
graph Service Corporation, 712 Kennedy Building, Tulsa 3. The chairman of the host 
city program committee is L. Y. Faust, Geophysical Research Corporation, Box 2040, 
Tulsa 2. Other host city committee members and regional members of the program com- 
mittee will be listed in the January number of GEopHysics and in the complete an- 
nouncement of the meeting to be circulated to all members. 

Members of the Society should begin to think about the preparation of papers for 
presentation in person or by title at the annual meeting, and for subsequent publication 
in Gropuysics. If there is a possibility that you may be able to prepare a paper, or if 
you know another member who is likely to be able to prepare a paper, kindly communi- 
cate this information to the General Chairman of the Program and Arrangements 
Committee. 
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PERSONAL ITEMS 


Bitty E. Ricwarps, formerly with the United Geophysical Company, is now serv- 
ing as a Lieutenant in the U. S. Army Air Forces. He may be addressed at A.P.O. 713, 
Unit 1, c/o Postmaster, San Francisco, Calif. 


J. Carr is being transferred from the Shell Co. of Australia to the Shell Co. of 
Egypt Ltd., Cairo, Egypt. 

OvELL C. Otson is a Captain in the U. S. Marine Corps Reserves. His mailing ad- 
dress is P. O. Box 349, Brady, Texas. 


CuHaRLEs C. WILLiAMs, formerly with the National Geophysical Co., is now asso- 
ciated with the General Geophysical Co., 2513-14 Gulf Bldg., Houston 2, Tex. 


Jack F. Jupson, of the United Geophysical Co., has been trasnferred from San 
Felipe, Yaracuy, Venezuela, to Caracas. He may be addressed in care of the company at 
Aptdo. 1085. 


H. G. Dott is Chairman of the Board of Directors and Director of Research of the 
Schlumberger Well Surveying Corp., 2720 Leeland, Houston, Tex. 


Lt. Ir. J. J: H. JANSEN M., who has been serving with the Dutch Air Forces, may 
now be reached through the Netherlands Military Mission, 2010 Massachusetts Ave., 
NW., Washington, D.C. 


G. F. KaurMann, who has been serving overseas with the Standard-Vacuum Oil 
Co., is now located in New York for an indefinite period and may be addressed at Room 
1421, 26 Broadway, in care of that company. 


Pau FarrEn is now residing at 807 Sixth St., Ballinger, Tex. 


WHEATON C. CrarK has left the Continental Oil Co. to accept a position as District 
Geologist with the British-American Oil Producing Co., 312 Con Roy Bldg., Casper, 
Wyo. 


VINCENT J. MERCIER, of the Lane-Wells Co., has notified the Society of his transfer 
to Wichita 2, Kansas, from Pratt. The Lane-Wells office is at 806 Central Bldg. 


RIcHARD WILLIAMs, Supervisor with the National Geophysical CoZ is located i in 
Palestine, Texas, with the postal address P. O. Box 487. 


A. A. Hunzicxer, of the Texas Co., has been transferred from Nashville, Ill., to 
Colby, Kansas. 


Lr. S. C. SrONEHAM may be addressed at P. O. Box 375, Waxahachie, Tex. 


S. H. Wrttston, Vice President of Horse Heaven Mines, Inc., and the Cordero 
Mining Co., has moved from Portland, Oregon, to 703 Mechanics’ Institute Bldg., 57 
Post St., San Francisco 4, Calif. 


H. C. Bicxet, Chief Geophysicist with the Tropical Oil Co., is located at Bogota, 
Colombia, with the following addresses: Aptdo. Aereo 3533 and Aptdo. Nacional 335. 


R. Maurice Tripp has left the Colorado School of Mines to serve as Assistant to 
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PERSONAL ITEMS . 579 
the President of the Geotechnical Corp. His address is 149 Sidney St., Cambridge 39, 
Mass. 


Lewis H. Boyp, of the Brown Geophysical Co., has moved from Bakersfield, 
Calif., to Tallahassee, Fla., where his postal address is Box 724. 


J. W. Tuomas, of Geophysical Service, Inc., is supervising work with the Depart- 
ment of Exploration, Petroleos Mexicanos, Av. Juarez 95, Mexico, D. F. 


V. L. JONnEs is now Physicist with the Stanolind Oil and Gas Co., 1321 S. Carson, 
Tulsa, Okla. 


LAWRENCE A. GOEBEL, of the Carter Oil Co., may be addressed at P. O. Box 801, 
Tulsa, Okla. 


Davip M. WEBER has returned from foreign duty with the Mene Grande Oil Co. 
and may be addressed in care of the Gulf Res. and Dev. Co., Box 352, Laurel, Miss. 


GerorcE B. Stone, after serving in the Naval Ordnance Laboratory for the past 
year, has joined the Phillips Petroleum Co., 521 Esperson Bldg., Houston 2, Texas, 
as Geologist and Geophysicist. His residence address is 2234 Tangley Rd. in that city. 


LEONARD J. NEUMAN has offices at 943 Mellie Esperson Bldg., Houston 2, Tex. 


C. Russett Hate, Party Chief with the United Geophysical Co., and recently 
transferred from Brazil to Venezuela, is having his mail forwarded to him from his home 
address, 502 Meander St., Abilene, Tex. 


Paut H. Boots, of the Gulf Res. and Dev. Co., is now serving with the Colombian 
Gulf Oil Co., Aptdo. Aereo 40-14 or Aptdo. Nacional 1339, Bogota, Colombia, S. A. 


HARDEE CHAMBLISS, retired Professor of Chemistry, and sometimes Dean of the 
School of Engineering, Catholic University of America, has joined the staff of the Geo- 
physical Instrument Co., Arlington, Va., as Consulting Chemist. In the Graduate School 
of the U. S. Dept. of Agriculture, Dr. Chambliss will give a general, introductory course 
in geochemistry. Geochemical prospecting for oil and minerals will be among the topics 
discussed. This is probably the first course of the kind to be offered in the U. S. 


A. L. WILKinson, of the Petty Geophysical Engineering Co., has been transferred ; 
from Homestead, Fla., to Sour Lake, Tex. His postal address is Box 415. 


Captain Lynn A. SMITHERMAN, U. S. Army Signal Corps., should be addressed at 
P.O. Box 67, Athens, Tex. 


Lr. Cot. DREXLER DANA should be addressed at Box 949, Route 7, Bakersfield, 
Calif. 


The mailing address for F. W. Oupt and T. Rozsa, of the Shell Oil Co., Inc., has 
been changed from Box 1007 to Box 1049, Saginaw, Mich. 


Lr. D. G. Ecan, U. S. Army Signal Corps, may be addressed at P. O. Box 366, 
Crowley, La. 


. H. G. Mace, of the United Geophysical Co., recently elected to associate mem- 


| 
| 


580 PERSONAL ITEMS 


bership in the Society, is working in the mid-continent area. He may be addressed at 
804 Thompson Bldg., Tulsa 3, Okla. 

O. A. SEAGER has joined the staff of the Standard Oil Co. of Egypt, 22 Sharia Kasr 
el Nil, Cairo, Egypt. 

Joun A. LEsTER, Seismic Supervisor with the Magnolia Petroleum Co., has been 
transferred from Dallas to Houston, Texas, with the postal address Box 111. 

Don B. HEANEY, formerly Computer on Western Geophysical Co. Party 20, is now 
serving as Radio Technician Second Class in the U. S. Navy. He may be reached 
through the Fleet Post Office, San Francisco, Calif. 

LELAND W. Jones, of the Anderson-Prichard Oil Corp., advises the Society of the 
change of offices from Oklahoma City to 609 Orpheum Bldg., Wichita 2, Kansas. 

Captain MERLE C. Bowsky, of the U. S. Army Signal Corps, may be addressed at 
A.P.O. 5541, c/o Postmaster, New York, N. Y. 

Wa ter D. Barrp, Party Chief for Geophysical Service, Inc., has been transferred 
to the Mid-Continent area, and may be addressed at 1311 Republic Bank Bldg., Dallas, 

- Texas. 

James A. Lone, Party Chief with the Stanolind Oil and Gas Co., should be ad- 
dressed in care of the company at Box 591, Tulsa, Okla. 

R. G. Sontserc is Geophysicist with South Mediterranean Oilfields, Ltd., 43 
Sharia Kasr el-Nil, Cairo, Egypt. 

RoseErt G. NISLE has been advanced to the rank of Lt. (j.g.), and is serving as 
Project Officer, Bureau of Ordnance, U. S. Navy Dept. His residence address in Wash- 
ington, D. C., is 332 Raleigh St., SE. 

J. B. Souruer, Consulting Geologist and Geophysicist, has returned to San An- 
tonio, Texas, from Lincoln, Nebraska. He may again be addressed at 239 Luther Drive, 
San Antonio 1. y 

G. W. Carr announces the establishment of the Carr Geophysical Co., with offices 
at 1321 Commerce Bldg., Houston 2, Tex. 

G. W. Gutmon, Geologist with the California Co., may be reached in Ferriday, La., 
at Box 546. 

Lt. R. G. Paterson, R. C. N. V. R., is Electrical Anti-Mining Officer with the 
Royal Canadian Navy. His mailing address is c/o Fleet Mail Office, Vancouver, B. C., 
Canada. 


C. H. Acueson has returned to Canada from South America. He is located at 301 
_ Examiner Bldg., Calgary, Alberta, Canada. 


Roy LEE GALttEway, of the Texas Co., has left Lafayette, La., for Plaquemine, 
where his posta] address is Box 150. 


The PETROLEUM ADMINISTRATION FOR Wak reported in September that there were 
481 geophysical units operating in the United States, as follows: seismic, 290 units; 
gravimeter, 124 units; core drills, 47 units; magnetic, 17 units; torsion balance, 1 unit; 
electrical, 2 units. The number of units operating increased by 72 since January 1, 1944. 
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INDEX TO VOLUME IX OF GEOPHYSICS FOR 1944 


No. Page 


Agocs, W. B. A Method of Determining the Time Break on tne Sea Seis- 
mic Records from the Water Sound Arrivals . 
. Contributors. . 
Alien Property Custodian, Patents Vested in the.O.F.Ritzmann . 
—— ~ Hydrocarbons in the Presence of Nitrous Oxide. Monroe W. 
Kriege 
Analytical and Experimental Data Concerning : a Solid Hydrocarbon Survey 
of the Fort Collins Anticline. R. Maurice Tripp . 
— Meeting, Minutes of the. Hart Brown. The Society Round 
Annual Meeting, 1944, ’ Announcements of. The Society Round Table . 
, Fourteenth, of the Society. The Society Round Table . . 
; the F ourteenth, Members Registered at. The Society Round Table. 
; the Fourteenth, Papers Presented or Read by Title at the Technical 
Sessions. The Society Round Table . 
—_e of the Secretary-Treasurer. Hart Brown. The Society Round 
able. . 
Arkansas Bauxite R ion, ’ Geophysical Survey of the. Mark C. Malamphy 
and James L. Vallely. 
Association of Magnetic and Density Contrasts with Igneous Rock C Clas- 
sifications. L. L. Nettleton and T. A. Elkins 


> 


Ballot, 1944, Official. The Society Round Table . . 
on Editor’s Term of Office. The Society Round Pie 38 3 423 
on Editor’s Term of Office, Results of. The Society Round Table. 4 577 


Barnes, Virgil E., and Frederick Romberg. Correlation of Gravity Observa- 
tions with the Geology of the Smoothingiron Granite Mass, Llano 
County, Texas I 79 

Bauxite Region, the Arkansas, Geophysical Survey of. Mark a Malamphy 
and James L. Vallely . . ‘ 3 

Beacham, Clarence C. Contributors. 

: and Thomas S. West. Measurement of Deep Electrical Anomalies . 4 

By-Laws, Constitution and, of the Society. . . pis 


Cayuga Field, Anderson County, Texas, The Geophysical History of er 
H. B. Peacock ‘ 3: 290 

Charge Size, The Effect of, on Reflection Records. Joseph A. Sharpe 2 

Communications, Discussions and. Appointment of Special Velocity Survey 


Editor. Joseph A. Sharpe . 3 408 
. “Education and the Industry. Professor V. C. Iling. 
Reprint of section titled “Exploration” . 2 217 ‘ 


. Letter, D. R. Knowlton, Director of Production, Petroleum Ad- 
ministration for War, to W. H. Taylor, Secretary-Treasurer, Society of 
Exploration Geophysicists, relative to “Classification of 
Exploratory Drilling for Petroleum” . . 3 409 

Constitution and By-Laws of the Society . . 2 243 

Contributors. Dart Wantland, R. W. Lawrence, T. T.R. - Shugart, Stanley W. 


Wilcox, A. E. Lockenvitz . . I 
. William L. Russell and W. B. Agocs . 2 230 
. Mark C. Malamphy, James L. Vallely, R. Maurice Tripp . : ¥ - 43 
. Clarence C. Beacham, Ervand Kogbetliantz, Monroe W. Kriegel, 
John Sloat, Thomas S. West . 4 §72 


Correlation of os ao County, Te with the Geology of the Smoothingiron 
ames . 79 
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Data, Analytical and Experimental, Concerning a Solid Hydrocarbon Sur- 
vey of the Fort Collins Anticline. R. Maurice Tripp 

Deep Electrical Anomalies, Measurement of. Thomas S. West and Clarence 
C. Beacham . 

Sea Seismic Records from the Water Sound Arrivals, A Method of 
Determining the Time Break on. W. B. Agocs 

Density Contrasts, Association of Magnetic and Gravity, with Igneous 
Rock Classfications. L. L. Nettleton and T. A. Elkins . 

Detonation in Explosives, Mechanism of. Robert W. Lawrence. 

Discussions and Communications. Appointment of Special Velocity Survey 
Editor. Joseph A. Sharpe. 

_ “Education and the Petroleum Industry. » Professor V. C. Illing. 

Reprint of section titled “Exploration.” . 

. Letter; D. R. Knowlton, Director of Production, Petroleum Ad- 

ministration for War, to W. H. Taylor, Secretary-Treasurer, Society 

of Exploration Geophysicists, relative to er titled, “Classification 

of Exploratory Drilling for Petroleum.” ; 


Earth a Method, Sand and Gravel Prospecting by the. Stanley W. 
Wilcox. 

Effect of Charge Size on Reflection Records, The. "Joseph A. Sharpe 

Elastic Earth, The Equation of Motion of a Geophone on the Surface of an. 
Alfred Wolf 

Electrical Anomalies, Deep, Measurement of. Thomas S. West and Clarence 
C. Beacham . 

Elevation Factor, A Proposed Least ‘Square Method for the Determination 
of the. J. A. Legge, Jr. . 

Elkins, T. A., and L. L. Nettleton. Association of Magnetic and Density 
Contrasts with Igneous Rock Classifications. . 

Equation of Motion of a Geophone on the Surface of an Elastic Earth, The. 
Alfred Wolf 

Executive Committee, Nominees for 1944-45. W. M. “Rust, ‘Jr., ean J 
Bazzoni, Roland F. Beers, Henry C. Cortes, W. Harlan Taylor, I Dart 
Wantland. The Society Round Table . : 

Explosives, Mechanism of Detonation in. Robert W. Lawrence. 


Fort Collins Anticline, Analytical and Experimental Data Concerning a 
Solid Hydrocarbon Survey of the. R. Maurice 7a : 
Frequency Discrimination in the Low Velocity Zone. T. R. Shugart 


Gallie, J. F. History of the Society of Exploration Geophysicists. : 

Gamma Ray Activity, The Total, of Sedimentary Rocks as Indicated 1 by 
Geiger Counter Determinations. W. L. Russell. 

Garrett, Lovic Pierce. Memorial by Paul Weaver. 

Geiger Counter Determinations, The Total Gamma Ray Activity of Sedi- 
mentary Rocks as Indicated by. W. L. Russell . . 

Geology of the Smoothingiron Granite Mass, Llano County, Texas, Corre- 
oe of Gravity Observations with the. Frederick Romberg and Virgil 

Barnes. . 

Geophone on the Surface of an Elastic Earth, ‘The Equation of Motion of 
a. Alfred Wolf. . 

Geophysical Exploration i in Poland. S. M. Wir 

History of the Cayuga Field, Anderson County, Texas, The. H. B. 

Peacock 

Survey of the Arkansas Bauxite Region. Mark C. >. Malamphy ‘and 
James L. Vallely (es 


Geophysics Looks R. D. Wyckoff 
Gravel — Sand and s by the Earth Resistivity Method. Stanley W. 
Wilcox . 


131 


494 


238 


180 


180 


79 


29 
143 


299 


324 
287 


36 


3 367 
4 404 
163 
I 60 
— I I 
2 217 
409 
| 
2 
I 29 
4 
7 2 175 
60 
I 122 
2 
: 3 3 
2 
a > 


INDEX TO VOLUME IX 


Gravitational Force, The Periodic Variations of the. I. Arthur E. Lockenvitz 
Gravity Maps, Magnetic and, Quantitative Interpretation ‘of. Ervand 
Kogbetliantz . 4 447 
Observations, Correlation of, ‘with the Geology of the Smoothingiron 
Llano Texas. Frederick and E. 
arnes 


Harris, Sidon. Memorial to Lester Harold Wise (1905-1943) . . . . 2 
History of the Society of Exploration Geophysicists. J. F.Gallie . ~. . 2 238 
, The Geophysical, of the Cayuga Field, Anderson os Texas. 

H. B. Peacock 3 299 
How Not to Find an Oil Field. John 
Horton, C. W. A Note on the Reading of Seismograms : 2 160 
Hydrocarbon Survey, Solid, of the Fort Collins Anticline, Analytical ‘and 

Experimental Data Concerning a. R. Maurice Tripp 2°... 367 i 


Hydrocarbons, Analysis for, in the Presence of Nitrous Oxide. Monroe W. 
Kriegel. 


Igneous Rock Classifications, Association of Magnetic and Density Con- 
trasts with. L. L. Nettleton and T. A. Elkins .* . I 

Index of Wells Shot for veneer | B. G. Swan, Special Velocity Survey 

Editor . 


Kearby, John Gallatin. Memorial 2: 
Kogbetliantz, Ervand. Contributors. 
. Quantitative Interpretation of Magnetic and ‘Gravity Maps 4° 463 
Kriegel, Monroe W. —— for asain in the Presence of Nitrous 
. Contributors . 4 


Lawrence, Robert W. Contributors. . . 117 


. Mechanism of Detonation in Explosives as I I 
Least Square Method, A Proposed, for the Determination of the Elevation 
Factor. J. A. Legg é, Jr. . 2. 


Legge, J. A., Jr. A Proposed Least Square M Method for the Determination of 


the Elevation Factor. a. 
Lockenvitz, Arthur E. Contributors. : 
. The Periodic Variations of the Gravitational Force.I . I 04 
Low Velocity Zone, Frequency Discrimination in the. T. R. Shugart. I 19 
Malamphy, Mark C. Contributors 
, and James L. Vallely. Geophysical Survey of the Arkansas Bauxite 
Region 
Magnetic and Density Contrasts, Association of, with Igneous Rock Classi- 
fications. L. L. Nettleton and T. A. Elkins .. I 60 
and Gravity Maps, of. Ervand Kogbet- 
liantz : 4 463 
Interpretation. “Dart Wantland I 47 
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Shock mounted in a 6x 6x4 inch hard- each temperature compensated for precise 
wood case and weighing only 214 pounds, _ observations. 


the instrument incorporates the portability Write for illustrated technical bulletins. 
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ibrary.’ 
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154 pp., 83 illus., 6 x 9 inches 
Repaged from the August, 1941, A.A.P.G. Bulletin, in paper cover 
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The following reprints from early publications of the Society of Exploration Geophysi- 
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of print), which are offered in a group at $1.00. As the supply of some articles is very 
limited, those interested should place their orders promptly. 


THE JOURNAL OF THE SOCIETY OF PETROLEUM GEOPHYSICISTS 
Volume VI, Number 1 July, 1935 
Explosives and Electric Blasting Caps for Geophysical Prospecting 

: G. H. Loving and G. H. Smith 
Note on the Theory of Seismic Prospecting C. H. Dis 
Notes on the Early History of Applied Geophysics in the Petroleum Industry E. DeGolyer 
Portable Dynamite Storage Magazines J. W. Flude 


GEOPHYSICS 


Volume I, Number 1 January, 1936 
(Reprints listed not sold singly; must be purchased as unit at $1.) 

Magnetic Prospecting in Santa Catharina, Brazil M. C. Malamphy 

The Dip Needle as a Magnetometer C. O. Swanson 

*A New Reflection System with Controlled Directional Sensitivity F, Rieber 

Fundamental Photographic Processing Operations Influencing Production of Seismo- 
graph Records F. A. Tompkins 

*The Place of the Vertical Gradient in Gravitational Interpretations H. M. Evjen 


Summary of Some Pending Patent Litigation Relating to Seismic Exploration 
C. R. Hrdlicka 


Exploration and Production in the Gulf Coast through 1935 (charts only; article 
out of print) E. E. Rosaire & M. E. Stiles 


The Geophysics of the Tomball Oil Field Harris County, Texas- J. Brian Eby 
Regional Geophysical Activity Reports for 1935 J. H. Wilson, E. McDermott, 
E. L. Caster, H. G. Patrick, J. B. Eby, J. P. Schumacher, E, E. Rosaire, K. Ransone 


GEOPHYSICS 
Volume I, Number 2 ; July, 1936 


The Organization of An Effective Exploration Department B. B. Weatherby | 


A Proposed Geophysical Program of Exploration for Nebraska and the Dakotas 
J. H. Wilson 


*Explosives for Seismic Prospecting N. G. Johnson & G. H, Smith 
*Note on Reflections from Steeply Dipping Beds C. A. Heiland 
*The Amplitudes of Waves to be Expected in Seismic Prospecting B. Gutenberg 
*Relation between Firing Current and Performance in Seismograpli Caps JL. A. Burrows 


GEOPHYSICS 
Volume I, Number 3 October, 1936 


Some Aspects of Multiple Recording in Seismic Prospecting P. W. Klipsch 
GEOPHYSICS 
Volume III, Number 4 : October, 1938 


Three Dimensional Reflection Control S. M. Rock 
Nomogram for Dip Computations Reed Lawlor 


* Supply very limited. 
Send orders and accompanying remittance to the 


SOCIETY OF EXPLORATION GEOPHYSICISTS 
P.O. Box 410 El Dorado, Arkansas - 


Please mention GropHysics when answering advertisers 


— 
a 
3 
| 
j 
4 
ig 
‘ 
4 
: 
i 


26 


Barret Magnetic Surveys 


Where to Look and Where Not 
to Look for Oil | 


Barret Magnetic Surveys, because dependable and com- 
paratively inexpensive, are widely used for the pre- 
liminary investigation of large areas. They separate the 
wheat from the chaff. They define acreag~ barren of 
promise and disclose areas which justify more search- 
ing geophysical investigation. 


Because of marked economy, speed of operation, con- 
sistent dependability, and disclosures of untold value 
to clients, Barret Magnetic Surveys have grown steadily 
in favor for almost 17 years—and with a most impres- 


sive clientele. 


Please mention GropHysics when answering advertisers 


f the Societ 3 
y of Explorati 
on Geophysicists 
4 
M. BARRET: {NC- 
CONSULTING 
| LA. 


PEED PORTABLE ORILLING UNIT: 


HOUSTON, 
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THOMPSON HARDWARE COMPANY 
HOBBS, NEW MEXICO 


PHONES 166, 167, 19 
Carload Explosives Distributor 


WHEN IN WEST TEXAS OR NEW MEXICO, CALL US FOR 
YOUR EXPLOSIVES REQUIREMENTS. 
ANY BRAND, ANY TIME, ANY PLACE. 


GEOPHYSICS 


Volume VI OCTOBER, 1941 Number 4 


A Direct Reading Phase Shift Meter Joseph D. Eisler 
The Relation Between Depth, Lithology, and Seismic Wave Velocity in Tertiary 
Sandstones and Shales N. A. Haskell 
Seismic. Velocities in the. Southeastern San Joaquin Valley of California £E. J. Stulken 
The Probable Errors of Delta-T Velocities Frederick Romberg 
Comparison of Well Survey and Reflection ‘“‘Time-Delta Time’’ Velocities 
; W. E. Steele, Jr. 
Notes on Refraction Prospecting Charles Hewitt Dix 
Current Penetration in Direct Current Prospecting M. Muskat and H,. H. Evinger 
Prospecting Effectiveness Esme Eugene Rosaire 


Considerations on the Vertical Migration of Gases Robert G. Nisle 


This issue is $1 to members, $2 to non-members 
' (add 20¢ for foreign postage). 


SOCIETY OF EXPLORATION GEOPHYSICISTS 
P. O. Box 410 — El Dorado, Arkansas 
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Now you can correct driller’s logs— 
core logs and old electrical logs to check 
depths and thicknesses of all formations 


RADIOACTIVITY WELL LOGGING 


Many fields were developed before the days 
of electrical logging. In these wells the geo- 
logical information has been derived chiefly 
from driller’s records, and sample logs of 
dubious reliability. Potential pay horizons, 
either oil or gas, have been cased-off. In 
many of these wells electrical logs have 
been made of deepened portions of the bore- 
hole. «Radioactivity Well Logging 
provides theinformation necessary 
to plan and execute a successful 


Write for Bulletins 


remedial program by locating accurately the 
depths and thicknesses of potential zones for 
subsequent gun perforating. The curves can 
be correlated and checked against old rec- 
ords and electrical logs made of wells drilled 
later in the same field. « Lane-Wells Radio- 
activity Well Logging Service is designed 
to help you get the information you need to 
obtain the maximum ultimate pro- 
duction from each well. « Call 
your nearest Lane-Wells branch. 


Proncers tn Of Elcctrontes 
LOS ANGELES HOUSTON OKLAHOMA CITY 


General Offices, Export Office and Plant: 
5610 So. Soto Street, Los Angeles 11, California 


24-HOUR SERVICE «30 BRANCHES 


Radioactivity Well Logging licensed by 
Well Surveys, Inc. * Tulsa, Oklahoma 


Please mention GeopHysics when answering advertisers 
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SEISMIC REFLECTION 


EQUIPMENT, reconditioned 


SEVEN YEARS 


Noted for its ability to stand up under 
rugged field conditions, Heiland equip- 
ment maintains its high standard of 
performance after long years of use. 


RESEARCH CORPORATION 
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BAROID 


WELL LOGGING 


content while drilling 


BAROID WELL LOGGING SERVICE provides a successsful method 
of determining the presence or absence of oil or gas in the vari- 
ous formations encountered while drilling, and of logging the 
depths and approximate thicknesses of such formations. The 
method is simple and also gives the operator additional infor- 
mation which heretofore has never been available and which 
enables him to drill a well more intelligently and with greater 
freedom from drilling difficulties. 

BAROID WELL LOGGING SERVICE has minimized coring in many 
fields and has been highly successful in the drilling of wildcat 
and exploratory wells. 

Units, with trained operators, are available in oil fields through- 
out the United States. 


BAROID SALES DIVISION 


NATIONAL LEAD COMPANY 


BAROID SALES OFFICES: LOS ANGELES 12 « TULSA 3 * HOUSTON 2 


Please mention GropHysics when answering advertisers 
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You Get This 
EXTRA SAFETY 
at No Extra Cost! 


Faced with labor shortages, many operators—even though they have 
already adopted many safety measures—are keeping a weather-eye 
open to any unusual and additional methods that offer protection 
for their men. 


At no extra cost, Atlas Manasite detonators—because they are less 
sensitive to impact and friction—put extra safety on the blasting 
Seb: They provide reliable action, too—millions have been sold since : 
Atlas Manasite detonators were introduced five years ago. 


Put this extra safety to work on YOUR blasting job. Call the Atlas 
Representative—he'll gladly give you full details. 


MANASITE—Reg. U. S. Pat. Off. 


ATLAS MANASITE DETONATORS 
ATLAS 


Powder Company 
Wilmington 99, Delaware 


Please mention Gzopxysics when answering advertisers 
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|Seismograph Service Corporation of Delaware 
CONSULTING EXPLORATION GEOPHYSICISTS 
‘races, Venezuela BUILDING  Bogot, Colombia 
TULSA, OKLAHOMA, U.S.A. 
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Whatever your exploration problems may be, the 
world-wide experience of G.S.1. is available to help 
you solve them. We invite you to consult us on 
your present and post-war programs. Call in G.S.1. , 


EOPHYSICAL SERVICE INC. 


McDERMOTT, President 


BRANCH OFFICE: HOUSTON, TEXAS | 
DALLAS. TEXAS 
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